PRODUCING ALL-INERTIAL GUIDANCE 


Many successful full-range test flights have proved the reliability 
and accuracy of the Arma inertial guidance system now in full 
production. Other advantages—salvo firing, immunity to jamming, 
low cost, a minimum of ground support equipment—are inherent 
in inertial guidance. Currently at Arma, company-funded re- 
search programs are studying smaller, super-sensitive devices 
for navigation and satellite instrumentation. 


ARMA, Garden City, New York, a division of American Bosch 
Arma Corporation . . . the future is our business. 
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DME « TACAN VorTAC 
ITT’S ANSWER To ANY 
OF TODAY’S TRAFF 
CONTROL PROBLEMS 
Just five years after the famous, 
Wright brothers’ Kitty Hawk 
flight, ITT started in avionics, 
with Such a heritage ig 
air Gaigation it is no wonder 
IFT developed DME 
distance Measuring equipment 
Which provides pilots with 
| “Gistance acCuracy of pilus or 
 miius two-tenths of a mile, 
yet Weighs only 35 pounds. 


DME aboard today’s aircraft 
tunes to any DME, VORTAC 
or TACAN ground station : 
4or-continuous distance information, 
The evelopment of VORTAC 
TACAN, the Civilian and 
“military rho-theta systems, is an 
_ TET achievement typical of 
its great capabilities in avionics, 
Now even greater 
is attainable through the . 
_ fusing of two divisions to form 
HTT Federal Laboratories. in a 
“single company are Research and 
_ Development plus Manufacture, 
Maintenance and Service... 
_teady to serve the military and 
industry with the shortest 
2 possible cycle between jnitia! 
_ Concept and delivered system. 


EDERAL LABOR 
WASHINGTON AVENUE NUTLEY, NEW 
PIVISION OF INTERNATIONAL TELEPHONE AND TELEGRAPH CORPORATION 


_ SPACE AND MISSILE 
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CAREER OPPORTUNITIES AT JPL 
IN THESE FIELDS — NOW 


Electronic Engineers 

. for component and system design of deep 
space communications, instrumentation, and auto- 
matic control equipments. 

. for microwave and RF solid state circuit design 
and flight evaluation. 

. for project management assignment on ad- 
vanced development and contracted effort in space 
communications. 

Physicists 

. for analysis in communications theory, orbital 
mechanics, guidance and control, and systems 
performance. 

. for analysis of digital communication and con- 
trol systems; real-time digital computer and 
closed-loop systems. 

. for research and development of servo and 
control mechanisms for large ground based and 
spacecraft antenna systems. 


Other opportunities exist for electronic engineers 
and physicists in many areas at JPL which has 
been assigned the responsibility for the nation’s 
Lunar, Planetary and Interplanetary unmanned 
exploration programs. 
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MOON BOUNCE...a collaborative project of the 
National Aeronautics and Space Administration, 
the Jet Propulsion Laboratory, and the Australian 
Ministry of Supply to link two continents by radio 
signals bounced off the Moon 


ELAPSED TIME 
2.44 ‘SECONDS 


On February 10, 1961, California and Australia 
were linked in the first international space com- 
munication experiment that bounced voice mes- 
sages between the two points via the Moon. The 
words were beamed at the Moon from the Jet 
Propulsion Laboratory transmitter at Goldstone, 
California to the receiver at Woomera, Australia. 


Principals in the conversation were Dr. Hugh L. 
Dryden, NASA Deputy Director, whose voice was 
relayed from Washington by telephone; Dr. Lee 
DuBridge, President of California Institute of Tech- 
nology, who spoke directly from Goldstone; and Alan 
Hulme, Australian Minister of Supply at Woomera. 


The occasion tested the new Australian station, 
the second of three Deep Space Instrumentation 
stations developed and directed for the National 
Aeronautics and Space Administration by the Jet 
Propulsion Laboratory. 


CALIFORNIA INSTITUTE OF TECHNOLOGY 


JET PROPULSION LABORATORY 
PASADENA, CALIFORNIA 
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Five layers 


A number of new oxide compounds with the perovskite 
structure have been developed by one of our research 
groups. Investigations of related three-layered and five- 
layered structures are also in progress. To discover the 
relationships between crystal structure and physi- 
cal properties, intensive studies are being carried on to 
determine ferroelectric, ferromagnetic and dielectric 
characteristics. From this basic research may come ad- 
vanced materials superior to any currently available. 
® Our corporate-sponsored research programs into the 
fundamental nature of matter offer unusual opportunity 
to any scientist interested in doing original work. In- 
quiries are invited particularly from those with experi- 
ence in— Solid State Physics * Plasma Physics * Particle 
Physics * Nuclear Engineering * Gaseous Electronics * 
Direct Conversion * Surface Chemistry. 


Please write to MR. W. G. WALSH, Personnel Department: 


RESEARCH LABORATORIES 


UNITED AIRCRAFT CORPORATION 
400 Main Street, East Hartford 8, Conn. 


April 1961 + Aerospace Engineering 


| 
3 
q 
AF 
Perovskite 
SN 
— 
i 3 = 
4 


= 
: 
se 

poses 
5525 


of engine parts reduces produc- 
tion lead time at Rocketdyne. 


Computer-listed master catalog 
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i Behind the thundering performance of Rocketdyne’s en- 


statetetstet' gines, a significant reduction in the cost of power for Amer- 

stores ica’s missiles has been quietly achieved. Today, Rocketdyne 
engineering skill and efficient production methods make it 
possible to power two missiles for the cost of one in 1957. 

Rocketdyne, the pioneer in rocket science, was first with 
power for America’s long range ballistic missiles —first with 
power for outer space. In establishing this technological lead- 
ership, Rocketdyne developed new management concepts at 
every level of operation, from early design through final test- 
ing. The result is outstanding technical achievement at the 
lowest possible cost. 

In data processing alone, advanced techniques are saving 
engineers hundreds of hours of experimentation and testing 
and have contributed to a 37 percent reduction in Atlas 
engine costs for the Air Force. An intracompany communica- 
tions network links test stands and research laboratories in 
Missouri, Texas and California; gives management the daily 
status of every program—whether it’s on schedule, what parts 
are in short supply, how the production line is performing. 

Through research, engineering, and management, 
Rocketdyne is constantly at work not only to increase thrust 
performance and develop new propulsion techniques, but at 

. the same time to reduce costs all along the line. 


J 34 of America’s 38 successful satellites and space probes have been 
L launched by Rocketdyne engines. 


FIRST WITH POWER FOR OUTER SPACE 


ROCKETDYNE tz 


DIVISION OF NORTH AMERICAN AVIATION 


Ls Canoga Park, California; Neosho, Missouri; McGregor, Texas 
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This portable, lightweight package keeps 
the pilot cool and comfortable on the ground 
in his flight pressure suit from one to two 
hours. Utilizing cryogenic liquid oxygen as 
the coolant, it requires no electric power or 
other connecting supplies which might 
hinder the pilot’s mobility. 

Easily carried by hand or slung from the 
shoulder, the AiResearch unit can cool a 
pilot wearing either full or partial pressure 
suit during travel to and from his aircraft, 
preflight checkout and while seated in the 
cockpit. Pure oxygen for prebreathing can 
also be provided as a simultaneous function. 

This extremely simple and reliable cool- 
ing unit has no moving parts. In operation, 
ambient air vaporizes a supply of liquid 


provides maximum pil 0. 4 a” oxygen to pressurize the system. Cooling air, 
comfort on the ground 4 y made up of stored oxygen and ambient air, 


». is then circulated through the suit. 


AiResearch is also in production on self- 
contained life support systems inside fully 
enclosed protective suits. These suits allow 
the wearer to work safely in hostile environ- 
ments such as toxic missile fuel handling 
and fire fighting. Research for modification 
of these systems is now being conducted for 
space use. 


¢ Please direct inquiries to Los Angeles Division 


cc 
AiResearch Manufacturing Divisions 


Los Angeles 45, California « Phoenix, Arizona 
Systems and Components for: AIRCRAFT, MISSILE, SPACECRAFT, ELECTRONIC, NUCLEAR AND INDUSTRIAL APPLICATIONS 
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EDITORIAL 


The Forward Look — 


A: with most nonprofit organizations in the United 
States (and also with a lot of others which are not 
intended to be run that way!), your Institute has been 
concerned with the problem of ‘‘making both ends 
meet” in a period of continually rising costs. This 
has been a matter of considerable interest to the 
Council and a subject of close and continuous study 
by Management over the past several years. As 
those of you know who have studied the Operating 
Statement for Fiscal 1960 (published recently in 
the JA S News), this has led to an operating deficit for 
the year ended September 30 last—the first really 
“red ink”’ figure in over 10 years. 

This is basically a reflection of the inflationary 
trends of the times. Everyone knows that a dollar 
buys less today than it did 5 or 6 years ago... As 
many people have discovered who depend on rela- 
tively fixed incomes, the ‘‘outgo’’ side begins to out- 
run the ‘“‘intake.”” Management of any enterprise is 
well aware that the costs of all goods and services 
have been mounting steadily since World War II. 
We are feeling the effects of this trend at a time when 
the services of the IAS are being augmented to meet 
the shifting professional needs of our members. 

The day before the opening of the Annual Meeting 
(January 22) the Executive Committee met to 
consider, among other things, the budgets for Fiscal 
1961 and 1962. All aspects of our situation were 
discussed, and further studies in detail were re- 
quested from Staff as a basis for future Council 
consideration. The principal question to be re- 
solved concerned when, where, and in what amounts 
IAS reserve funds are to be committed to help tide 
over the current period of “‘recession.”’ 

Somewhat to our amazement, the following day 
turned up a gaggle of rumors in the Astor’s cavernous 
corridors to the effect that we had “‘gone broke’’— 
that the Institute was bankrupt !—that we were, like 
the Arabs, about to fold our tents and as quietly 


+ 


steal away. We brushed off the first such report as 
too ridiculous for a second thought, but when other 
similar bits of intelligence fell off the grapevine it 
seemed clear that someone was really “passing the 
word.” After that, whenever the subject was 
raised (including the meeting of IAS Section Chair- 
men and Regional Advisory Committees) we took 
pains to explain carefully, and in as much detail as 
anyone wanted, that, although we do face a some- 
what unprecedented budgetary imbalance which 
deserves the most serious study and consideration by 
Council and by Management, there is nothing 
in the picture that is critical or remotely cata- 
strophic! 

The situation was again reviewed by the Director 
and the Executive Committee of the Council at a 
meeting on February 14. Ways and means of cor- 
recting the unbalanced condition of current budgets 
were discussed in great detail. But no one suggested 
that we deviate in any way from the programs and 
operating policies laid down and approved by the 
Council during the past 2 years for any aspect of 
IAS programs. The consensus was clear—that we 
should continue to improve and expand all essential 
IAS services to our membership as their professional 
needs require. Management has a clear-cut man- 
date to proceed along lines already established, as 
effectively and as efficiently as possible. We know 
that we cannot expect a balanced budget operation 
in the current, or in the next, Fiscal Year, but in the 
thinking of your Officers and Council, the main- 
tenance of IAS services is the paramount issue. 
They are confident that ways and means will be 
found to close the budget “gap” in the not-too- 
distant future. Meanwhile, the word is to Carry 
On at full throttle—to strengthen in every possible 
way the Institute’s position as the outstanding pro- 
fessional society in the aerospace field. 

—SPJ 
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Vehicle Design—Satellites 


The Inflated Satellite 


Characteristics in Sunlight and Darkness 


Herman Mark, AFIAS, and Simon Ostrach, AFIAS 


Lewis Research Center, NASA 


— THE VARIOUS TYPES Of satellites being 
considered are those made of lightweight plastic 
which can be inflated pneumatically. It has been 
determined! that an aluminized hollow plastic sphere 
with a diameter of 100 ft in a terrestrial orbit at an 
altitude of approximately 1,000 miles would be 
easily visible to the naked eye, would be a good re- 
flector of radar waves, and would be protected from 
deterioration caused by radiation. On the basis of 
this kind of information, the NASA communications 
satellite (Project EcHo) was developed. 

The performance and life of an inflatable satellite 
also depends on its retaining its shape, which is 
governed by the pneumatic process. This process is 
temperature-dependent; and, since Earth satellites 
are subject to a range of thermal conditions in the 
course of their orbits, an independent study was made 
of the behavior of an inflated satellite as it passes 
from sunlight to darkness and back into sunlight. 
The essential considerations of that investigation 
are presented in this paper. Although the Project 
Ecuo satellite was used as a model for the present 


Dr. Mark received a B.S. degree in me- 
chanical engineering from CCNY in 1944. His 
U.S. Army Air Force assignment was to the 
Aircraft Engine Research Laboratory, NACA, 
Cleveland. Upon discharge, the Laboratory 
kept him engaged in combustion research 
and on jet engine comp t development. 
Receiving a Case Institute of Technology 
M.Sc. degree in June 1953, his entry into 
full-time study at Cornell University Graduate 
School of Aeronautical Engineering culmi- 
nated in a Ph.D. in 1957. Dr. Mark is, at 
present, head of Space Environment Section, 
Lewis Research Center, NASA. 
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study, the analysis presented herein is of more 
general application. 


Destructive Forces 


The satellite can be destroyed if either the internal 
pressure becomes so large that it will burst or so 
small that it will collapse under external forces. The 
pressure required to burst a spherical satellite is 
calculated from the relation 


P = (26/R)S (1) 


where P is the bursting pressure, 6 the sphere thick- 
ness, XR the sphere radius, and S the tensile strength 
of the plastic material. For a Mylar sphere 1 /2-mil 
thick, with a 50-ft radius and tensile strength of 10* 
Ib/sq in., the bursting pressure is 0.86-mm Hg. 

The two most important external forces which 
could collapse the sphere were found to be the stag- 
nation pressure and the radiation pressure. On the 
basis of Newtonian flow (differing by a factor only 
from free-molecule flow), the stagnation pressure P; 
acting on the leading surface of the sphere is deter- 
mined from the relation 


Dr. Ostrach is a Professor of mechanical en- 
gineering at the Case Institute of Technology 
in Cleveland, Ohio. He received a degree 
of B.S. in mechanical engineering in 1944 
and a degree of M.E. in 1949 from the Uni- 
versity of Rhode Island. Further studies 
brought him a Sc.M. in 1949 and a Ph.D. in 
1950 in applied mathematics from Brown 
University. Before joining the faculty of the 
Case Institute, Dr. Ostrach was Chief, Fluid 
Physics Branch at the Lewis Research Center 
of NASA. He is the author of a number of 
publications in the fields of fluid mechanics 
and heat transfer. 
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Investigation results are presented of an inflated satellite as it 
passes from sunlight to darkness and back 

again into sunlight. Although the Project Echo 

satellite was used as the model, the analysis obtained is 


applicable in more general categories. 


nore = Pi + pV? (2) 


where P; and p; are atmospheric values (extrap- 
olated from available data’). The velocity of a 


satellite in a circular orbit is computed from 
ma 


r so V=Rz Vg (Re + h) (3) 
The where Ry is the Earth’s radius, g is the acceleration as | | 
te 1S due to gravity, and h is the height of the satellite rn © | 
above Earth. Stagnation pressures computed at | 
(1) 186, 500, 900, and 1,000 miles from Eqs. 2 and 3 are sel | | 
presented in Table 1. | 
hick- The radiation pressure due to the sun at the 340}- | : 
ngth earth’s distance from the sun is given in reference 3 een 
2-mnil as 4.5 dynes per em’, or 3.375 X 10~*-mm Hg. 320}- 
f 108 This will also be the average collapsing pressure on‘a ae \ 
me reflecting sphere when it is in the sunlight (Appendix —. : | x 
vhich I). In Appendix I there are developed explicit ex- ¥ i * 
stag- | pressions for the radiation pressure on the sphere re- ™ -/ Re 
n the sulting from influences of the sun and the earth. ee tae 
only The expression for the latter shows that the radia- 
re P, tion pressure on the sphere caused by the earth is Poo = 
leter- about '/, that of the sun and exceeds the stagnation 
pressure at 900- and 1,000-mile altitudes (Table 2). i 
Internal Condition 
200) 
The satellite is inflated by the vapor pressure of 
some material such as water, benzoic acid, or anthra- 180}— 
quinone which is stored inside the plastic skin. To ae eee te 
determine the internal pressure, it is necessary, first 
of all, to find the satellite skin temperature under vo, oe 
conditions of radiation equilibrium in the sunlight Fig. 1. Temperature-time history of anthraquinone vapor in inflated 
and darkness (earth’s shadow). For this purpose Solid 
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Operations 


T HE U.S. Nava Arr Test Facility (Ship Installa- 
tions) [NATF(SI)] is a field activity created by the 
former Bureau of Aeronautics for the specific purpose 
of evaluating catapults and arresting gear and related 
components. It was initiated in 1956 when facilities 
and space at the old Naval Aircraft Factory (N.A.F.) 
site, Mustin Field, became inadequate for the speeds 
and weights of modern jet aircraft. A test and 
evaluation activity independent of the research and 
development functions at N.A.F. had long been de- 


a a Fig. 1. Layout of the U.S. Naval Air Test Facility (Ship Installations), a tenant activity of the U.S. Naval Air Station at 
Lakehurst, N.J. 
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Aircraft Launching and Recovery Systems 


Recent Advances and Their Potential Application to Commercial Aircraft—Part | 


The first in a series of two articles describing the U.S. Naval Air Test Facility (Ship Installations) 
referring to the overall description of the test sites and details of launching systems. 


Capt. William C. Fortune, AFIAS, U.S. Naval Air Test Facility (SI), Lakehurst 


sired. The test facility is under the military control 
of the Naval Air Research and Development Ac- 
tivities Command, Johnsville, Pa., and under man- 
agement control of the recently formed Bureau of 
Naval Weapons, Washington, D.C. 

Largest user of the various test sites (see Fig. 1) 
is the Naval Air Engineering Facility (SI) [NAEF 
(SI)], Philadelphia (formerly called the Naval Air- 
craft Factory mentioned above), which has the ma- 
jority of components in development status. All 
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American Engineering Company, E. W. Bliss Com- 
pany, Van Zelm Associates, Hydro-mill, and others, 
also have equipments under assessment. A Swedish 
barricade built by Borg Fabrik Aktiebolog (BEFAB), 
for possible Marine and NATO applications, will be 
evaluated soon. Emergency overrun devices are 
being tested for the Air Force. The Federal 
Aviation Agency is expected to join in with related 
efforts in this field involving jet transport type air- 
craft (KC-135). The Flight Safety Foundation, 
Inc., has provided certain operational data which 
has been of value in appraising the applicability of 
the various systems to commercial usage. 

A different utilization for the jet car track sites is 
being planned by the Air Crew Equipment Labora- 
tory at Philadelphia, for investigation of the medium- 
speed portion of their test envelope for ejection 
seats, escape capsules, etc. Shots up to 230 knots 
can be made with little modification; study has been 
initiated to see if 400-500 knots might be attained. 
The NATF(SI) has the only such potential on the 
east coast, and can provide normal service accelera- 
tions, normal escape speeds, and low test run charges, 
having jet-engine-driven sleds instead of supersonic 
rocket vehicles such as are used at west coast tracks. 

To accomplish the various tasks within the scope 
of its mission, the test facility operates the following 
test sites: 

(1) Recovery Systems Track Sites. Five double- 
rail jet car tracks are used for testing various arrest- 
ing gears or components, plus barriers, barricades, 
and their components. In addition to arresting gear 
designed for shipboard use, shore-based emergency 
arresting systems are installed, development and 
evaluation programs being conducted with one or 
the other to keep site productivity high. The Mark 7 
Mod 2-3 hydraulic arresting gear, now operational in 


Fig. 2. Four J-33 jet engines ready to propel a dead load at track No. 1. 
the kinetic energy through fluid or mechanical means. (Official photograph, U.S. Navy.) 


Arresting gear stops travel of deadload by absorbing 


the U.S.S. “Independence,” was evaluated here. Re- 
cent improvements in capability of the U.S.S. “Bon 
Homme Richard,” “Midway,” and “‘Ranger’’ have 
also been perfected at these tracks. Other recovery 
systems designed for operation under the energy- 
absorption principles of metal deformation, fluid 
friction, mechanical braking, and exchange of 
momentum are being evaluated. Large nylon nets 
for trapping aircraft, of NAEF(SI) design, as well as 
BEFAB, are being proved at the sites. Track 1 is a 
mile long, the other four 1'/2 miles in length. As 
noted before, light deadloads or aircraft sleds can be 
accelerated down these tracks to over 250 mph. 
Loads up to 100,000 Ibs can be tested at nominally 
slower speeds. Four J33 engines normally power 
each jet car (Fig. 2); however, one of the group has 
been converted to J48’s. For the high-energy tests 
scheduled—i.e., over 70 million ft-lbs—JATO boost- 
ers are added to the deadloads. Trailing behind the 
jet car, on each rail, looking like waffles stacked on 
their edges, are brake assemblies. These normally 
roll along the top of the J (Continued on page 38) 


Capt. Fortune was graduated from the 
Naval Academy as a commissioned ensign 
in 1933. In Naval Aviation since 1935, he 
has served on U.S.S. “Texas,” “Ranger,” 
and “San Francisco.” Capt. Fortune was 
awarded an M.S. in aeronautical engineer- 
ing (M.LT.) in 1942. He served on Staffs of 
Commander Air Forces Pacific Fleet, during 
World War Il and Commander Fleet Air, 
Japan, during Korean incident, and was 
awarded the Commendation medal (both 
theaters). Capt. Fortune was Chief Engineer 
at Naval Aircraft Factory, 1946-48, and 
has had various R&D assignments, including Aerodynamics Lab. at the 


_David Taylor Model Basin; Air Branch, ONR; BurAer General Repre- 


sentative, Western District and Comp t Development Officer (direct- 
ing Power Plants, Avionics and Airborne Equipment Divisions), in the same 
Bureau in Washington; and at present is Commanding Officer at the 
NATF(SI), Lakehurst, N.J. He is a member of the "N” Club and the 
Caterpillar Club and is an Associate Fellow of the IAS. 
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I. THE ENVIRONMENT EXTERIOR of the earth’s atmos- 
phere, man must have protection against the hazards 
of tissue-damaging radiation from electrons, protons, 
neutrons, photons, and other high-energy particles. 
The degree of protection in the form of material 
shielding will be a function of the residence time in 
space, and the type, quantity, and energy of the 
radiation encountered in the trajectory or orbit de- 
sired. 

At present the available information regarding this 
extraterrestrial radiation is sketchy—in many areas 
of interest, information does not exist at all. A tre- 
mendous effort in data gathering and compilation 
will have to be made before we know with reasonable 
certainty the average values of the types and ener- 
gies of this radiation as a function of altitude, lati- 
tude, and longitude. The availability of this infor- 
mation will facilitate the selection of orbits, shield 
designs, and payloads for specific missions. A need 
exists for the breakdown of the number of particles 
per unit energy interval as a function of energy. 

In addition to the data described, access to a ref- 
erence, or handbook, showing the attenuation, deg- 
radation, and transmission of the several types of 
radiation through a variety of structural materials 
(as well as the common elements) as functions of 
thickness would be of tremendous assistance to the 
design and shielding engineer. This information 
plotted in the form of curves, together with data re- 
garding the kind and intensity of the radiation ex- 


This work was performed in connection with Air Force 
Contract AF 33(616)-6641 under the auspices of the Flight 
Dynamics Laboratory, Wright Air Development Division, 
Wright-Patterson AFB, Ohio. 
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Fig. 1. Range-energy relations for electrons in aluminum. 


Mr. Barbieri is a Research and Development 
Scientist in the Space Experimentation Sec- 
tion at Aeronutronic, coordinating space ex- 
periments for the USAF TS 609A program. 
He has a B.S. (physics) from George Wash- 
ington University and did graduate work at 
University of Maryland. While at the 
Radiation Division of NRL, he worked in the 
fields of gamma-ray and neutron spectros- 
copy and helped develop techniques _in- 
volving heat transfer from metals to liquids 
for passive and turbulent boiling conditions. 
Mr. Barbieri has developed an analytical 
method for determining the relative concentration of uranium and plu- 
tonium isotopes as a function of power output of a nuclear reactor; he 
coauthored a TID document on the physics and economics of plutonium 
production and has applied these techniques in other studies. 


Dr. Lampert is Manager of the Aeromech- 
anics Dept. at Aeronutronic. He has a B.S. 
degree from Georgia Institute of Technology 
and a Ph.D. from California Institute of 
Technology. He joined Aeronutronic in the 
spring of 1957 and has been involved 
principally in problems related to aerospace 
vehicles. Prior to joining Aeronutronic he was 
affiliated with NACA (Ames Lab.), JPL, 
NOTS, and University of Southern California 
in various capacities. He has taught mathe- 
matics at Georgia Tech and engineering 
courses at University of Southern California. 


Table 1. 


Electron Absorption in Aluminum 


Thickness to 


Electron —Absorb Completely— 
Flux Material 
(cm~2- Density Aluminum 
Energy Range sec!) (ing/cm?) (cm) 

30 kev—190 kev ~10! 40 0.015 
190 kev—460 kev ~10? 170 0.063 
460 kev—1 mev ~10 400 0.15 
Above 1 mev None 
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Vehicle Design—Manned Spacecraft 


The Interdependence of Manned Spacecraft 
Design and Radiation Shielding 


“There is no reason why a 4-week flight is not possible 


from a radiation hazard viewpoint.” 


L. J. Barbieri and S. Lampert, Aeronutronic Div., Ford Motor Co. 


pected, would enable the designer to select an opti- 
mum satellite wall and shield design. 


Shielding for the 500-Km (310 Miles) Orbit 


The shielding requirements for an orbit 500 km 
above the earth’s surface are determined from data 
presently available. Alterations may be necessary 
if subsequent detailed investigations reveal signifi- 
cant deviations from the meager data so far obtained. 


Electron Shielding 


General Discussion—The energy ranges and fluxes 
shown in Table 1 were supplied by the Air Force 
Special Weapons Center (AFSWC) in a private 
communication (March, 1960). From the table we 
see that 0.15 cm of aluminum will completely shield 
against 1-mev electrons. Even if the satellite walls 
total considerably less than this thickness, any 
padding on the inside walls, air within the capsule, 
and the clothing worn by the occupants will provide 
ample shielding against electrons. Any brems- 
strahlung produced by the incidence of these elec- 
trons on aluminum will be negligible. 

Interactions With Matter—Electrons dissipate their 
energy by two distinct mechanisms!~* when passing 
through matter—ionization and _ bremsstrahlung 
radiation. Jonization is the process by which a 
high-energy charged particle collides with the atomic 
electrons of the substance through which it passes 
and thus distributes its energy among them. Ioniza- 
tion involves the planetary electrons and as such 
is rather insensitive to the atomic number of the 
stopping material. Bremsstrahlung radiation is elec- 
tromagnetic energy released by the high-energy 


charged particle through its interaction with the field 
of an atomic nucleus. It is a function of the mass 
and energy of the particle, as well as of the atomic 
number of the target nucleus. 

Range-Energy Relationship—Several empirical re- 
lations between the practical range and energy of 
an electron have been proposed! by various work- 
ers. From 0.01 to ~3 mev (Continued on page 41) 


Fig. 2. lonization and radiation losses for electrons in air, Al, and Pb. 
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Vehicle Design—V/STOL Aircraft 


V/STOL Logistic Transport 


Lifting Fan Design Selection and Performance 


Elliot Kazan, AIAS, and Wilbur Bergen 


Republic Aviation Corporation 


F THE MANY SYSTEMS and schemes being con- 
sidered to provide aircraft with the capability of ver- 
tical take-off and landing (VTOL), the lifting fan 
design, in combination with turbojet engines, offers 
the advantage of high-speed cruise with minimal 
VTOL ground erosion problems. This advantage is 
particularly important in logistic or assault opera- 
tions where maximum utilization from austere sites 
is a paramount consideration. 

This paper discusses the highlights of a parametric 
design and performance study of a 2-ton payload 
logistic transport with lift fans installed in the wings. 
Analysis is made of mission capability with various 
propulsion systems and VTO weights, as well as the 
trade-offs possible with a selected power plant design. 


Lifting Fan System 


The fan system considered in this study if a self- 
contained convertible propulsion system of the type 
currently being tested by the General Electric Com- 
pany (GE). The basic components of a typical fan 
system are illustrated in Fig. 1. The gas generator 
is equipped with a two-way diverter valve to enable 
utilization of the exhaust gases from a conventional 
jet nozzle and to provide coupling between the gas 
generator and lift fan for the vertical lift mode of 
operation. In the lifting mode, the highly loaded 
fan (300 to 400 lb/ft?) is driven by a tip turbine 
which receives generator exhaust gases through a 
nozzle diaphram in the scroll. To provide vectored 
horizontal thrust for the transition from vertical to 
the horizontal flight mode or for STOL (short take- 
off and landing) capability, exit louvers are provided 
at the outlet of the fan. 

With this propulsion system, the basic gas gen- 
erator thrust has been shown to be augmented by a 
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Fig. 1. X353-5 propulsion system. 
factor of almost three. Coupled with this augmenta- 
tion ratio, the tip turbine-fan system achieves ver- 
tical lift with lower exhaust gas downwash velocity 
and at a considerably lower temperature than would 
result from the use of direct-action turbojet engines. 
A comparison of these characteristics between several 
currently popular turbojet approaches to VTOL is 
shown in Fig. 2. The lift-fan system with an average 
exhaust gas temperature of 180°F, an exhaust veloc- 
ity of 550 ft/sec, and an augmentation ratio of 
almost three, appears potentially to be an optimum 
propulsion system for high-performance VTOL air- 
craft. This system also provides a growth potential 
which minimizes the possibility of early obsolescence. 
Since the tip turbine lifting fan concept is a rela- 
tively new approach to achieving VTOL, it is felt 
that conventional aircraft design techniques must 
be re-analyzed to ensure optimization of designs 
using this propulsion system. The following discus- 
sion presents the highlights of a detailed parametric 
study in which basic design variables were re- 
analyzed for application to a V/STOL transport. 
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Combined with turbojet engines, the lifting fan design offers 


@ High-speed cruise 
® Minimal VTOL ground erosion 


Highlights of a parametric design and performance study of a 
2-ton payload logistic transport with wing-installed 


lift fans are detailed. 


V/STOL Logistic Transport Parametric Study 


The optimum design of a V/STOL subsonic trans- 
port requires a complete knowledge of the integrated 
role played by the power plant system, the aero- 
dynamic configuration, and the payload-mission re- 
quirements for the aircraft. The degree to which 
each of these factors influences the size (weight) of 
the aircraft will indicate the direction for design 
compromise without significantly sacrificing perform- 
ance capability. If the primary capability is VTOL, 
emphasis must be concentrated on maximizing the 
fuel/gross weight ratio for a minimum-gross-weight 
aircraft performing the desired mission. Depending 
on the mission requirement, and even the mission 
calculation rules, some or all of these factors play a 
different role in designing a VTOL transport as com- 
pared with designing a conventional transport with 
the same mission capability. The low-drag, high- 
aspect ratio wing used in conventional transport 
design, must be reconsidered in light of the VTOL 
weight requirement. Again, a VTOL design based 
on a power plant with high thrust-to-engine-weight 
ratio may not necessarily yield an optimum aircraft. 
High specific fuel consumption, together with a 
heavier structure and a higher drag configuration 
required to house the system, could easily negate 
the advantages of the power plant with the best 
cruise specific thrust. 

Based on these considerations, the parametric ap- 
proach taken to design an optimum two-ton payload 
V/STOL transport attempts to determine the rela- 
tive role played by the power plant, structural 
weight, and aerodynamic configuration as each af- 
fects the mission capability. The study considers 
configurations designed for VTOL as the prime op- 


erational requirement with STOL a secondary capa- 
bility resulting only from a need for extended range 
or increased payload requirements. The VTOL de- 
signs considered are limited to configurations having 
lifting fans installed in the wing panels. The number 
and size of these fans, together with the gas genera- 
tor requirements to provide the VTOL and cruise 
capability as a function of VTO gross weight, deter- 
mine the extent of the parametric study. 

The approach used in the parametric study is 
summarized in three analysis (Continued on page 47) 


Mr. Kazan studied at the Polytechnic Institute 
of Brooklyn, where he received his B.S. 
(1949) and M.S. (1950) degrees in Aero- 
nautical Engineering. He is also a graduate 
of the Oak Ridge School of Reactor Tech- 
nology (1957-1958). As an aerodynami- 
cist at Boeing Airplane Company (1950-52), 
he conducted stability and control, drag, 
and load analysis in the initial development 
of the Bomarc interceptor. At Republic Avia- 
tion he has worked in the field of aero- 
dynamic stability and control, drag, aero- 
elastic load analysis, and design require- 
ments for preliminary design proposals. Mr. Kazan has participated 
in numerous ballistic missile and V/STOL aircraft studies for the Armed 
Services. At present, he is group leader in charge of aircraft and 
missile aerodynamics in the Applied R&D Aerodynamic Section of 
Republic Aviation. 


Mr. Bergen, the coauthor of this paper, is a 
Group Aerodynamics Engineer with Re- 
public Aviation Corporation's Research and 
Development Division. He became affiliated 
with the company in 1946 and since that 
time has been engaged in work on such 
projects as the F-84 fighter bomber and the 
XF-91 and XF-103 interceptor aircraft. He 
has had broad experience in the field of 
preliminary aircraft design, including 
subsonic as well as supersonic transports. 
His more recent efforts have been 
directed toword the development of the 
V/STOL aircraft for utilization in various military applications. 
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Fig. 1. Elliptical hole in a biaxiaily stressed plate. 


Wis ONE THINKS of fabrics, he does not ordi- 
narily think of structures; yet we are living in an 
era when the structural uses of fabrics are multiply- 
ing. Already they include airships, tires, inflatable 
shelters, life rafts, radomes, and such unlikely things 
as dams, airplanes, and radar antennas. Generally, 
fabric structures require another structural element, 
such as air, which maybe you hadn’t thought of as 
structural either. Actually, inflated fabric structures 
have been with us for a long time. The era began 
in 1783 with the first balloons: the Montgolfier 
brothers’ balloon was linen lined with paper and in- 
flated with hot air. The Robert brothers made a 
great advance by coating silk cloth with a rubber 
varnish so that it would hold hydrogen—this nearly 
forty years before Mackintosh made a decent rain- 
coat. But there seems to have been little scientific 
study of fabric structures until after the advent of 
the pressure airship. In 1907, the German Society 
for the Study of Airships! laid out a program for the 
study of the relative merits of single-ply, parallel- 
doubled, diagonal-doubled, and various 3-ply cotton 
fabrics, also linen, silk, and goldbeater’s skin. By 
1909, papers were published in France and England, 
and in 1912 there appeared two landmarks in the 
history of fabric stress analysis: one by Haas and 
Dietzius? in Germany and one by Booth and Hyde*® 
in England. The problems they faced are still our 
problems today; that is, strength and stiffness. 


This hitherto unpublished work was presented at the 
Rubber and Plasics Conference of the American Society of 
Mechanical Engineers, Erie, Pa., October 10-12, 1960. 
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Structures and Materials 


An Introduction to 


A. D. Topping, MIAS 


Goodyear Aircraft Corporation 


Why Fabrics? 


Why does anyone use fabrics structurally in the 
first place? There are several reasons. Historically, 
the first reason was the lightness of a fabric balloon 
as compared to one made of thin sheet metal, for 
example. (Sheet metal balloons were proposed 
more than once and at least one was built; it never 


Symbols 


S,, Sy = average principal stresses in a sheet or fabric 
with a hole 

a,b = major and minor semiaxes of an ellipse normal to 
the x and y axes respectively 


oz, 0, = local stresses at the edge of an elliptical hole; 
warp and fill stresses 

& = thread count per unit length along the other thread 
set 

m = number of threads cut by a slit 

n = effective number of threads at either end of a slit 
taking the load released by cutting m threads 

01, 02 = principal stresses 

y, 8 = angles made by thread sets, after loading, with the 

normal to the irrotational axis 

0 = angle between the normal to the irrotational axis 
and the nearest thread set before loading 

Tzy = Shearing stress along warp and fill threads 

Yzy = change in angle between warp and fill, shearing 
strain 

Gr, |= secant shear modulus for warp-fill axes 

u,v = dimensions of an element 


T = tension in a single thread 

R = radius of test cylinder 

A = beam deflection due to transverse shear 
V = transverse shearing force 

A = enclosed cross-sectional area 

p = inflation pressure 


Note: All stresses and moduli are in terms of force per unit 
length. 
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Biaxial Stress Problems in Fabric Structures 


The superiority of fabrics to isotropic sheet and films 
in tear resistance is given an analytical basis, 


and problems of shearing stiffness are investigated. . 


flew, but in 1929 an aluminum-alloy pressure-airship 
envelope was built and successfully flown). In tires, 
rubberized fabric’s ability to absorb repeated flexing 
and to recover from overloads has been the decisive 
factor. In airship envelopes, all these assets are 
valuable, and fabric structures have the additional 
advantage of being able to redistribute local over- 
loads. The property attracting the most attention 
just now is fabric’s ability to be folded or collapsed 
into a small space such as that available in a rocket- 
boosted vehicle, then expanded at will into a full- 
fledged structure by a gas, or, to mention a recent 
development, by a foam which solidifies and forms a 


Dr. Topping received B.S. degrees in science 
and civil engineering from Purdue Univer- 
sity, an M.S. in civil engineering from Missouri 
School of Mines and Metallurgy, and a 
Ph.D. in theoretical and applied mechanics 
from the University of Illinois. He also has 
taught engineering mechanics at each of these 
schools. Since 1951, he has been at Good- 
year Aircraft Corporation where he is now 
in charge of structural analysis research and 
development. He has published papers and 
reports on teaching techniques, deep well 
stresses, multiaxial strength criteria, spot- 
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Wires SPACECRAFT TRAVEL as commercial carriers 
between the planets and the trip that goes to all 
the Galilean satellites of Jupiter has been given the 
name ‘‘The Milk Run,” the crews of those fantastic 
ships will doubtless still be arguing the role of man 
in space. When the adaptive behavior and the crea- 
tive capacities of man-made machines surpass the 
capabilities of both Edmund Hillary and Pablo 
Picasso, the men who built them will still point out 
that no man-made machines can surpass the crea- 
tive, adaptive, and imaginative abilities of their 
creators. And, when genetic control has produced a 
race of beings whose capabilities are so far beyond 
ours that even the word “human” no longer really 
applies, these superhumans will still argue that 
machines can do anything they can do in space 
better, cheaper, and more reliably. 

Of course, things may turn out differently than 
we expect, so that this fanciful prediction does not 
come true at all. However, if we use man’s past 
history as a guide in attempting to predict his future 
actions, we find ample precedence for the endless 
continuation of similar arguments. The weight of 
history is strong enough, at least, to prevent me 
from attempting to give a final solution to the ques- 
tion of the role of man in space. 

The best that I can hope to do is to illuminate 
the problem—to bring into some sort of perspective 
the various aspects of the question which may some 
day lead toward (if not to) a solution. 


Why Explore Space? 


What are the reasons why man should go into 
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Space Man vs. 


A. R. Hibbs, Jet Propulsion Laboratory 


space? And, on the basis of the reasons, should 
his machines act as his temporary predecessors or 
his permanent substitutes? More fundamental still, 
why are men interested in space at all? 

If we proceed into space one step at a time, we 
come first to the domain of the artificial earth satel- 
lite. In this volume of space, the utility of our 
space effort is comparatively easy to establish. 
Satellites capable of predicting the weather or relay- 
ing communications can, by their services, more than 
pay for themselves, as has been repeatedly shown in 
analyses of these commercial applications. The 
military utility of surveillance satellites is also 
commonly accepted. 

The picture is somewhat different if we extend 
our range to include the moon. Here, we have far 
exceeded the useful altitude for orbits of weather 
observation or communication satellites. We have 
reached a range from earth where few military ac- 
tivities would be either feasible or desirable. But 
let us not argue this question of the utility of lunar 
bases. 


Dr. Hibbs has a B.S. degree in physics 
(1945) from California Institute of Tech- 
nology, an M.S. degree in mathematics 
(1947) from the University of Chicago, and a 
Ph.D. in physics (1955) from CIT. He joined 
Jet Propulsion Laboratory in Feb. 1950, and 
r d on completing his Ph.D. The Divi- 
sion of the Space Sciences, under his direc- 
tion, is responsible for development of in- 
struments to carry out scientific aspects of 
moon and planets explorations as assigned 
to JPL by NASA. 
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Space Machine 


Must man actually land on distant plaents in order to 


fulfill our intentions in space? 


Or is this just an indulgence of our national 


urge toward human adventure and excitement ? 


Although the military utility of the moon will 
likely develop as the technique of space flight de- 
velops, this probably does not apply to the planets— 
even the near ones—for a long time to come. How- 
ever, we anticipate sending both unmanned and 
manned spacecraft to the planets without waiting 
for some possible far distant commercial or military 
utility to justify these undertakings. 

So, in order to get more to the heart of our pres- 
ent problem, suppose we confine the discussion to 
planetary flights. Then we will rot be bogged down 
with arguments about commercial desirability or 
military necessity. Instead, we can look directly at 
the question, “‘Can machines do all of the things that 
need doing on a planet, or must live human beings ac- 
tually land on planetary surfaces in order to fulfill 
our intentions in space?” 


Machines Only? 


Early in any discussion of this problem, one cus- 
tomarily hears two typical points. First, there is 
the point that automatic machinery can make all 
of the necessary measurements required in plane- 
tary exploration. Furthermore, they can be made 
somewhat adaptive to their environment and can 
certainly respond to command from earth to change 
routines in the face of unexpected developments. 
Machines generally demand much less complexity of 
their carrying spacecraft and launching rockets as 
well as less auxiliary equipment to maintain them 
in functioning condition. Furthermore, of course, 
they can be abandoned on the planetary surface. 
Although in principle man, too, could be abandoned, 


once he had served his mission, it is doubtful that a 
tax-paying public would willingly finance a mission 
wherein this outcome was anything other than a 
highly improbable‘accident. 

The argument in support of the machine continues 
by pointing out that, even in an earth-based labora- 
tory, man has by now far exceeded those areas of 
investigation where he can with his own eyes and 
ears perform the necessary observations. Even on 
earth man needs a whole battery of machinery to 
carry out any sort of an exploratory mission. Cer- 
tainly, the same would be true in space and on the 
surface of a planet. Thus, even if man should go to 
the planets, machinery would still be necessary. In 
that case, why not leaye the man at home, where he 
is not only safer but also does not use up valuable 
weight and power allotments? 


Manned Flights 


The other side of this argument is, of course, in 
support of manned flights. The point is often made 
that no matter how adaptive and creative automatic 
machinery is made it can never surpass in creativity 
and adaptivity the abilities of its makers. Cer- 
tainly, in this new and completely unexplored en- 
vironment, on the surface of an only partially known 
planet, the unexpected will be routine. Only a man 
on the spot will be able to adapt sufficiently well to 
preserve the usefulness of the mission. Although re- 
mote control is possible, the round-trip time for a 
signal (a sensed occurrence from planet to earth and 
a response command from earth to planet) takes 
several minutes. There will (Continued on page 59) 
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Aerodynamics and Fluid Mechanics 


Effects of Boundary Layer and Geometry on 
Characteristics of Perforated Walls for 


Transonic Wind Tunnels 


J. Lukasiewicz, AFIAS, ARO, Inc. 


W ane THE USE of mixed-boundary (slotted or Symbols 
perforated) walls effectively eliminates wind-tun- At = perforated test-section entry area (= nozzle 
nel blockage at high subsonic’ Mach numbers, exit area) 
the elimination of wall interference at low supersonic Ay = surface area of perforated walls 
speeds presents a much more difficult problem. Cp OF (Pe — Ppienum)/o» pressure 
coefficient 
Ideally, the shock and expansion-wave cancellation Cop = BC p80 — slope of (Cp, 6 wall characteristics (per 
properties of the walls would have to be such as to radian) 
match exactly the requirements of any flow field D = model diameter 
produced by a model. h = test-section height (= width) 
In Fig. 1, the required ideal characteristics are l = length of. perforated walls 
M = free-stream, or test-section, Mach number 
indicated for the particular case of a 2-percent - ‘ ; 
a m = mass flow entering perforated test section 
blockage, 20-degree, cone-cylinder model at / - p = static pressure 
1.2, as calculated by Goethert (reference 1, Figs. Pplenum = pressure in the plenum chamber 
10 and 11). The theoretical flow field on the center- bo = stagnation pressure 
line of the tunnel wall is described in terms of pressure Po = free-stream, or test-section, static pressure 


ll 


porosity (based on @), percent 

t = plate thickness, in. 

x = distance downstream from model nose 

a = inclination of holes from vertical, degree (Fig. 


coefficient C, = (p — p..)/g. and direction of stream- 
lines 9. The ideal wall characteristics would ex- 
actly match the C,, @ plot of the flow field at the wail 


—i.e., would maintain the direction @ by allowing 5) (a = 0 unless otherwise stated ) 
the corresponding amount to flow through the wall 5 = boundary-layer displacement thickness 
6*,, 65% = boundary-layer displacement thickness at entry 
under pressure difference C,. It is evident from ; : : 
x : and exit to perforated section, respectively 
Fig. 1 that a single wall could not satisfy the theo- Stay = average 8+ 
retical requirements and, in practice, a wall char- aa? ee 
acteristic corresponding to what is considered an Am = net mass flow between test section and plenum 
= diameter 1oles in perforate ates (dri 
case considered (Fig. 1), two such alternatives are 
diameter for inclined holes), in. 
indicated : ” ch teristic and 
. > a linear average characteristic and a pu = mass flow per unit test-section cross-sectional 
area 
pv = mass cross-flow per unit perforated wall surface 
Publication of this paper for the first time in the open litera- 0 = inclination of streamlines, positive for flow from 
ture has been made possible by recent declassification. It test section into plenum chamber, radians 
was issued originally by the National Aeronautical Estab- Ow = slope of movable perforated walls, positive for 
lishment, Ottawa, Ont., Canada, in August 1957. diverging walls, minutes 
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Data reviewed suggest favorable interference characteristics 
obtainable with perforated walls having linear flow 
characteristics, provided boundary-layer displacement 


thickness does not exceed hole radius. 


Considering perforations 


caused flow disturbances, a satisfactory 60° inclined-hole 
perforated test-section design can be obtained with 
boundary-layer displacement thickness not exceeding 

0.5 percent of the smaller test-section dimension. 


linear ‘‘differential’’ characteristic, having different | 


slopes for outflow and inflow. 

In the case of a negligible boundary layer on a 
perforated wall, it has been shown theoretically and 
experimentally by Maeder (reference 9) that the 
perforated-wall characteristics are linear. On the 
other hand, with a thick boundary layer and a small 
flow velocity near the perforated surface, the per- 
forations may be treated as orifices yielding a para- 
bolic wall characteristic, C,a6?. 

It is therefore immediately apparent that, in prac- 
tice, the presence of the boundary layer would tend 
to make the perforated-wall characteristics non- 
linear, in a sense opposite to the desired one, indi- 
cated in Fig. 1. Generally, it will be seen that the 
wall characteristics can depend as much on the 
boundary-layer flow as on the wall geometry. The 
wall geometry needs then to be related to both the 
boundary-layer thickness and to the working section 
size, so that acceptable cancellation characteristics 
without excessive coarseness of the walls, which 
would cause appreciable disturbance at the model 
location, are obtained. 

It is with these aspects of transonic-tunnel wall 
design that this paper is concerned. 

The results are particularly applicable to cases in 
which the boundary layer is relatively thick at the 
entrance to the perforated-wall section, a situation 
which may be encountered in a blowdown tunnel 
working over a wide Mach number range. In this 
type of tunnel it is convenient, for mechanical rea- 


(Continued on page 62) 
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Fig. 1. Theoretical flow field of cone-cylinder model, 2 percent block- 
age, M = 1.2 (reference 1, Fig. 10, 11). 


Mr. Lukasiewicz is Chief of the von Karman 
Gas Dynamics Facility, ARO, Inc., at the 
USAF’s Arnold Engineering Development 
Center at Tullahoma, Tenn. He holds a 


the University of London and a Diploma of 
the Imperial College, London, in aeronautical 
engineering. From 1945 to 1948, he 
worked at the Royal Aircraft Establishment, 
Farnborough, in the field of supersonic aero- 
dynamics. Later, until 1958, he was Head 
of the High Speed Aerodynamics Labora- 
tory, National Aeronautical Establishment, 
Ottawa, Ont., Canada. Mr. Lukasiewicz is an Associate Fellow of the 
IAS, a Fellow of the Canadian Aeronautical Institute, and an Associate 
Member of the Institution of Mechanical Engineers. 
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Instrumentation 


Compliance of Nonuniform Gimbals 


Ira Cochin, Kearfoot Division, General Precision, Inc. 


Fur PURPOSES of analysis, gimbals are defined as 
either rigid or perfectly elastic. If an elastic gimbal 
is in the form of a short cylinder, and if it is loaded in 
tension or compression across a diameter, it flattens 
out somewhat. We can neglect loads normal to the 
plane of the ring since the cylinder is relatively rigid 
in that direction. Thus, the problem reduces to 
tensile or compressive loads, F,, or F,, on a circular 
ring, as shown in Fig. 1. Now refer to Fig. 2. Note 
that two deflections are produced by the load. The 
first is called 6,,, and is in the direction of the load; 
the second, 6,,, is at right angles to the first. (The 
compliances are, correspondingly, S,, and S,,.) 
The double subscript is used to define the direction 
of the deflection with respect to the load causing it. 
Thus, the first subscript defines the direction of de- 
flection, while the second gives direction of the load. 
The tensile or compressive load F may be due to 
reactions from adjacent gimbals, thermal expansion, 
or other sources. Note that the force F,, which acts 


Mr. Cochin is a unit head in the Systems 
Branch of Kearfott Division, General Pre- 
cision, Inc. His group has been responsible 
for such equipment as an inertial auto- 
navigator and a mirror landing system, and 
for new, unconventional products such as a 


rate gyro. Mr. Cochin, who has both the 
B.M.E. and M.M.E., also lectures at the 
Graduate Engineering School of The City 
College of New York on Guidance and 
Inertial Navigation. He holds, or has pend- 
ing, six patents, has published four articles, 
and will have a book on gyros published next year. His name ap- 
pears in “American Men of Science” and he is a member of Toast- 
master’s International. He has been elected to Pi Tau Sigma for scholar- 
ship and service, and to Sigmo Xi for an original contribution in the 
field of mechanics. 
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superconducting gyroscope and a vibratory ° 


Symbols 


A 


half the length of side normal to direction 7 of non- 
circular ring (Fig. 8) 


B = half the width of side in direction 7 of noncircular 
ring 

E = modulus of elasticity 

F; = tensile load on gimbal in direction 7 

Hy) = statically indeterminate tensile force at cut 

I, = moment of inertia of cross section of ring in region a 
(containing the bearing) 

I,, = moment of inertia of cross section of main portion 
of ring 

K, = coefficient of nonuniformity = J,,/J, in region a 

M = bending moment at any point of ring 

M; = bending moment due to load in direction 7 

M. = statically indeterminate moment at cut 

P.,, = bearing load resulting from P; 

P; = 1/2 bending load in direction 7 

r = (A+B +2p)/2 

R- = radius of circular gimbal, measured from its centroid 

S;; = compliance in direction 7 due to load in direction j 

U = strain energy ina beam = f(M?*/2EI)ds 

V, = statically indeterminate shear force at cut 

a = variable angle in analysis (Fig. 6) measured from 
vertical. Limits on a boundaries at which points 
a= 

B = (B + p)/(A + p) = noncircular geometric co- 
efficient 

6» = deflection resulting from force F 

6,4, = deflection resulting from unknown force H, 

6;, = deflection in direction 7 due to load in direction j 

¢ = p/r = noncircular geometric coefficient 

¢, = angle subtended by region a in circular ring 

% rotational deformation of ring at point 0 (Fig. 6) 
due to unknown moment M, 

p = radius of corner of noncircular ring 

Subscripts 
a = direction 1 
a = region a 
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It is relatively simple to determine deflections of a uniform circular gimbal. 
This paper discusses, derives, and tabulates formulas for 


evaluating the deflections of nonuniform, 


circular and noncircular gimbals in terms of the gimbal geometry. 


along the y axis, produces a symmetrical elongation 
6,, (1/26,, on each side of the center) and a sym- 


metrical contraction —6,; (1/2 6,, on each side of 


the center). 

If the same gimbal were supported on two bearings 
and a bending load (2P,) were applied in the plane of 
the equivalent ring and normal to the bearing axis, 
the gimbal would become slightly triangular, as 
shown in Fig. 3. The bearing reactions (P,,) act 
tangentially to the ring (defining the payload as 2P, 
simplifies the arithmetic). Now refer to Fig. 4. 
Note that the y deflections are not symmetrical 
while those in the x direction are. Curiously enough, 
the total diametric changes 6,,, and 6,, turn out to be 
exactly equal to those for the force F,, shown in Fig. 
2, when F, = P,,, that is, when the force F, is half 
the total payload 2P,. 


Nonuniform Circular Ring 


Now consider a ring whose cross section is not 
uniform, as shown in Fig. 5. The nonuniform ring 
is treated as a mutation of a uniform ring of thickness 
t. In several regions, the moment of inertia J, at 
the cross section differs from that at the main body 
In. In Fig. 5, the mutation (nonuniformity) is 
graphically illustrated by a thinner or thicker cross 
section in the region of nonuniformity. 

The ring is assumed to be vertically and horizon- 
tally symmetrical—i.e., = $6, = $5, 3 = 
= Ts = Ts, and J; = 

Briefly, the method of finding gimbal compliance 
is as follows. Cut the ring at point 0 and replace the 
separation by statically indeterminate internal loads 
Vo, Ho, and Mo, at the separation, as shown in Fig. 6. 


The unknown moment My, and the two unknown 
forces Hy and Vo, are determined for closure by using 
Castigliano’s theorein. Vo is equal to zero due to 
symmetry; Mp, and A are statically indeterminate, 
but can be found by determining the resulting de- 
flection bx, and the rotation @) as functions of My 
and Ho, and then solving these two equations simul- 
taneously. Once these are known, the remaining de- 
flections can be determined. 


(Continued on page 69) 


TENSION 


COMPRESSION 


Fig. 1. Ring in tension or compression. 
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Structures and Materials 


a stainless steel is anisotropic; it has 
different properties in different directions. Be- 
cause of this and since it has no definite yield point, 
it is difficult to predict accurately the load capacity 
of a stainless steel structural member. Use of the 
tangent modulus in formulas for predicting column 
load capacity gives results which, in general, are too 
low. 

Inasmuch as stainless steel derives much of its 
strength from cold rolling, it seems likely that when a 
stainless steel structure member is formed the mate- 
rial in its rounds, having been further cold worked, 
would become still harder. Although the effect of 
forming on rounds occasionally is mentioned in the 
literature, nothing has ever been done to evaluate 
it, and it was therefore considered worthwhile to 
conduct a limited test program to measure this 
effect. 

In addition to hardening the material in the bent 
portion of a member, the forming process will also 
introduce residual stresses in the corner material. 
Determining the exact stress at any one point in the 
bend under a load is quite difficult if not impossible, 
but it is possible to determine an average stress in 
the rounds. By conducting compression tests on 
fiat strips (representative of the unformed, fiat 
material in formed columns) and formed short 
columns, and comparing the stress-strain curves, 
it was possible to determine the hardening effect of 
cold forming. 


Effect of Forming 


Specimen Preparation 


The flat strip compression specimens, representa- 
tive of the unformed material in formed columns, 
were 2!/, in. by 1/2 in., and of thicknesses as shown 
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Effect of FORMING and WELDING on 
STAINLESS STEEL COLUMNS 


Robert E. Petersen and Axel O. Bergholm, The Franklin Institute Laboratories 


Mr. Petersen, who is a graduate of the 
Pennsylvania State College, joined the staff 
of The Franklin Institute Laboratories, Phila- 
delphia, for R&D in 1952. He is, at 
present, engaged in working on nuclear 
power reactors. He has been occupied in 
conceptual design, analysis of stresses re- 
sulting from physical loading, and analysis 
of distortion and stresses caused by steady- 
state and transient thermal gradients. Mr. 
Petersen has been concerned most with the 
methods of trying to prevent a change in 
the core volume resulting from fuel element 
subassembly bowing, and has been making experimental investigations 
into piping reaction stresses in a nozzle in a thin wall vessel. Before 
this work on the nuclear power reactors, Mr. Petersen—along with Mr. 
Bergholm—was engaged in the establishment of specifications for 
structural use of cold-rolled Type 301 stainless steel. Work on this 
project included the theoretical and experimental investigation of 
column loads, flat-plate buckling loads, and methods to be used in the 
analysis of stability of related structures. 


Mr. Bergholm's varied and numerous en- 
gineering experiences include many years in 
artillery and rockets development; stress 
analysis of the P-35 airplane wings and 
planning, supervision and construction of 
static test equipment for testing stainless- 
steel P-35 wing; design of test equipment 
for a railcar test plant, supervision of struc- 
tural testing of numerous railcar com- 
p ts; and pl d and tested structural 
components of the RB-1 and B-36 airplane. 
His work also has included design, construc- 
tion and testing of a stainless-steel model of 
a floating dry dock for the U.S. Navy to verify the design and analysis 
of the 1,200-ft prototype. As construction engineer for the Tri- 
borough Bridge Authority of New York, he was in charge of construc- 
tion of anchorages, steel erection and road construction of a suspen- 
sion bridge over Hell Gate in New York Harbor. With Mr. Petersen 
and others of the Franklin Institute staff, he developed specifications 
for structural use of cold-rolled Type 301 stainless steel. He is cur- 
rently on the staff of the Radio Corporation of America. 


in Table 1. The ends of each strip were ground 
parallel and checked to within 0.0003 in. on an opti- 
cal comparator. Identical strips were prepared 
from annealed, 1/4-hard, 1/2-hard, and full-hard 
steel. Table 2 shows the yield point and tensile 
strength of the steel. 
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The formed specimens of annealed, 1/4-hard, and 
1/2-hard steel consisted of two hat sections joined to 
form a closed column 3 in. long (Fig. 1). They were 
glued together to avoid the effect of welding. For 
full hard material, the formed specimens consisted 
of two U-shaped sections glued together to form a 
closed section 2!/2 in. long (Fig. 2). The formed 
specimens were so designed that, under load, they 
would not fail due to local crippling of rounds or 
buckling of flats. The ends were ground parallel to 
within 0.0003 in. 

For joining formed sections, Eutecto-Fillal glue 
proved very satisfactory. In several specimens 
tested, no failure of the bond was detected until 
after column failure. To use the glue properly, the 
surfaces to be joined must be heated to between 250 
and 300°F for several minutes. The glue is then 
applied, the surfaces clamped together, and the 
specimen baked in an oven for 24 hours at a tem- 
perature of 250 to 260°F. 

To learn whether this temperature treatment af- 
fected the stress-strain properties of the material, 
the unformed strips were given the same heat treat- 
ment as the glued formed specimens. These tests 
disclosed that for 1/4- and full-hard material the 
modulus of elasticity was improved but for the 1/2- 


hard material it was unaffected. This improvement - 


in strength results from stress relieving. 

It was not possible to determine the cross-sectional 
area of the formed specimens by direct measure- 
ment because they were machined after being 
formed. So the cross-sectional area of all specimens, 
including the flat strips, was computed from their 
weight, their measured height, and a density factor 
of 0.288 Ib./in.*. 

In the specimens made of hat sections, the cor- 
ners accounted for from 33 to 39 percent of the total 
cross-sectional area. This variation is due to differ- 
ences in thickness and bend radii. The percentage 
of corner material was computed from measure- 
ments taken from photographs of the ends of the 
specimens. The proportion of corner material for 
the U-sections of full hard material, measured 
directly, was 47 percent. 


Test Equipment 


The tests were performed on a Universal testitig 
machine with a maximum loading of 120,000 Ib. 
The flat strips were supported and prevented from 
buckling by steel guides (Fig. 3) made from draw- 
ings prepared by the U. S. Bureau of Standards. 
The flat sides of the specimens were lubricated so 
that no load was carried by the guides during the 
test. Two Tuckerman optical strain gages were 
used to record strain during the tests on both flat 
strip specimens and formed specimens (Fig. 4). 


Test Procedure 


In testing these compression specimens it was 
found that at higher stresses the transient creep 
was more pronounced; that is, a longer time had 
to be allowed for the strain gages to come to rest. 
It was not considered satisfactory to allow a set 
time period for every reading. Instead, readings 
were taken only when the creep rate fell below 
0.00001 in./in. change in two min. There was no 
need to wait for the first few readings on each speci- 
men, but as the test progressed longer periods of 
time were required between readings. For some of 
the last readings, it was necessary to wait 30 to 40 
min. for the gages to stop moving. (Cont. on page 71) 


Fig. 1. Column hat section. 


Research conducted by the Laboratories under American Iron and Steel Institute 
sponorship had as its objective the obtaining of data for preparation of 4 


design specification for Type 301 stainless steel. 


A method was derived 


for predicting stability of this type steel columns considering the 


nonlinearity of the stress-strain relationship, the effects of forming, spot welding 
and the direction of rolling. 
that are too low; 


Present column-strength formulas give results 
this method gives a more correct prediction. 
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Nuclear Applications 


Cosmic Rays, Nuclear Reactors, and 


Manned Space Systems 


David Reitz, The Martin Company, Denver Division 


ENERGETIC PARTICLES—cosmic rays— 
pervade all space. When they enter the earth’s at- 
mosphere they generally either penetrate it quite 
deeply, or produce secondary radiations when they 
interact with its constituents. Ionization caused by 
the primary cosmic rays or their products extends 
downward to all atmospheric levels, and even quite 
far underground. This ionization was first discovered 
and attributed to extraterrestrial sources in the 
years immediately preceding World War I. The 
interval between world wars was marked by wide- 
spread research, as the period following World 
War II also has been. Many fundamental insights 
have been obtained'~* concerning the nature of the 
cosmic radiation and its products. However, as vast 
and diverse as activities have been toward furthering 
descriptions of basic particulate phenomena, only a 
relatively few scientific workers have attempted to 
assess the biological implications of cosmic rays. 
As a result, even though the age of space travel is 
imminent, the susceptibility of man to these radia- 
tions is still unclear. It is certain, though, that 
ionizing radiations have harmful effects upon human 
tissue, and long-term ionizing cosmic particle ex- 
posures will be a part of prolonged space travel. 
Indeed, it appears that man will be making frequent 
flights above the atmosphere before adequate bio- 
logical experiments can be performed and results 
ascertained. In addition, the situation becomes 
somewhat more complicated if nuclear units are 
used in space applications as sources of electrical 


The author wishes to express his gratitude to Robert B. 
Demoret for his enthusiasm and encouragement, to Robert H. 
Edgerley for his helpful criticisms and suggestions, and to Jane 
Blizard for her assistance in the preparation of the final 
manuscript. 
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Mr. Reitz, 26 years old, received a Physics 
degree from the University of California, 
Berkeley, in 1955. He continued studies in 
physics at the University of Chicago from 
1955 to 1956. He joined the Douglas Air- 
craft Company's Nuclear Technology De- 
partment in 1956, and conducted investiga- 
tions associated with radiation conditions 
caused by nuclear power sources. Since 
1958, he has been with The Martin Com- 
pany as an Astronautics staff member and a 
Research Department member, where he 


4 has been analyzing space systems in re!a- 


tion to their space environment. 


energy or propulsive power. These have a relatively 
high power output, compared to more conventional 
systems, and thus can appear quite attractive. An 
unfavorable factor, however, is that nuclear reactors 
are radiation emitters and their radiations must be 
added to the hostile cosmic ray background for ex- 
tended flights above the atmosphere. It is the sum 
of these radiation fields which constitutes the total 
radiation intensity in space, and it is this very sum 
which is considered in the present paper. The de- 
velopment proceeds in such a way as to allow the de- 
termination of the most advisable orbital configura- 
tions. To begin, there is a statement of the maximum 
permissible radiation intensities to which man 
should be exposed. This is followed by descriptions 
of cosmic ray time variations and steady-state 
values. The allowable radiations from nuclear power 
sources are next given, and a method of obtaining a 
minimum weight reactor shield is outlined. Finally, 
this shielding method is applied to a specific 300-kw 
electrical auxiliary power reactor. 


Permissible Radiation Intensities 


When ionizing radiations enter human tissue, they 
often cause injurious changes in molecular structure. 
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Time variations and steady-state cosmic radiation conditions are presented. Steady-state biological 
intensities are given both for low-altitude orbits as functions of orbital inclination to the plane of 
the earth’s equator, and for paths lying far above the earth’s surface. 


Total allowable 


nuclear reactor radiations in the cosmic-ray background are nezt calculated as functions of satellite 


orbital inclination and exposure lime. 


These reactor radiation intensities give no 


more than the recommended maximum doses as stipulated by the U.S. Committee on 
Radiation Protection. A method of oblaining minimum weight reactor shields also is outlined, 
and a hand computation performed for a nuclear 300-kw electrical system. 


The harmful effects of radiation exposure were 
realized only after a number of injuries were observed 
among early workers in the fields of radiation and 
nuclear physics and others who worked with radio- 
active materials. In 1929, the Advisory Committee 


on X-Ray and Radium Protection was formed under ~ 


the sponsorship of the National Bureau of Standards 
to set safety regulations. In December 1946, with 
the advent of the atomic age, the radiation protection 
problems had become so numerous 
that the Committee had to enlarge 
its scope and membership. At 


that time, the name was changed : 


Fig. 1. Allowable average radiation intensity as a function of exposure time. 


equal to 83.9 ergs absorbed by each gram of air. 
When the term roentgen is applied to particles of 
energy greatly above 3 mev, as in the case of cos- 
mic rays, it is commonly taken to be a measure of 
their total ionization. The rad, quite simply, is the 
ionization energy absorption of 100 ergs in one gram 

of any material. 
The rem is a biological unit of ionization measure- 
(Continued on page 77) 


to the National Committee on Ra- . 

diation Protection (NCRP) and ten . at 

subcommittees were formed. Sub- 

committee 1 was responsible for in- . Dad 

vestigating the permissible doses 3 N 

from external radiation sources. At ae 

the meeting of the International Yearly Doge Notirepeated ‘Emergency Dosel of (25 [rem 
Commission on Radiological Pro- 


tection in 1953, the recommenda- 


tions of Subcommittee 1 were dis- 


cussed and its general principles 
adopted. Since that time some re- 


visions have been made. 


Common physical units of ioniza- 


tion measurement are the roentgen 
and the rad. Strictly speaking, the 
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inflated Satellite (Continued from page 11) 


sented in Appendix II, expressing the 
conservation of energy is used; the 
equilibrium condition requires that the 
radiation leaving the sphere equals the 
incoming direct solar radiation, earth 
radiation, and solar radiation reflected 
from the earth. 

The sunlight equilibrium temperature 
was found to be 394°K in the manner 
outlined in Appendix II. The values 
used in Eq. II-2 were obtained from 
reference 4 or the best available tables. 
The equilibrium temperature in the 
darkness was found to be 168°K. 

It is important to know the time re- 
quired to go from the equilibrium condi- 
tions in sunlight to the equilibrium 
conditions in darkness. Accordingly, 
the energy Eq. II-1 presented in Ap- 
pendix II was integrated in the manner 
shown in Appendix III. The time calcu- 
lated for the dark side gives a conserva- 
tive estimate of the time required to 
change conditions from darkness to 
sunlight because the effects of solar 
radiation are not involved on the dark 
side. It was found that darkness equi- 
librium conditions are established ap- 
proximately 6 min. after the sphere 
leaves the sunlight. Thus, equilibrium 
conditions prevail for the major part of 
the time that the sphere is in darkness 
and in sunlight. 

The internal equilibrium pressures 
corresponding to the equilibrium tem- 
peratures determined in the foregoing 
are computed analytically from Clapey- 
ron’s equation 


P)/d(1/T) = — AQsuvi/® (4) 


where AQwupi is the heat of sublimation 
and ® is the universal gas constant. 
The slope of the In P against (1/T) 
curve is determined from data taken 
from the ‘“Handbook of Chemistry and 
Physics” for a given material, and the 
equilibrium pressure is calculated by 
means of the point-slope formula 
anthraquinone is the subliming material 
specifically considered herein. 

If it is assumed that the anthra- 
quinone sublimes at the skin equilib- 
rium temperature, conservatively low 
values will be obtained for the vapor 
pressure on the dark side, as well as a 
maximum pressure to compare with the 
bursting pressure on the sunlight side. 


At the sunlight equilibrium tempera- 
ture (394°K) the internal equilibrium 
pressure is found to be 9.6 X 107*-mm 
Hg. The subliming solid equilibrium 
pressure corresponding to the darkness 
equilibrium temperature of 168°K is 
7.06 X 107~7!-mm Hg. 

Summarized in Table 2 are the signifi- 
cant results obtained under the as- 
sumption that the anthraquinone tem- 
perature instantly follows the skin 
temperature. 

From this table it can be seen that the 
internal pressure exceeds the maximum 
external pressure in sunlight and yet is 
sufficiently below the bursting pressure. 
In the darkness, however, the internal 
pressure is below the stagnation pres- 
sure and below the earth radiation pres- 
sure; this raises the question whether 
the sphere will collapse under the 
external pressure. To answer this ques- 
tion, a calculation is made in Appendix 
IV to see how much the sphere would be 
deflected by the external pressure in the 
time it would be in darkness (approxi- 
mately 25 min). The deflection was 
found to be 77 ft. This value is proba- 
bly fictitious, since no restraining forces 
were considered in the calculations, and 
a computation of the pressure difference 
required for buckling a sphere in Ap- 
pendix IV shows that 0.3 X 1074mm 
Hg is required, which is several orders 
of magnitude greater than that en- 
countered at 1,000 miles. Thus, it ap- 
pears likely that no deflection at all will 
occur. 

When conditions are such that the 
satellite surface strength will not with- 
stand the external pressure, as, for 
example, in the cases of a lower orbit 
or weaker skin material, it would be 
necessary to maintain the internal pres- 
sure above the external throughout the 
entire orbit. This problem is not neces- 
sarily resolved by simply utilizing a ma- 
terial of higher vapor pressure, because 
in the sunlight the bursting pressure 
of the skin might be exceeded. One 
possibility for moderating the tempera- 
ture extremes at which the solid sub- 
limes is to make use of the heat capa- 
city of the subliming solid. If the sub- 
liming solid is exchanging heat with the 
satellite skin by radiation, the vapor 
pressure variation will actually differ 


from the foregoing estimate which as- 
sumed that the temperature of the solid 
corresponds to the skin equilibrium tem- 
perature. If this assumption is not 
made, then the heat which is absorbed 
by the solid in sunlight, because of the 
heat capacity of the solid, will cause it 
to sublime at a higher temperature than 
the skin equilibrium value in darkness, 
In a like manner the sunlight sublima- 
tion will occur at lower than the sun- 
light skin equilibrium temperature. In 
Appendix V expression for the surface 
temperature (at which sublimation 
takes place) of a solid sphere are derived 
as a function of time. Although this 
configuration is arbitrary, calculations 
made from these expressions for 20 lbs 
of anthraquinone are compared in Fig. 1 
with those previously determined here. 
For convenience in this calculation, the 
sunlight and darkness times are taken 
to be equal (45 mineach). It is thus ap- 
parent that the vapor pressure is 
moderated from the extreme predicted 
earlier, and, furthermore, that the new 
extremes occur for only a fraction of 
each satellite period. The highest 
vapor pressure is approximately 1.34 X 
10°-mm Hg and the lowest about 
10-!!-mm Hg. Comparing these values 
with the values given in Table II, we 
see that the minimum internal pressure 
still goes below the maximum external 
pressure on the dark side, and, for a 
short time in sunlight the internal 
pressure (107!!-mm Hg) is lower than 
the solar radiation pressure (3.375 X 
10-*-mm Hg) acting on the outside of 
the sphere surface. Although this 
would still indicate for the present case 
the danger of collapse on these con- 
siderations alone, a comparison with the 
previous maximum and minimum in- 
ternal pressures shows the highest pres- 
sure to be reduced 3 decades, and the 
lowest to be increased over 10 decades, 
by making use of the heat capacity of 
the solid. For the present problem (100- 
ft-diam. sphere), there is no need for 
concern about the actual collapse so 
long as the sphere can withstand a 
buckling pressure greater than the 
radiation pressure. 

The decrease in solid mass with time 
is not included in these calculations. 
The amount lost will be only that 
vaporized from the solid at its maximum 
temperature plus that which sublimes 
during the dark period. This vapor con- 
denses on the cold satellite skin in the 


Table 1. Stagnation Pressure as a Table 2. Summary Table 
Function of Satellite Altitude Sunlight Darkuess 
Stagnation Equilibrium temperature, °K 394 168 
Height, pressure, Equilibrium internal pressure, mm Hg 9.6 X 1073 1.06 xX 10-** 
miles mm Hg External pressure, mm Hg 93.375 X 10-8 3.675 X 10~-° 
Bursting pressure, mm H 0.86 —e 
186 2.139 X 103 
500 3.827 X 1078 * The maximum external pressure in sunlight is due to solar radiation. 
900 5.028 XK 107 > The radiation pressure from the Earth is the largest external force in darkness fora 
1,000 2.526 X 107% 1,000-mile orbit. 
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One of a series of advertisements showing the breadth of Douglas experience being devoted to commercial, military and space-exploration projects. 
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into versatile Navy Skyhawk 


The new Douglas A4D-5 Skyhawk — a powerful ace-of-all- 
trades for the Navy — is the result of the same imaginative 
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commercial aviation and the push into outer space. 
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Sketch A. 


darkness, but it will re-evaporate from 
the warm skin when the satellite enters 
the sunlight, and will condense on the 
(relatively) cold solid eliminating further 
material loss from the solid. The present 
configuration was chosen only for the 
calculation and to show, easily, the ef- 
fect of the heat capacity of the solid on 
internal pressure. Other configurations 
would probably be chosen for an actual 
case. 


It might be worth while to use another 
material (rather than anthraquinone) 
with a higher vapor pressure and by 
utilizing the heat capacity of the solid 
keep the internal pressure from bursting 
the satellite in the sunlight. In this 
way, an internal pressure exceeding the 
external pressure could be maintained 
throughout the entire orbit. Extended 
use could be made of the heat capacity 
of solids to alter the internal pressure 
of inflated satellites in other ways, as, 
for example, by using a material with 
higher heat capacity and surrounding it 
with the subliming material. 


Miscellaneous Considerations 


In this section the possible influences 
of electrostatic, electrohydrodynamic, 
and electromagnetic forces on an in- 
flated satellite will be considered. 

It has been reported that satellites 
become charged with from 10 to 60 
volts as they orbit. Such an electro- 
static charge would cause a uniform 
outward force to act all around the sur- 
face of the sphere. Calculations (not 
presented herein) were made to deter- 
mine whether this force would be of the 
same order as the external pressure so 
that it could counteract it in the dark- 
ness to eliminate the positive external 
pressure. It was found that for 10 
volts the electrostatic pressure would be 
2.86 X 107'4-mm Hg, which is negli- 
gible. To balance an external pressure 
of approximately 10~"-mm Hg, it was 
determined that a potential of 1,259 
volts would be required. It appears 
unlikely that such a voltage would be 
obtained. Thus, this force is most 
likely insufficient to be either beneficial 
or detrimental to a 100-ft-diam. spheri- 
cal satellite. 


The electrohydrodynamic force is that 
resulting from the interaction of a con- 
ductor (the satellite) with the iono- 
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sphere. The electromagnetic force is 
that caused by the interaction of a 
conductor with the Earth’s magnetic 
field. No detailed calculations of these 
forces were made because they are es- 
sentially drag forces acting uniformly 
over the surface of the sphere (if the 
sphere is not so large as to encounter 
large field variations) and, thus, will 
neither deflect nor collapse the sphere 
surface. 


Although the life of the satellite in 
orbit would be difficult to calculate 
exactly, it is clear that all drag forces 
are low enough to permit the satellite 
to stay in orbit for years or until a 
catastrophic collision occurs. All col- 
lisions, even those puncturing the satel- 
lite, may not end its usefulness, and 
thus it is probable that once in orbit 
the satellite can be expected to remain 
there usefully for almost as long as a 
rigid satellite of the same size. 


Conclusions 


On the basis of the study presented 
herein it is concluded that a 1/2-mil 
thick Mylar satellite 100 ft in diameter 
inflated with anthraquinone will retain 
its shape during the entire course of an 
orbit at 1,000 miles altitude. This 
conclusion is based on estimates of the 
buckling strength of the sphere which 
appears to be several orders of magni- 
tude greater than the deformation 
forces at an altitude of 1,000 miles. 
Thus, the internal pressure is only 
necessary to inflate the sphere, but plays 
no important role in keeping it inflated. 

Under other conditions, however, 
when the skin strength is insufficient to 
withstand the deformation forces, the 
satellite shape and performance will be 
determined by the internal pressure. 
Since this pressure is very temperature- 
dependent and since the satellite is 
subject to large temperature variations 
as it passes through darkness and sun- 
light portions of its orbit, it is extremely 
difficult to keep the internal pres- 
sure below the bursting pressure in the 
sunlight and above the external pressure 
in the darkness if the vaporization oc- 
curs at the skin equilibrium temperature. 
A method of obviating this problem is to 
utilize the heat capacity of the sub- 
liming solid. Calculations presented 
herein show that the heat capacity 


causes the subliming solid temperature 
to lag the skin equilibrium temperature 
and, thus, the pressure extremes are 
moderated. It appears possible then to 
use a solid with a high vapor pressure 
to provide an excess of internal pressure 
always and still stay within bounds of 
the bursting limits. 


Appendix |. Radiation Pressure 


The radiation pressure is defined in 
reference 6 as 


Praa = (1/c) I cos? @ (I-1) 


where J is radiation energy permit time 
per unit area perpendicular to the direc. 
tion of propagation of J, c is the velocity 
of light, and 6 is the angle the direction 
of propagation of radiation makes with 
the normal to the surface under con- 
sideration. From this definition, Eq, 
I-1 gives for solar radiation pres- 
sure on a black disk perpendicular to the 
solar radiation 


Praal oAgunT (1-2) 


where R is the distance between earth 
and sun, Asunis total area of the sun, ¢ 
is the Stefan-Boltzmann constant, and 
Teun is the temperature of the sun. To 
obtain Eq. I-2, it was assumed that the 
sun is a point source and a black body. 
In calculating the average collapsing 
pressure on a reflecting sphere, we ob- 
tain a factor of 1/2 for the average value 
of the cos? @ in Eq. I-1 (6 is the angle 
around the sphere) and a factor of ap- 
proximately 2 to account for the reflec- 
tion. These two factors just cancel and 
the value of the average collapsing pres- 
sure due to solar radiation shown in the 
paragraph following that of Eq. 3 may 
then be calculated directly from Eq. 
1-2. 

For the radiation pressure due to the 
earth, Eq. I-1 leads to 


Praay = (1-3) 


for the radiation pressure on a small 
black disk perpendicular to the line 
joining earth and disk, where ez, Tz, 
and A, are the emissivity, temperature, 
and area of the earth, respectively, and 
Fy, is a coupling factor between the 
visible cap of the earth and the satel- 
lite. We may define Fz_, (sketch A) as 


AgFs-, = 


cos 6; cos? dA dA sat 
Asat Ag aR? 
Re 2 
I 4) 


where h is the height of the satellite 
above earth, and is the earth's 
radius. 


Combining Egs. I-3 and I-4, we have 


| 
| 
we 
Be 
Te 
i 


ApriiryOr Aerospace engmecning we 


‘ature 
ature 
S are 


-ssure 


“ssure 


| | the railroad” 


1ed in 


= Hard basing is one way to protect 
(1-1) America’s force of retaliatory 
RE . : = ICBM’s. The problem was to find an 
t tine alternate means of accomplishing 
direc. the same mission. The Air Force 
locity solution was a new ICBM mobility 
ection concept—railroad car-mounted Min- 
s with utemen, utilizing the nation’s vast 
con. track mileage for numerical and 
pres- difficult target for enemy attack. 
tote To put the Minuteman, its support 
systems and associated equipment 
on rails was a completely new prob- 
(1-2) lem in missile handling. The first 
requirement assigned by Boeing to 
t, and a feasibility study of the existing 
— limitations of roadbeds, rails, rail- 
body. Unique tactical cars are being de- 
— signed within these limitations to 
we ob- carry the Minuteman—cars that can 
evalue handle the missile and its operating 
om. roadbed shock and ready for immedi- 
3 may Whether for conceptual problems 
m Eq. such as this one, or for challenges in 
design or manufacturing, AMF has 
-to the ingenuity you can use. AMF people 
are organized in a single operational 
unit offering a wide range of engi- 
Its purpose—to accept assignments 
sm — at any stage from concept through 
he link development, to production, and 
ex, Ta j : service training...and to complete 
erature, * them faster in 
ly, and * Ground Support Equipment 
een the * Weapon Systems 
e satel- Undersea Warfare 
A) as Radar 
* Automatic Handling & Processing 
Range Instrumentation 
* Space Environment Equipment 
¢« Nuclear Research & Development 


GOVERNMENT PRODUCTS GROUP, 


AMF Building, 261 Madison Avenue, 
New York 16, N. Y. 


) (14) 


satellite 
earth’s 


we have 


Seectine and manutacturing AMF has ingenuity you can Use...» AMERICAN MACHINE & FOUNDRY COMPANY 


ta 
ure 
| 
te 


Assembled at the Columbus Division of 
North American Aviation are the facilities 
and the proven technical intellect to bring 
original concepts swiftly to practical produc- 
tion by the most economical and efficient 
methods. Here, in one of the most complete 
centers of advanced systems technology in 
the world, many of the important advances 
in electronic, electromechanical, and envi- 
ronment systems, as well as other areas, have 
been made. This is true systems capability 
... this is the Columbus Division. 


THE COLUMBUS DIVISION OF 2: 


wind 


NORTH AMERICAN AVIATION, INC. 


Columbus, Ohio 


for the study of transition from vertical to level flight 


is a part of the facilities at Columbus. Other V/STOM 
equipment includes six-degree-of-freedom flight simu 
lator, and zero-altitude, zero-speed escape systems 


HIGH OR LOW. The Mach 2 A3J Vigilante, now being 
built for the Navy at Columbus, operates effectively at 
deck level or up in the stratosphere. The versatile A3J 
can perform either attack or interception missions in 
any weather, at any attitude or altitude, day or night. 


34 Aerospace Engineering «+ April 1961 


“Gq” sEAT. Columbus developed unique “G” seat to studyge 
human tolerance to a high degree of vertical accelera™ 
tion, and to mate human factors to the machine. Ad Fé 
vanced low-level high-speed escape systems and pilot” 
response in these flight regimes will be studied with it. 


fi 


tunnel 
flight a 
/STOL 
simu 
stems 


DIO/RADAR ANTENNAS. An. advanced type of antenna, constructed to such 
o studyfetise tolerances that it will be tuned to exact’ pitch, will be built by NAA- 
ecelera- blumbus for the Air Force on “Haystack Hill” near Tyngsboro, Mass. Unique 
ne. Adpsign developed at Columbus allows the rigid specifications to be met with 
id pilotfnificant reductions in dead weight, and in control power requirements. 
with it 


April 1961 + Aerospace Engineering 35 


Re 2 
c Re + h 


ecoA 
(1-5) 
4ac(Re + h)? 
Since this is the pressure on a small 
black disk, we would again have a fac- 
tor of 2 for a reflector and a factor of 
approximately 1/2 for the reduction in 
the average pressure around a sphere. 
Thus, to calculate the average collapsing 
pressure on a reflecting sphere, we may 
use Eq. I-5 directly, and obtain the re- 
sults presented in Table 2. A compari- 
son of the values of radiation pressure 
indicates the solar radiation pressure to 
be approximately 9 times that due to the 
earth at a height above the earth of 
1,000 miles. 


Appendix Il. Conservation of 
Energy Equation and Equilibrium 
Temperature Calculations 

The equation expressing conservation 
of energy for a satellite in sunlight can 
be written as 
dT sp 

dt 


= ass(Bs + 


0.8 Bez) —4éspoTs,(t) (11-1) 


Bg = heat flux from the earth 


B, = heat flux from the sun 

B,z = heat flux from the sun reflected 
by the earth 

b == sphere wall thickness 

Cs = specific heat of Mylar 

Tsp = sphere temperature 

= time 


ase = absorptivity of the sphere for 
earth (low-temperature) ra- 
diation 

ss = absorptivity of the sphere for - 
solar (high-temperature) ra- 


diation 
&p = emissivity of the sphere 
ps = Mylar density 


The following assumptions are implicit 
in Eq. II-1: 

(1) The temperature of the satellite 
surface is uniform. 

(2) The angle factor between the 
earth and sphere? is 0.8. 

Under equilibrium conditions, the 
sphere must have a constant energy. 
Therefore, the left side of Eq. II-1 must 
be set equal to zero in order to compute 
equilibrium temperatures; this yields 


1 


ass (Bs + (08)B.») | (11-2) 


(a) Sunlight equilibrium temperature. 
The following values were used in Eq. 
TI-2 to determine the sunlight equilib- 
rium temperature: 


4 


= = 


0.0534 cal/(sec) 


0.0333 cal/(sec) (cm?) 

= 0.006 cal/(sec) (cm?) 

300°K 

= 0.1 

0.00127 cm 

1.4 g/cm* 

0.3 cal/(g) (°K) 

1.55 X 107! cal/(sec) (cm?) 
(°K*) 

= €,» = 0.03 (assumed gray be- 

cause of low temperature) 

Under these conditions the sphere 

equilibrium temperature, 7,,,. is 394°K. 
(b) Darkness equilibrium temperature. 

To compute the equilibrium tempera- 

ture in the darkness, all the afore- 

mentioned values were used in Eq. 

II-2 except that B, = Bye = 0; then 

Tip, = 168°K. 


Appendix Ill. Time Required for 
Equilibrium Conditions to Be 
Established 


Equation II-1 is made dimensionless 
by substituting 


4a,oT,* 
Te psbcs 


into it, where 7; is the earth’s tempera- 
ture. The equation can be written as 


dn/dr = — 


where B, = B,z = 0 on the dark side. 
The solution to Eq. III-1 using the data 
in Appendix IT is 


0.54 
E 


0.54 — » 
3.18 | 


(III-1) 


7 

(III-2) 
where = and T,, is the sun- 
light equilibrium temperature. For the 
temperature to become 168°K (the 
darkness equilibrium temperature, 7 = 
0.56, 7 = 3.788, and t = 6.76 minutes. 


Appendix IV. Sphere Deflection 
and Buckling due to the Net 
External Pressure 


(1) The maximum deflection of the 
sphere during the darkness period, when 
the external pressure exceeds the inter- 
nal pressure, can be estimated from 


d?x/dt? = P/(m/A) (IV-1) 


where restraining forces and the internal 
pressure are neglected, and m/A is the 
mass per unit area of the sphere; x is 
the deflection, ¢ the time, and P the ex- 
ternal pressure. Eq. IV-1 gives the 
maximum acceleration of a section of 
the sphere relative to the rest of the 
sphere in the absence of restraining 
forces. Integration of Eq. IV-1 with 
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the initial conditions dx/dt = x = 0 
at t = 0 yields 


x = [P/2(m/A)]t (IV-2) 


At an altitude of 1,000 miles and a mass 
per unit area of 1.75 X 107% g/cm?, the 
deflection is 77 ft after 25 min in the 
darkness. 

(2) The critical buckling stress, o.,, 
of the sphere can be estimated from the 
relation 


Gor = Eb/RV/3(1 — v) (IV-3) 


taken from reference 8, where E is 
Young’s modulus, } is the shell thickness, 
R the sphere radius, and » Poisson’s 
ratio. An order of magnitude for the 
pressure difference for buckling, P.,, 
is then determined from 


Pe 


to be 0.3 X 10°-*mm Hg. The actual 
buckling pressure may be as low as 1/4 
this value because of wrinkles in the 
skin of the sphere.* 


Appendix V. Heat-Sink 
Temperature 


Consider a sphere of solid anthra- 
quinone of radius a (to weigh 20 Ibs.). 


At time zero (say as it enters sunlight - 


from darkness) the solid has the tem- 
perature 7;“ and is subjected to the 
environment temperature 7, (sunlight 
equilibrium temperature). The conduc- 
tion equation 


2ar\ _ ar 
“\ Or? ror) a 


and boundary conditions 


(V-1) 


oT 
—k— =h(T—T\)atr=a_ (V-2) 
or 


together with the initial conditions will 
give the temperature of the solid. The 
thermal diffusivity is denoted by a, k 
is the thermal conductivity, and fp is a 
measure of the heat losses (due to radia- 
tion) from the solid. Let 


v= — T, and h = ra (V-3) 
Eqs. V-1 and V-2 become 

ov 

— = — V-1 
o( Ore r =) ot 
dv 
+hv=Oatr=a_ (V-2a) 
or 


The solution® at r = a to this problem 
is given as 


a’B,2 + (ah — 1) 
Bn2[a7Bn? + ah(ah — 1)] 


sin?aB, (V-4) 


where £, is obtained from 


aB, cot a8, +ah —1=0 (V-5) 


As the sphere enters darkness (say 
t = 0 again, arbitrarily), the environ- 
ment temperature is 7, (darkness equi- 
librium temperature) and the solid tem- 
perature is Then 

T — A, 
= 2- —aBnrt 

a? B,? + (ah — 1)? 
+ ah(ah — 


Now consider 


sin*ag, (V-6) 


a =115cm 
T,; = 394°K 
T2 = 168°K 


0.113 cal/(sec) (em?) (°K) 


k 10~8 cal/(sec) (em) (°K) (guess 
like naphthalene) 

c = 0.258 cal/(g) (°K) 

p = 1.419g/cm* 

a = 6 X 107% cm?/sec 


Eq. (V-5), therefore, becomes 
11.58, cot 11.58, + 0.5 = 0 
which has the solutions?® 


11.56, = 1.8366 11.58; = 7.9171 


11.562 = 4.8158 11.58, = 11.0409 


Substituting these values and 2.7 xX 
10° sec (45 min) into Eq. (V-4) yields 


T = 394 + 0.441(7; — 394) = 
T;® (V-7) 


and from Eq. (V-6) 


T = 168 + 0.441(7; — 168) = 7; 
(V-8) 


Simultaneous solution of Eqs. (V-7) and 
(V-8) gives 


= 937°K T; = 325°K. 


and now the complete solid temperature 
variation with time can be determined 
from Eggs. (V-4) and (V-6). The qualita- 
tive behavior is shown in Fig. 1. 
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jualita- Said J. Stefan and L. Boltzmann: “The total radiation from a black body is proportional 
to the fourth power of the absolute temperature of the black body.” 


Radiation is usually associated with high temperatures. Yet very cold bodies emit a radiation which can be highly 


Pre. significant in missile and space applications. The problem faced by infrared scientists, trying to detect variations in radia- 
nflatable tion from low temperature atmospheres, can be likened to detecting a one-foot cube of ice from a distance of five miles. 
h Satel- 


Lockheed Missiles and Space Division scientists are deeply engaged in studying the problems of infrared emission 


siete from the earth and its atmosphere, as seen from orbital altitudes. Although the earth resembles a black body at 300° Kelvin, 
re the emission from its atmosphere, under some circumstances, is much colder. To make measurements under these cir- 
: 3 cumstances, Lockheed has evolved radiometric equipment with one of the most sensitive detection systems yet conceived. 
.S. W, Scientists and engineers must also take careful measurements of a potential employer. Lockheed Missiles and Space 
| Atmos- Division in Sunnyvale and Palo Alto, California, on the beautiful San Francisco Peninsula, invites this close scrutiny. As 
rics, No. Systems Manager for the DISCOVERER and MIDAS satellites and the POLARIS FBM, Lockheed preeminence in Missiles and 
9. , Space creates positions in many disciplines for outstanding engineers and scientists. 

S., an 

ind Co, Why not investigate future possibilities at Lockheed? Write Research and Development Staff, Dept. M-13A,962 West 


El Camino Real, Sunnyvale, Calif. U.S. citizenship or existing Department of Defense industrial security clearance required. 
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Launching and Recovery (continued from page 13) 


beam, brake pucks free of the top and 
bottom surfaces of the flange, until the 
car comes to a predetermined point in 
the tracks. Here a finger-switch trips a 
fuel cutoff valve shutting down the jet 
engines, then the flange of the J beam 
thickens so that the brake pucks wedge 
against the top and bottom surfaces, 
grinding the jet car to a stop. The 
deadload, which is pushed along in 
front, goes sailing into the gear at what- 
ever speed has been determined by the 
rpm setting of the jet engines. The jet 
car and the deadload ride on rubber- 
tired wheels rather than directly on the 
track J beams, thus tight tolerance is 
not required for track alignment nor do 
the joints have to be exactly mated. 
This is fine for moderate speeds, but 
what will happen to the tires beyond 250 
knots is not yet clear. Firestone Rubber 
Co., who makes most American racing 
tires, may be able to help. After the jet 
car is brought to a stop, the deadload is 
engaged by a cross-deck pendant. For 
ease in catching the wire, the hook is 
usually in an up-raised position, in- 
verted from the angle at which it nor- 
mally trails from an aircraft. 

(2) Runway Arrested Landing Site. 
This site houses the latest shipboard- 
type arresting engines, as well as older 
ones on which service problems might 


(center). Both catapults are d 


(Official photograph, U.S. Navy.) 


38 Aerospace Engineering 


April 1961 


arise. Actual carrier layout is dupli- 
cated here—i.e., sheave locations, deck 
pendant length, purchase cable length, 
terminal impact pads, engine location 
and arrangement. Test engagements 
with the latest service aircraft are made 
under controlled conditions, completely 
instrumented. Off to the side there is 
an expeditionary control tower and in- 
strumentation center for recording data. 
An SPN-12 radar, similar to shipboard 
equipment, is used here to control the 
approach speeds of aircraft and deter- 
mine their engaging velocities. Speeds 
are checked with police radar which is 
fairly close even at 100-120 knots, giv- 
ing accuracies within 2-7 percent of 
other readings. An infrared velocity- 
sensing device has been provided by the 
Naval Air Development Center, Johns- 
ville, and extensive photographic cover- 
age is utilized. Located at this site, 
also, is an optical landing-aid mirror 
with the light arrays outstretched, in 
the center of which the pilot sees his re- 
flected image, called the ‘‘meatball,”’ 
when simulating preset carrier-glide 
path approaches. Emergency overrun 
gears are placed at other points along 
the runway for test. It is only after all 
components have been thoroughly 
checked out at the Recovery Systems 
Track Sites, with deadloads first, and 


Fig. 3. Pant (at right) supplies power for flush-deck catapult TC14 (left) and elevated-deck TC13 

igned to launch high-performance fleet aircraft. TC13 simulates air- 
craft carrier environment in providing elevation to study aircraft sink effect after leaving deck-edge. 
Beveled ramp at forward end of TC13 facilitates launchings with deadload, seen in battery position. 


then with stricken aircraft hulks, that 
they are brought over to this site for 
evaluation with live aircraft. 

(3) Catapult Site. Two high-energy 
shipboard-purpose catapults (Fig. 3) 
are installed—one steam-powered, the 
other configured to operate on either 
compressed air or the gases from an in- 
ternal combustion powerplant. The 
former is an elevated installation that 
provides means to conduct development 
and evaluation programs of the latest 
catapult equipment, TC-13, yet also 
permits simulation of the rotational be- 
havior of different aircraft types a they 
are launched from the deck-edge of an 
aircraft carrier. The latter, designated 
TC14, is a flush-deck installation. Both 
catapults have a 250-ft power stroke 
and use the same types of launching 
cylinders, shuttles, water brakes, re- 
traction engines, and other components 
as practicable. The TC13, of NAEF 
(SI) design, is an improved and more 
powerful version of the steam catapults 
which are now in fleet carriers (provid- 
ing 63 million ft-lbs of energy versus 
present fleet capability of 42 million 
ft-lbs). It will go into the USS. 
“Kitty Hawk,” “Constellation,” ‘En- 
terprise,’” and CVA-66. The TC14 
originally included an _ experimental 
powerplant design by Reaction Motor 
Division, Thiokol Chemical Corp., which 
had been under development in the 
event that future aircraft weights and 
launching speeds would require greater 
energies than were available from steam. 
It now operates on compressed air. It 
is a useful tool for research, and for 
component evaluation, when the TC13 
is so often committed to higher priority 
programs. Any limited end-speed or 
off-center launches that would be dan- 
gerous on the elevated platform are also 
conducted from the flush-deck site. It 
is planned to modify this to an extended 
stroke (350 ft) capability in 1963. 

The TC13 catapult incorporates fea- 
tures which give it higher reliability, 
longer service life, and rapid-launching 
cycle capability of approximately 30 
sec for even the heavier weight aircraft 
approximating 70-80,000 Ibs. The 
TC13 derives its launching energy from 
high-pressure steam. The accumulators 
are charged before making each launch, 
the rate of steam flow from these to 


‘the launching tubes being controlled 


for each aircraft weight and desired 
end speed, by the launch valves (Fig. 4). 
A two-stage-orifice hydraulic system 
controls the opening rate of the launch 
valves and provides the means for limit- 
ing the jolt factors; that is, the rate of 
acceleration rise (10g/sec normal, 16g/ 
sec max.), and the sustained accelera- 
tion during a launch (3 to 4g normal, 
6g max.). The aircraft is attached to 
the shuttle by a wire-rope bridle or 
pendant. When the catapult is fired, 
high-pressure steam enters the cylinders 
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hat behind the pistons, the shuttle assembly 
for is driven forward towing the airplane 
with it. Near the end of the piston 
rgy stroke, the steam is vented through ex- 
3) haust valves and the piston assembly 
the is brought to a stop by a water brake. 
her The entire assembly is stopped in ap- 
in- proximately five feet, experiencing a 
lhe deceleration of 200 to 300 G’s. The re- 
hat traction engine is used to return the 
ent shuttle to battery position after each 
test launch. 
also As of December 31, 1960, more than 
be- 7,000 launchings with deadloads and 
hey aircraft had been made on this catapult 
an (Fig. 5). This large number is signifi- 
ited cant in that it represents the greatest 
oth number of launchings ever made on a 
‘oke new catapult before fleet use. Three 
ling thousand additional launchings are 
re- scheduled to be made before mid-1961. Fig. 4. Cutaway of a steam-type catapult illustrating: (1) steam accumulators, (2) launching valves, 
ents This large number of launchings is (3) launching pistons and shuttle, (4) water brakes, (5) exhaust valve, (6) retraction engine, (7) main 
EF indicative of the degree of effort that is control console. (Official photograph, U.S. Navy.) 
nore being exerted to assure high reliability, 
ults easy maintenance, and long service life 
vid- for our fleet catapults. 
lion 
Ss. Expeditionary Catapults 
‘En- * Short-field and expeditionary launch- 
C14 ing systems are also under development 
ntal and evaluation by the Navy, primarily 
otor for U.S. Marine Corps requirements. 
hich The purpose of these is to provide 
the means for launching high-performance, 
and close-support, or air-cover jet aircraft 
ater at advance bases where runway lengths 
eam. may be limited. They are charac- 
It terized by simplicity in design, installa- 
| for tion, and operation; all-weather capa- 
C13 bility; suitability for various terrains; 
ority low packaged weight and minimum 
d or logistic support requirements. Three 
dan- typical expeditionary catapults are the 
also M-1B Turbo-Cat, the XRE-1 Cata- Fig. 5. The Navy's latest supersonic attack bomber, North American Aviation Vigilante, being . 
It | port, and the Jet Car Catapult. Other launched from the U.S. Naval Air Test Facility (SI) TC13 steam catapult. (Official photograph, U.S. 
nded designs are being considered as well. Navy.) 
fea- 
ility, M-1B Turbo-Cat 
hing The All American Engineering Com- 
y 30 pany M-1B Turbo-Cat is basically a re) 1000 FEET 
craft 50,000-hp gas turbine, driven by the ex- k. 
The haust of six conventional jet engines. = 
from The main components of the system 
ators are six J33 jet engines, turbine wheel, 
inch, flow gates, nozzle vanes, capstan, brakes, 
se to cable, sheaves, shuttle, and track. 
-olled The turbine wheel is directly con- 
sired nected to the capstan which drives 1,800- ~~ 
g. 4). ft length of wire cable. A detachable _~ 
stem shuttle, connected to the towing cable, SS 
vunch transfers the accelerating force to the air- NYLON TAPE cenicgdlh™ ee 
limit- craft. The power stroke is approxi- 
ate of mately 600 ft. 
16g/ The profile of the capstan is so de- aa: 
elera- signed to prevent vertical motion of the \ 
rmal, cable during launch. The brake system 
ed to is designed as a parking brake and 
me "4 ppoping nena. The parking brake = Fig. 6. This is a sketch of the E. W. Bliss Co. cataport, a rotary-type catapult driven by a 134 
fired, used to prevent rotation of the turbine engine. Launching stroke is 1,000 ft Kinetic energy is 25 million ft-lbs. Note that the cataport is a 
inders wheel when the jet engines are being dual system—i.e., launching and arresting system. 
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Fig. 7. Take-off considerations for four-engine propeller and jet transports.* 


brought up to speed. The launch is 
initiated by release of the parking 
brake. The stopping brake is used to 
stop rotation of the turbine wheel after 
completion of the launch. The use of a 
capstan as the friction drive to ac- 
celerate the cable requires that some 
means must be provided to prevent the 
slack-side cable tension from dropping 
to zero. If this occurs, the driving power 
of the capstan is lost. To prevent this, 
a pretension system is used in the M-1B 
catapult. 


Approximately 100 deadload launch- 
ings have been made with the M-1B 
catapult at the All American Engineer- 
ing Company’s test site at Georgetown, 
Del. Deadload weights ranged from 
10,000-32,000 Ibs and speed ranged 
from 100-150 knots. A maximum 
energy output of 30 million ft-lbs was 
attained during the deadload test pro- 
gram. Incidentally, the power drive of 
this catapult has potential as a drive 
for a high capacity centrifuge, or as an 
emergency field power source for light- 
ing, pumping, etc. 


Cataport 


The E. W. Bliss Company “XRE-1 
Cataport”’ is basically similar to a cable- 
drive catapult powered by a T34 engine. 


The main components of the system 
are the T34 engine, clutch, tape reel, 
control console, shuttle, track, nylon 
tapes, sheaves, and anchors. The reel, 
nylon tapes, sheaves, and anchors are 
also used as part of the TE-15 arresting 
system. The launching stroke of the 
Cataport is 1,000 ft. The catapult 
launching energy is controlled by set- 
ting enginerpm. Complete performance 
data for this catapult is not available at 
the present time, however, results of 
initial tests show promise up to 25 
million ft-lbs. The Cataport is pres- 


ently undergoing development tests at 


the E. W. Bliss Company’s test site in 
Woodbine, N.J. Deadload develop- 
ment tests are expected to be completed 
around January 1961. If these tests 
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prove successful, a Cataport is to be in- 
stalled at NATF(SI) for evaluation with 
aircraft (Fig. 6). 


Jet Car Catapult 


This proposed All American Engineer- 
ing Company Jet Car Catapult will use 
a two-engine jet car and a 2 to 1 reeved 
cable system to provide necessary 
energy for launching aircraft. 


The main components of the system 
are a manned, steerable, jet car powered 
by two J-57 engines, shuttle, cable, 
sheaves, anchors, and cable tensioning 
system. Wheel brakes and a steering 
system will be provided for positioning 
of the jet car. Engine thrust reversers 
will be used to decelerate the catapult 
system after each launch and to return 
the system to battery position. The 
launching stroke for the catapult will 
be 2,000 ft. Available energy is ap- 
proximately 40 million ft-lbs. 


Adaptability to Commercial Use 


Until recently, only the Navy made 
extensive use of launching equipment 
in aircraft operations. The Navy’s 


*See Pilots Safety Exchange Bulletin 
59-105/106, Flight Safety Foundation, 
Inc., New York. (1) Required length 
must not be less than the take-off run 
(start of run. to 50 ft height reached) nor 
less than the accelerate stop distance 
(start of run to V; distance plus distance 
required to come to a stop). (2) Vii = 
maximum speed for which the aircraft 
can be safely braked to a stop or continue 
take-off even with one engine. (3) Vo = 
minimum safe take-off speed. Equals to 
1.15 X Vs (stall speed) for propeller air- 
craft and 1.2 X Vs stall speed for jets at 
zero thrust and idle engine power. (4) 
Vior = minimum lift-off speed attained 
prior to Vo. (5) Vp = speed at which air- 
craft is rotated in order to obtain desired 
angle of attack for maximum lift. 


need resulted from the limited operat- 
ing space aboard its carriers. The Air 
Force has done some experimental work 
with zero-length launchers for small 
field application directed toward NATO 
needs. VTOL/STOL aircraft have not 
come along as fast as predicted. Space 
problems face a great many activities, 
Not too many airports can be adapted to 
jet operations. To maintain any degree 
of safety, high-performance military 
and commercial jet aircraft require run- 
ways in excess of 9,000 ft. Special 
penetration or approach corridors and 
climb-out patterns are necessary be- 
cause of noise and safety considerations, 


Although the FAA and commercial 
airlines do not have, at this time, any 
projects under way for development of 
land-based catapults, planning for 
launching requirements is in the think- 
ing state. Commercial aircraft are 
structurally more limited in permissible 
acceleration loads than military air- 
craft. The FAA and commercial airlines 
are also more concerned with passenger 
comfort. The FAA’s present concepts 
relating to utilization of launching 
equipment are as follows: 


(1) for take-off assist only 

(2) 1,000-2,000 ft launching stroke 
(3) 1/2 G maximum acceleration 
(4) 100-200 million ft-lbs of energy 
(5) 1-min launch interval 


With development and adaptation of 
assisted take-off equipment, it would be 
possible for the FAA to reduce runway 
length requirements for jet transports 
while increasing the margin of opera- 
tional safety. Vi, he go-no-go decision 
or refusal speed, as it is called for take- 
off, could be achieved much earlier, al- 
lowing a great deal more runway ahead 
and psychological breathing room for 
the pilot (Fig. 7). For a no-emergency 
condition, rotation of the airplane at 
the optimum time would be less critical, 
and climb-out speed easier to attain. In 
case of engine failure or flame-out, 
starlings on the run-way, or some other 
emergency condition, distance for abort- 
ing take-off would be ample. Use of 
the assisted take-off technique would 
also be reflected in increased payload 
capability. 


Rockets or JATO of present design 
are too expensive for repeated opera- 
tions at any air field, and add to the 
nonpayload weight that must be carried 
by the airplane. A steam rocket- 
powered carriage or other low-operating 
cost device to push against the main 
landing gears, or other strength point, 
might be a solution. The water-jet 
principle, used in NACA landing car- 
riage tests at Langley Field some years 
ago, could be applied. Present Navy 
shipboard or land-based catapults are 
not suitable for immediate adaptability 
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to commercial airline use. A hybrid, 
possessing the high-reliability and rapid- 
cycle capability of the shipboard gear, 
the long-stroke and low-acceleration 
characteristics of the land-based gear, is 
required. It could be either nose wheel 
tow, as is being perfected for new car- 
rier aircraft, or pusher design. How- 
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ever, a wide range of aircraft launching 
weights and configurations must be ac- 
commodated by any method, and much 
more study seems necessary before we 
can define the mechanical specifica- 
tions for a universal airport launcher for 
the next generation of commercial air- 
craft. 
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The Interdependence (continued from page 15) 


R,(mg/em?) = 412E” (1) 


where = 1.265 — 0.0954 1n From 
~1 to ~20 mev 


R,(mg/cm?) = 5380E — 106 (2) 


where in each case E is in mev. 
Bremsstrahlung Radiation—High- 
energy electrons convert part of their 
energies to bremsstrahlung. An expres- 
sion which is approximately valid for 
electron energies up to about 2.5 mev is 


I/B ™ 7(10)-* ZE (3) 


where J is the amount of bremsstrahlung 
produced in mev per incident electron, E 
is the incident electron energy in mev, 
and Z is the atomic number of the target 
nucleus. Thus, the fraction of incident 
electron energy converted to brems- 
strahlung is about 1 percent for 0.5-mev 
electrons absorbed in copper. 

A more general expression for the 
ratio of radiation to ionization loss is 
given by 


(4T/ds)raa. Z (= y T 


My (4) 


(dT/ds)ion 1,600 moc? 
where Mp is the rest mass of any charged 
particle. For electrons, My = mo and 
the radiative and ionization losses are 
equal for a kinetic energy T = 20moc? = 
10 mev in lead (and for T ~ 100 mev in 
water or air.) 

Absorption of the Electron Energy— 
Fig. 1 shows the range-energy relation- 
ship for electrons absorbed by alumi- 
num. Relations (1) and (2) given 
earlier fit the curve very well in the 
energy range 10 kev to 20 mev. Only 
the absorption curve for aluminum is 
shown since the range for electrons is 
not strongly dependent upon the atomic 
number of the stopping material so 
that the curve given for aluminum can 
apply within limits for other materials 
as well. 

Fig. 2 shows the mass absorption 
energy losses for electrons in air, alumi- 
num, and lead. The losses due to both 
ionization and radiation are presented. 
This graph points out that the better 
Shield for high-energy electrons is a 
low-density material such as aluminum 


since a considerably smaller portion of 
the electron energy is converted to 
bremsstrahlung with its correspondingly 
higher penetration ability. Clearly, a 
wider range of materials needs to be 
dealt with and included in a more 
comprehensive study. 


Applications to an Orbit—Since we are 
specifically considering an orbit in the 
vicinity of 500-km above the earth’s 
surface, reference to the data furnished 
by AFSWC for electrons reveals that 
there exists a flux at this altitude of 
less than 10‘ electrons in excess of 30 
kev of which less than 10 range between 
460 kev and 1 mev. Reference to Fig. 1 
shows that 400 mg/cm? of aluminum 
will suffice to stop all the electrons. 
The amount of bremsstrahlung produced 
is insignificant. Consequently, for these 
data we can discount the effects of elec- 
trons as a radiation hazard for this 
type of vehicle orbiting at 500 km. 


Units—Mention should be made at 
this time about the units used. Above, 
it was stated that 400 mg/cm? of alumi- 
num will effectively stop all the elec- 
trons expected to be encountered. Also, 
it has been mentioned that, although 
only the curve for aluminum is pre- 
sented, the error resulting in using this 
curve for other materials is not exces- 
sive. Consequently, it could have been 
said that 0.4 g/cm? of steel, copper, or 
lead would also have sufficed. Thus, we 
have a consistent system of units which 
is easily translated into a required 
thickness of whatever material is to be 
used merely by dividing by the density 
of that material. For an example, see 


Table 2. 
Table 2. Electron Absorption by 
Several Elements 
Thickness 
Density Required 
Material (g/cm?) (cm) 
Magnesium 1.74 0.230 
Beryllium 1.85 0.216 
Aluminum 2.70 0.148 
Titanium 4.54 0.088 
Iron 7.87 0.051 
Nickel 8.90 0.045 
Lead 11.34 0.035 
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Proton Shielding 


General Discussion—High-energy pro- 
tons, because of their high penetrability 
of matter, present the major shielding 
problem, and therefore the following 
material is presented as an aid to devis- 
ing shields for protecting the occupants 
of space vehicles. Perhaps the greatest 
unknown, at present, is the density and 
energy of the protons normally to be 
encountered in any space flight. Type 
3+ solar flares, however, complicate 
the problem. Flares of this type con- 
tain protons in the hundreds and thou- 
sands of mev and require special shield- 
ing consideration. 

High-energy protons penetrate, with 
decreasing energy, into a considerable 
thickness of matter. This property of 
matter to slow down protons is referred 
to as the ‘stopping power” of the ma- 
terial traversed. The quantity —dE/dx 
is a function of energy. The energy loss 
by protons as they traverse matter is 
due primarily to ionization and excita- 
tion of the constituent atoms. If the 
material is thick enough, the proton 
will lose all its energy and so remain im- 
bedded. However, the proton may 
have sufficient energy to penetrate the 
target completely, in which case it will 
emerge with a portion of its original 
energy. 

Range-Energy Relationship — The 
theory of energy loss by energetic pro- 
tons through matter has been developed 
and verified by experiment. Simplify- 
ing assumptions in the theory lead to 
the approximate range-energy relation.! 


R = 100 (E/9.3)"8 (5) 


for protons between 10 and 200 mev, 
where R is the range of protons in centi- 
meters of air at 15°C and 760 mm Hg, 
and E is the proton energy in mev. ‘ 

Range-energy relations have been 
plotted in Fig. 3 for protons of energies 
between 1.0 mev and 10 bev. In addi- 
tion to the several elements presented, 
it would be desirable to have graphs 
available for the more common struc- 
tural materials such as type 347 stain- 
less steel, type Ti-4Al-3Mo-1V ti- 
tanium alloy, or type 2024 aluminum 
alloy, to mention but three of a multi- 
tude. 

Stopping Power—The “stopping” 
power” of a material as it is traversed 
by a charged particle is defined as the 
differential space rate of loss in energy 
of that particle in the substance. This 
concept is an especially important one 
in connection with the penetration of 
protons through the walls of satellites. 
Although it may not be possible to 
shield against all protons (since those 
of high energy will penetrate consider- 
able thicknesses of lead), these data 
plotted for the designer will enable him 
to determine readily whether a given 
thickness of wall material will reduce 
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PERFORATED METAL TUBE: 


The aerial refueling system on the U.S. Air Force KC-135 
jet tanker delivers fuel at an extremely high flow-rate—yet 
must disconnect instantly when required. 


To take the tremendous shock pressure when flow is 
interrupted, B.F.Goodrich designed and built for Boeing a 
double-walled surge boot to carry fuel in the flying boom 
extension. A pressurized air cushion between the walls of 
the boot dampens shock loads, permits fast refueling 
with safety. 


B.F.Goodrich designs and produces special products for 
a wide range of fueling requirements—surge boots, flexible 
fuel line connectors, fuel cells, hose, vapor barriers, and 
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SURGE BOOTS — absorb hydraulic shock. 


FUEL HANDLING 


B.F.GOODRICH PRODUCES THE ANSWERS 


INNER BOOT 


FLEXIBLE CONNECTORS— 
for internal fuel lines. 


others for aircraft... motor cases, insulator liners, and solid 


propellants for rockets and missiles. For help in your appli- ) 
cations check B.F.Goodrich Aviation Products, a division of | 


The B.F.Goodrich Company, Dept., AS-4, Akron, Ohio. 


PRESSURE SWITCHES—BFG 
switches actuated by precision, 
Omega-design bellows perform 
many critical functions on jet 
engines such as indicating 
adequate fuel pressure, actuat- 
ing afterburner igniters and 
unlocking nozzle area controls. 


h 
Sp 
Ca 
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EXPULSION DIAPHRAGM—BFG supplies a 
special flexible diaphragm mounted in the fuel 
cavity of the Bomarc supersonic area defense 
missile. Diaphragm expels hydrocarbon fuels 
under pressure to external ram-jet engines. 


WEIGHT-SAVING FLEXIBLE FUEL CELLS —This T-38 twin-jet trainer, 
made by Northrop, is one of many current aircraft using BFG fuel cells. Constructed 
with special fuel-resistant compounds, these cells are strong and lightweight. BFG 
fuel cell manufacturing is backed by complete development and test facilities. 


SOLID-FUELED ROCKET 
MOTOR CASES—BFG offers 
complete capability in filament- 
wound rocket motor cases — 
also insulators. These glass 
fiber reinforced plastic struc- 
tures are lightweight, strong, 
and are produced in many 
sizes and configurations. 


BEGoodrich aviation products 
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the energy of a proton flux so as to 
make the interior environment of the 
vehicle sufficiently tolerable for a 
specified time and mission. A need 
exists for a comprehensive set of curves 
(similar to those presented here) showing 
the transmitted euergy of protons for 
a wide range of energies after having 
traversed increasingly greater thick- 
nesses of several elements and struc- 
tural materials. 

Figs. 4, 5, and 6, compiled from the 
data,* show graphically the remaining 
energy of protons after penetrating 
increasing thicknesses of lead, alumi- 
num, air, and water. Incident energies 
have been taken as 30, 50, 100, 300, 
500, 1,000, 3,000, 5,000, and 10,000 mev. 
Since the major portion of the range 
of the proton in each case is traversed 
by the time the proton energy has been 
reduced to about 0.19, where Ep is 
the incident energy, the graphs have 
been plotted to permit greater accuracy 
in determining material thicknesses. 
For example, a proton with energy 
Ey = 500 mev will have its material 
thickness shown on Fig. 5b up to the 
point where its energy is reduced to 
300 mev. Then, for greater accuracy 
between 300 and 100 mev, Fig. 5c should 
be referred to, and so forth. 

It is noted from the graphs that the 
higher Z material (higher atomic num- 
ber) makes the poorer shield for protons. 
For the same energy attenuation of a 
proton, the material thickness of water 
is about only half or less than that of 
lead. 

Interaction of Protons With ‘‘ Average 
Man’’—Shielding a space vehicle for 
human occupancy depends upon many 
factors such as orbit and length of time 
in that orbit. Other strong influencing 
factors are the weight of the vehicle 
and the allowable radiation exposure 
of the human occupants. The allow- 
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able dosage, after all, is of paramount 
importance and represents the principal 
reason for requiring a shield against the 
anticipated radiation. The physiolog- 
ical processes involved in radiation 
damage to living tissue are very complex 
and can be assayed only in terms of 
average effects. Certainly the micro- 
scopic biological effects of ionizing 
radiation depend upon the type of 
radiation and its energy and, very 
signfiicantly, where and by what process 
the energy is transferred to the living 
organism. Because of averaging 
effects over the entire body for ex- 
tended periods of time—and also for 
convenience—it is assumed in the follow- 
ing that the ‘average radiation path 
length” for a high-energy proton or 
heavy particle through an ‘“‘average 
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Fig. 5. Degradation of proton energy (E) = 
300, 500, 1,000 mev). 


man”’ is 50 cm, and that the man has the 
stopping power of water. These as- 
sumptions are not unreasonable in the 
light of present knowledge related to 
effects of proton radiation on any given 
individual. 

Using these assumptions and Fig, 7, 
which gives the energy loss of protons 
in water, Fig. 8 was constructed. This 
graph peaks very sharply at 300 mey 
and falls off almost as sharply after that. 
This peak is readily explained. A 
300-mev proton will completely lose 
its kinetic energy in passing through 
50 cm of water. This means that all 
protons of energy equal to or less than 
300 mev will be completely absorbed 
by the “average man.’”’ For protons 
of energies higher than 300 mev, the 
energy absorbed by the ‘‘average man” 
drops off sharply to about 100 mev since 
these protons now pass completely 
through the body. Also, —dE/dx has 
almost a constant value for energies 
above 400 mev and that is about 1/5 
the value at 60 mev, below which it 
begins to increase sharply. Conse- 
quently, if an incident proton has suf- 
ficient energy to pass through the ‘‘aver- 
age man” and emerge on the other side 
with 60 or more mev of energy, the 
amount of energy absorbed by the 
average man will be less than if the 
proton was able to penetrate, say, 
3/4 of the body. 

Applications to a 500-Km Orbit— 
Figs. 4 through 6 show graphically the 
energy possessed by a proton as it 
traverses increasing thicknesses of alu- 
minum, lead, air, and water. Thus, 
a 100-mev proton will have lost half its 
energy after having passed through 
5.2 g/cm? of water, 6.8 of aluminum, 
or 11.3 of lead. This example also 
illustrates how much more effective 
on a weight basis low atomic number 
materials are for proton shielding since 
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Table 4. 


Range-Energy Relations for Protons in Aluminum and Lead 


Thickness for 100% 


Proton Range—————~ 
Energy Aluminum Lead Aluminum Lead 
(mev) (g/cm?) (g/cm?) (cm) (cm) 

52 3.1 5.2 1.15* 0.457* 
100 10.5 16.3 3.89 1.44 
300 66 100 24.4 8.80 
500 150 225 55.5 19.8 

1,000 430 620 159 54.6 


could not practically hope to shield. 


* This thickness will stop all the protons expected to be encountered in the 500-km 
orbit, except for the ‘less than one’’ proton/cm?-sec whose energy is known only as in 
excess of 52 mev—and which could range into the several bev region against which one 


ary cosmic radiation are not expected 
to present a radiation hazard. 

Reactions of Photons With Matter— 
There exist primarily three distinct 
methods by which a quantum of elec- 
tromagnetic energy (photon) may react 
with matter and so consequently trans- 
fer all or part of its energy to the target 
particle. These are photoelectric ab- 
sorption, Compton scattering, and pair 
production. 

The importance of each of these means 
of attenuation varies widely, depending 
upon the energy of the quantum and the 
material with which it reacts. At low 
energies (below 0.1 mev for low Z 
materials such as aluminum and rising 
to about 1 mev for lead) the predom- 
inant reaction is the complete absorption 
of the quantum by an electron, de- 
scribed as the photoelectric effect. 
At somewhat higher energies the Comp- 
ton effect predominates, while some- 
where between 2 and 5 mev the process 
of pair production begins to play the 
major role and becomes essentially 
the only way by which electromagnetic 
quanta will vanish as such. 
Electromagnetic Energy Spectrum at the 
500-Km Level—The short wavelength 
of the solar spectrum lies at about 7A. 
This integrated energy has been given 
as the order of 0.1 microwatt/cm? or 
less. 

The significant electromagnetic 
energy spectrum at 100 km is given in 
Table 5. The average photon energy 
is on the order of one kv. 

Electromagnetic Radiation Hazard— 
Since the earth receives only soft X-rays 
from the sun and the interplanetary 
value of gamma rays is negligible, the 
hazard from electromagnetic radiation 
exists only if tissue is exposed directly 
to these soft X-rays. However, almost 
any thin sheet of material between an 
astronaut and the sun will suffice as 
an effective X-radiation shield. 

Using the lower wavelength limit of 
7A yields 1.77 kev as the energy of the 
most energetic photon. The mass at- 
tenuation coefficient for such X-rays 
is greater than 1,000 cm?/g in aluminum 
(actually it is closer to 10,000 cm?/g). 
Hence, the ratio of transmitted to inci- 
dent energy through 0.01 cm of alu- 
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minum, using the much more conserva- 
tive value of 1,000 cm?/g, is of the order 
of 107", 


Solar Flares 


The phenomenon described as a solar 
flare occurs so frequently that several 
thousand, perhaps, have been observed 
within the last decade. However, they 
vary so greatly in size, intensity, and 
in their terrestrial effects that a classi- 
fication has been set up which divides 
them into three main groups. Each 
type seems to be associated with a 
discrete range of velocities. Type 3 
flares are those which contain particles 
having the highest velocities—0.1 to 
0.3 the velocity of light. These limits 
correspond to proton energies in the 
range of 5-50 mev. The duration of 
these bursts is on the order of 5-10 sec 
although a series of bursts may increase 
the time to 1 or 2 min.® The solar 
flare which took place on February 23, 
1956 (the fifth flare observed since 1942 
which was accompanied by large in- 
crease in cosmic-ray intensity) is known 
to have contained primary particles 
in excess of 15 bev. Such a flare has 
been classified as type 3+. A number 
of these flares have been observed in 
recent years. Two of the more recent 
and closely spaced ones occurred on 
May 12, 1959," and July 10-17, 1959. 
Most of these flares have occurred 
during a period of increased solar ac- 
tivity with which the IGY was selected 
to correspond, and so their frequency 
and intensity may be expected to con- 
form to some solar cycle. 

With respect to protons, then, we can 
discount the effects of type 1 and 2 
flares since the proton energies in these 
range from a few kev to the order of 
1 mev. We may include in this cate- 
gory most flares of type 3 in which the 


Table 5. Solar Radiation Spectrum at 
100-Km Above the Earth 


Wavelength Energy 
(A) (ergs /cm?-sec) 
7-10 ~10-*10-3 
8-18 ~0.6 
10-60 ~0.1 


protons associated with them do not 
exceed in energy those which are nor- 
mally present, as indicated from the 
AFSWC data, and against which we 
must shield for the protection of the 
crew. Nevertheless, the possibility of 
occurrence of type 3+ solar flares, 
with their associated fluxes of high- 
energy protons, must be reckoned with 
when planning any space mission re- 
quiring a significant residence-time 
outside the earth’s protective atmos- 
phere. 

Perhaps the best way to avoid dam- 
aging radiation on flights of short dura- 
tion due to the exceptionally energetic 
type 3+ flares is to gather statistical 
data and predict these occurrences 
with some degree of assurance that one 
will not occur during a given period, 
However, since manned space stations 
are expected to become a reality within 
the near future, shielding will have 
to be provided for the occupants against 
these high-energy protons. It is likely 
that adequate shielding can be obtained 
in the form of a thick-walled box. 
Then, if the start of a large flare is 
observed on the surface of the sun by a 
constantly alert observer, the procedure 
for the space station operators will be 
to seek the safety of the orbiting vault 
before the flare reaches maximum 
intensity. 

A heavy, thick-walled enclosure such 
as the one mentioned above is, needless 
to say, impractical for the presently 
expected mission in a 500-km orbit. 
The probability of a large flare occurring 
in any 30-day period is quite small. 

Also associated with solar flares is an 
increase in the X-ray intensity from the 
sun. Though the energy density may 
be great to 10~? ergs/cm?-sec),® 
the photon energies are in the 1-10 kv 
region and so are easily absorbed by 
the satellite walls. 


Van Allen Radiation Belts 


The Van Allen radiation belts, though 
not of prime concern for a 500-km orbit, 
do influence the general nature of 
shielding techniques for manned space 
flight because of the high intensity of 
the radiation within these regions. 
Van Allen discovered the existence of 
two distinct regions with peak in- 
tensities at about 10,000 and 23,000 
km from the center of the earth perpen- 
dicular to the earth’s geomagnetic axis. 
Actually the influence of the inner belt 
becomes important at about 1,000 km 
above the earth’s surface, at which 
point the radiation intensity increases 
rapidly with increasing altitude. Not 
until a point is reached about 60,000 
km from the earth’s center does the 
radiation intensity drop to the constant 
interplanetary value of about 3 counts/- 
cm?*-sec. An indication of the particles 
encountered at the inner and outer 
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Table 6. Radiation Intensities in the Van Allen Belts 


Particle Energy Flux 
Outer Belt 
Electrons >20 kev 1 X 10!! elec/cm?-sec 
Electrons >200 kev 1 X 108 elec/cm?-sec 
Electrons >2.5 mev 1 X 108 elec/cm?-sec 
Protons >60 mev 1 X 10? pro/cm?-sec 
Inner Belt 
Electrons >20 kev 1.6 X 10” elec/cm?-sec 
Electrons >600 kev 6.3 X 107 elec/cm?-sec 
Protons >40 mev 2 X 104 pro/cm?-sec 


peaks is given by Van Allen’ in Table 
6. Shielding against protons in these 
regions appears to present the more 
severe problem although considerable 
complexity will be encountered in shield- 
ing against the more energetic electrons 
of the outer belt because of the brems- 
strahlung produced in their deceleration. 


Conclusions 


The problem in essence revolves 
around the ‘shielding against protons. 
Consequently, the satellite wall thick- 
ness is relatively unimportant in this 
respect unless it is made thick enough 
to absorb at least those protons of 
52-mev energy. However, since it is 
not necessary to shield the entire cap- 
sule, the better approach would be to 
provide an adequate shield in those 
quarters where the astronauts will 
spend most of their time. Whether 
this interior shield be made of lead or 
aluminum will depend largely on their 
respective fabrication properties, plus 
the very significant fact that, for the 
same shielding power, about 2.5 times as 
thick an aluminum shield is needed 
as lead, but the density of aluminum 
is only 0.238 that of lead. So, whereas 
a l-cm? slab of lead of the required 
thickness will contain 5.2 g, the same 
shielding by a 1-cm? slab of aluminum 
contains only 3.1 g. 

If sufficient shielding is provided 
to absorb 52-mev protons in those 
regions where the astronauts will spend 
most of the time while in orbit—and 
provided the radiation flux is no higher 
than is believed at present, then there 
is no reason why a 4-week flight is not 
possible from a radiation hazard view- 
point. Effectively then, the only radia- 
tion the inhabitants of the capsule will 
be continuously subjected to will be 
hard cosmic radiation against which 
shielding becomes prohibitive, but which 
fortunately is of such low flux as to 
present no great potential danger. 

Large type 3+ solar flares, of course, 
remain the greatest unknown. During 
periods of reduced solar activity, the 
potential danger is not too great. Since 
present knowledge of their occurrence is 
based upon limited experience, it must 
be realized that some finite probability 
does exist for one to occur during any 
30-day period. Nevertheless, at this 


incipient stage of space excursions, it 
seems reasonable that the potential 
hazard from a type 3+ solar flare is 
somewhat overshadowed by other prob- 
lems related to successfully orbiting a 
manned capsule and retrieving the 
crew in a separate re-entry device. 
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phases as follows: 

Phase I: Preliminary studies to de- 
termine important parameters and their 
effects on vehicle design 

Phase II: Vehicle size parametric 
study 

Phase III: Performance characteris- 
tics of a two-ton payload V/STOL logis- 
tic transport. 


Phase | 


Preliminary Studies to Determine 
Important Parameters and Their Effect 
on Aircraft Design 


In this phase, the relative effects of 
the Breguet range parameters (lift/ 
drag ratio, cruise SFC and weight 
empty) on a typical VTOL design mis- 
sion radius are evaluated and indicate 
the role each has in the configuration 
design. The results of this study are 
presented in Fig. 3 as percentage 
change in mission radius with percent- 
age change in design parameter. The 
data is based on a reference VTOL mis- 
sion radius of 375 nautical miles for a 
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55,000 lb VTO aircraft. Mission calcu- 
lations were based on MIL-C-5011A 
cargo transport rules with an added 15 
min. of hover at the radius of action. 
The configuration considered had an as- 
pect ratio four wing with a take-off 
wing loading of 55 Ib/ft.2 Changes in 
the drag level, power plant SFC, and 
vehicle weight empty were then applied 
to yield the results shown in the figure. 
Each parameter was varied with the 
other parameters having the same values 
as in the reference mission. The results 
shown indicate the importance and the 
relative role each design parameter has 
on the mission performance. For 
example, a five percent change in drag, 
SFC, and weight empty, changes the 
mission radius by 2 percent, 4 percent, 
and 17.5 percent, respectively. 

Because of the fixed VTO lift-off 
weight and the large fuel quantities re- 
quired for range-free MIL-C fuel al- 
lowances (take-off, hover, landing), the 
fuel available for climb and cruise is in- 
versely proportional to the aircraft 
weight empty. Changes in weight 
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empty must, therefore, be compensated 
by an equal change in cruise fuel for the 
same mission profile. If the fuel 
initially available for cruise is small, a 
small increase in weight empty will re- 
duce the cruise fuel by a large percent- 
age. The results of this preliminary 
study indicate therefore, that compro- 
mises in design should be directed to- 
ward minimizing the aircraft weight 
empty as the prime design variable. 
For wing-installed lifting fans, the 
minimum possible wing area is indicated 
by the number and size of the fans. 
Design studies have indicated that the 
minimum area thus obtained corre- 
sponds to a wing loading of 50 to 55 
Ibs/ft? and a wing aspect ratio of ap- 
proximately four. Increasing the wing 
span to increase aspect ratio will result 
in higher weight empty. The results 
shown in Fig. 3 indicate that the per- 
formance improvement due to the re- 
sulting drag reduction is more than 
offset by this weight penalty. High as- 
pect ratio wings familiar in conventional 
short-range design therefore do not ap- 


AV. EXHAUST AV. EXHAUST. 
TEMP (°F) VEL. (FT/SEC) 
TIP DRIVEN FANS 180 550 


ROTATED THRUST (DRY) — 1200-1400 ——— 2,000 


ROTATED THRUST (HIGH 
TURBINE INLET TEMPERA- 1600-1800 2,200 
TURE DRY) 


ROTATED THRUST 


Fig. 2. VTOL engine exhaust conditions. 
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Fig. 3. Design variable effect on mission radius. 
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Fig. 4. Lift per fan variation with fan diameter. 
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Fig. 5. Fan lift/weight variation with fan 
diameter. 
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Fig. 6. Radius variation with useful load -ratio 
(6 engines, cruise with 4). 
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Fig. 7. Radius variation with useful load ratio 
(8 engines, cruise with 4). 


pear optimum for VTOL wing-installed 
lifting fan configurations. 

Based on these conclusions the Phase 
II aircraft size parametric study was 
performed with the design emphasis 
placed on minimizing the weight empty 
for practical wing planforms having as- 
pect ratios of approximately four. 


Phase Il 


Aircraft Size Parametric Study 


The effect of aircraft size and power 
plant selection on the mission radius, 
ferry range, cruise speed, and altitude 
are determined in this study phase for a 
two-ton payload transport. This study 
integrates detailed propulsion and weight 
parametric studies to yield mission per- 
formance trade-off characteristics with 
aircraft size. 

Power Plant Parametric Study—Since 
it has been shown that the aircraft 
weight empty is the prime design 
variable, selection of the best (mex. 
thrust/weight) gas generator power 


plant compatible with the tip turbine 
lifting fan was then determined from a 
detailed parametric engine study. Pro- 
pulsion system weight, thrust, SFC, fan 
diameter, and pressure ratio charac- 
teristics were evaluated for some five 
possible gas generators. A partial sum- 
mary of this study for a typical gas gen- 
erator is presented in Figs. 4 and 5 in 
terms of lift per fan and system lift to 
weight ratio per fan as a function of fan 
diameter. Variation of these param- 
eters is shown as a function of fan 
pressure ratio and number of gas gen- 
erators per fan required to achieve a 
given lift per fan. Graphical repre- 
sentation of the powerplant data in 
such a manner enables the rapid selec- 
tion of the best VTOL propulsion system 
to produce a given lift with sufficiently 
high fan pressure ratio to enable suc- 
cessful transition from vertical to con- 
ventional flight. From this study it 
was concluded that the most efficient 
lifting fan VTOL transport designs would 
be those employing versions of the GE 
J-85 gas generator. Fan diameter, fan 


= .48 
|| 4 
8 
5 -40 
+ 36 
8 76 ENG 
32 +4 
38 42 46 50 54 58 62 66 70 
VTO WT.~ |IOOOLB 
Fig. 8. Useful load ratio variation with weight. 
600 
8 
4 
rapius 400 
N. ML. 6 ENC 
DAC 
Limit 
40 50 60 70 


VTO WEIGHT — 1000 LBS. 


Fig. 9. Radius trade-off with vehicle size. 
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Fig. 10, Ferry range trade-off with vehicle size. 
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With the firing of this unique hypersonic 
missile from Cape Canaveral in February 1959, 
a new word was born. The word, AEROBALLISTIC, 
was coined by McDonnell to describe the new 
aerospace vehicle concepts demonstrated 
by this firing. Level flight was achieved at over 
Mach 5 within the atmosphere. 


The principles proved by McDonnell’s Aeroballistic 
missile are the first major breakthroughs in the 
design of multi-mission spacecraft for 
earth-to-space and space-to-earth transport of men 
and payloads. Low lift-off weight, efficient orbital 
configuration, low heat re-entry and precise 
atmospheric maneuverability are provided in a 
re-usable Aeroballistic vehicle. The Aeroballistic 
space mission terminates with a 
conventional runway landing. 


MCDONNELL 


Project Mercury and Aeroballistic Spacecraft « 
Talos Airframes and Propulsion Systems « 
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Fig. 11. Artist's conception of V/STOL transport. 


pressure ratio and number of gas gen- 
erators required for VTO were deter- 
mined for a range of VTO gross weights. 
Configuration design studies were then 
performed using this gas generator for 
both cruise and lifting propulsion. To 
ensure sufficient thrust for vertical take- 
off unsticking, hover capability, and 
transition to level flight, the design 
propulsion system selected has a thrust 
to gross weight ratio of 1.08 and a fan 
pressure ratio greater than 1.15 at 
lift-off. 

Configuration Study—This study was 
made using, as variables, the ratio of 
fuel plus cargo weight to VTO gross 
weight, the number and size of lifting 
fans, and the numer of gas generators. 
The configurations considered covered 


INBOARD ENGINE INTAKE DOORS - 


CREW ESCAPE HATCH - 


CARGO COMPARTMENT 
82" xX 90" x'30! 


LIFT FANS 


CREW ENTRANCE DOOR 


a VTO gross weight variation from 
40,000 Ib to 70,000 Ibs. The fuselages 
were designed to carry a 2-ton payload 
having a density of 2.6 lb/ft.* The 
designs considered resulted in two basic 
power-plant categories; those requir- 
ing six and those requiring eight gas 
generators to achieve VTOL. 
Performance Parametric Study—For a 
given cruise power plant system, gen- 
eralized curves of the radius variation 
with VTO gross weight were then gen- 
erated as a function of the useful load 
to VTO weight ratio. Fig. 6 and 7 sum- 
marize these results for typical VTOL 
aircraft designs having four GE J-85 
engines coupled to two aft fans for 
cruise and useful load/weight ratios 
from 0.32 to 0.44. The results shown 
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Fig. 12. General arrangement of V/STOL transport. 
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are based on the MIL-C-5011A cargo 
and transport mission rules with the ad- 
dition of a 15 min. hot day hover at 
the remote base. ‘ 

At a fixed VTO gross weight these 
curves clearly indicate the mission 
radius penalty due to an increase in air- 
craft weight empty. As concluded from 
the Phase I study, any increase in 
weight empty must be compensated 
for by an equal decrease in cruise fuel; 
hence the radius degradation. With a 
given useful load/weight ratio these 
curves also indicate the VTO gross 
weight required to achieve a desired 
mission radius. 

The estimation of the actual useful 
load to weight ratios was then deter- 
mined for six designs within the 40,000 
Ib to 70,000 lb VTO weight range. The 
resulting variation of this ratio with 
VTO weight is presented in Fig. 8. The 
configurations considered have four 
GE J-85 engines coupled to two aft fans 
for cruise and six or eight gas generators 
feeding two lift fans for VTOL. The 
configurations utilizing six lift-off gas 
generators covered a design VTO weight 
range of 40,000 lb to 60,000 Ib, and the 
eight generator configurations covered a 
weight variation from 49,000 lb to 
70,000 Ib. Since the fan diameter re- 
quired for a given VTO weight is less 
with the eight generator system than 
with the six generator system, the result- 
ing decrease in structural weight yields 
the increase in useful load ratio shown 
in Fig. 8. The radius trade-off with 
aircraft size for the design power plant 


PITCH CONTROL JETS 


POWER PLANT 
Lift fan system G, E, No, 166 
6 G.E, J85 Phase II engines 
4 used with 2 aft fans for cruise 
2 G.E, 119" dia, lift fans 


WING 
Airfoil section--NACA 66(215)--416 
Thickness ratio --Inboard 16%--Tip 12% 
Area--980 sq, ft. 
Aspect ratio--4, 18 


HORIZONTAL TAIL SURFACES 
Thickness ratio--14% 
Area--244 sq, ft. 


VERTICAL TAIL SURFACES 
Thickness ratio--12% 
Area (2 at 85) Total--170 sq, ft. 


VTOL/STOL 8 HOVER CONTROLS 
Pitch-trim fan 8 Jet unit in tail 
Yaw-jet unit in tail 
Roll-jet unit in wing 
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Fig. 13. MIL-C cargo and transport mission 


profile. 


system can now be determined from 
the generalized radius useful load curves 
of Figs. 6 and 7 and the actual useful 
load curves of Fig. 8. For each lift-off 
gas generator combination, Fig. 9 sum- 
marizes this radius trade-off with air- 
craft VTO weight. The maximum diam- 
eter indicated on these charts corre- 
sponds to fan construction limitations. 
For the six gas generator system this 
limit corresponds to approximately 
58,000 lb VTO weight with a fan area 
to wing area ratio of 0.16. With the 
eight gas generator system, this limit 
corresponds to approximately 70,000 Ib. 
The results of Fig. 9 indicate a radius of 
up to 400 nautical miles can be ob- 
tained with the six generator system for 
VTO weights up to 58,000 lb. Increased 
radii above 400 nautical miles can be 
obtained only by the eight gas genera- 
tor system. For gross weights less than 
57,000 Ib, the six generator system has 
the greater radius capability although 
the useful load/weight ratio is less than 
for the eight generator system. This 
condition results from a net increase in 
mission range-free fuel allowances re- 
quired by the additional two gas genera- 
tors. At a given VTO weight the fuel 
available for cruise with the eight gen- 
erators is therefore less than that with 
the six generators; hence the reduction 
in radius. The mission profile yielding 
these radii is essentially invariant with 
VTO gross weight. For both lift-off 
gas generator combinations, the mission 
cruise conditions correspond to 350 
knots at 25,000 to 36,000 ft altitude. 

In a similar manner the corresponding 
ferry range trade-off with aircraft size 
can be determined and summarized as 
indicated in Fig. 10. The results shown 
in this figure are based on overload fuel 
weights corresponding to STO gross 
weights equal to 1.2 VTO weight. By 
using this weight ratio of 1.2, the STO 
distances are approximately 1,500 ft. 
The mission profile is calculated using 
MIL-C-5011A cargo transport rules 
excluding the 30 min sea level loiter and 
the 15 min hot day hover fuel allow- 
ances at the remote base. The 2-ton 
cargo is not carried in this mission. 

The results of this study phase im- 
mediately enables the selection of any 
2-ton payload VTOL transport weight 
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end the propulsion system required to 
achieve any mission radius up to 700 
nautical miles. In the third and final 
study phase, a more detailed design and 
performance study is presented for one 
such aircraft having a 240 nautical mile 
mission radius and a 2,500 nautical mile 
ferry range. 


Phase Ill 


Performance Characteristics of a Typical 
Two-Ton Payload V/STOL Transport 


Configuration—Based on the Phase II 
parametric study, a 50,000 lb VTO gross 
weight aircraft was selected to repre- 
sent the performance capability of a 
typical V/STOL logistic or assault 
transport. The basis for this particular 
weight selection is satisfactory mission 
radius performance and long range self- 
delivering capability. An artist’s concep- 
tion of this lifting fan transport is shown 
in Fig. 11, and a detailed general ar- 
rangement of the configuration is shown 
in Fig. 12. Except for the lift fan sys- 
tem, the configuration is conventional 
for cargo transports with ramp type 
loading. The aircraft is powered by six 
GE J-85 gas generators feeding two GE 
119-in diameter lifting fans for VTOL. 
Four of these generators are coupled 
with two aft fans for cruise propulsion. 


TOTAL CRUISE FUEL | 
33.5% 
| 

3 L TOTAL CLIMB FUEL 99% 
28 RESERVE 15.7% | 
= 

15 MIN. HOVER 25.1% 

TAKE-OFF ALLOWANCE 158% | 

Fig. 14. Basic mission fuel allotment. 
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Fig. 16. Radius trade-off with cruise speed and 
altitude. 
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Fig. 17. Radius trade-off with hover time, 


The four cruise generators are wing 
pylon mounted and the two additional 
gas generators are installed in the for- 
ward fuselage. For VTOL, the exhaust 
gas from the cruise generators is fed 
through diverter valves to the scroll of 
each fan where it mixes with the gases 
from the generators in the fuselage. 
The wing aspect ratio is 4.2 with a VTO 
wing loading of 51 Ib/ft.2 The esti- 
mated useful load to VTO weight ratio 
is 0.375. 

To provide V/STOL and hover con- 
trol, auxiliary reaction units are pro- 
vided for control about the three axes. 


A variable thrust pitch-trim fan and | 


pitch and yaw reaction jet units are in- 
stalled in the aft end of the fuselage for 
pitch and yaw control. Roll control 
jets are located at each wing tip. The 
design details and effectiveness of these 
units are beyond the scope of this paper. 

Design Mission—The design mission 
is the MIL-C-5011A cargo and trans- 
port mission with the addition of a 15 
min hover on a 90°F day at 3,000 ft 
altitude at the remote base. The 2-ton 
payload is carried outbound only. The 
240 nautical mile radius mission profile 
is depicted in Fig. 13. Mission cruise 
conditions vary from 350 knots at 
25,000 ft outbound to 365 knots at 
36,000 ft altitude inbound. 


An indication of the fuel allotments for | 


each phase of this mission is presented 
by the bar graph of Fig. 14. From this 
graph it is seen that only 33.5 percent 
of the total fuel is available for cruise, 
9.9 percent for climbing 25.1 percent for 
hovering, and the remaining 31.5 per- 
cent for MIL-C range-free fuel allow- 
ances. It appears from this fuel break- 


down that the use of conventional fixed | 


wing aircraft take-off and landing MIL- 
C fuel allowances unduly penalizes the 
aircraft mission capability. By pet 
forming the same mission using MIL-C 


rotary wing take-off and landing fuel 


~ 


allowances, the fuel available for cruise i 
is increased to 49 percent and the radius) 


is increased from 240 to 350 nautical] 


miles. If the rotary wing mission rule 
were applied to the aircraft size trade-oli 
study previously presented in Fig. ! 


the reference 50,000 lb VTO gross weighit 


would decrease to approximately 47,000 
Ib for the same 240 nautical mile radius 
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For VTO weights of the order of 58,000 
Jb, it is also indicated that the different 
mission rules could dictate the selection 
of the number of gas generators re- 
quired for lift-off. The need for a defi- 
nite, realistic set of military mission cal- 
culation rules specifically for VTOL air- 
craft is evident from this studys 


High Speed and Climb Capability— 
The variation with gross weight and 
altitude of the high speed and rate of 
climb capability for this V/STOL trans- 
port is summarized in Fig. 15. For alti- 
tudes below 30,000 ft the high-speed 
capability of this aircraft is approxi- 
mately 450 knots. Sea level military 
power rate of climb at a 44,000 Ib gross 
weight is approximately 5,500 fpm 
and the absolute ceiling is 37,000 ft. 
The best climb speed is seen to increase 
from 260 knots at sea level to 320 knots 
at 36,000 ft. 


Cruise Characteristics—Radius varia- 
tion with cruise speed and cruise altitude 
is presented in Fig. 16 for sea level, 
15,000 ft, 20,000 ft, and 36,000 ft cruise 
altitudes. Increasing the cruise altitude 
from sea level to 36,000 ft is seen to in- 
crease the best cruise speed from 240 
knots to 350 knots. The corresponding 
radius increases from 155 to 240 nautical 
mniles. 

Hover Capability—As previously indi- 
cated by the bar graph of Fig. 14, 25 
percent of the total fuel is consumed 
during the remote base 15 min. hover 
condition. If this hover fuel is utilized 
for cruise, a substantial increase in 
mission radius can be realized. A 
radius trade-off with hover time study 
was therefore performed and the results 
are presented in Fig. 17. The radius 
is seen to increase linearly from the 240 
nautical mile design mission hover time 
to 390 nautical miles for zero time. 
The rate of increase is approximately 10 
nautical miles per minute of hover. 


This paper has presented the high- 
lights of a detailed design-performance 
parametric study of the application of 
the lifting fan-turbojet concept to the 
design of a two-ton payload V/STOL 
logistic transport. A re-analysis of the 
telative importance ‘of the basic per- 
formance variables (drag, SFC and 


weight empty) has indicated the aircraft 


Weight empty as the prime design varia- 
ble for optimizing mission performance. 
The high aspect ratio and relatively 
heavier wing, familiar in conventional 
design, is not considered optimum for 
wing-fan, short radius VTOL transports. 

The VTOL propulsion system per- 
formance data is presented in a manner 
Which immediately indicates the best 
(max. thrust/weight) design lifting fan 
system for VTOL. The procedure used 
also enables the selection of the corre- 
sponding fan pressure ratio, fan diam- 
eter, and number of gas generators 
per fan required to produce a given lift. 


Parametric performance data obtained 
from configuration design studies using 
two lifting propulsion systems is shown 
for a VTO gross weight range of 40,000 
Ib to 70,000 Ib. The designs considered 
were based on a power plant system 
(considered optimum) using versions of 
the GE J-85 gas generator for both lift- 
ing and cruise propulsion. The results 
of the study enable the selection of any 
VTO gross weight and propulsion system 
required to achieve any mission radius 
up to 700 nautical miles. 

The parametric study revealed that a 
VTOL design is more sensitive to mis- 
sion rules regarding range free fuel 
allotments than are conventional air- 
craft designs. Large range free fuel 
allotments combined with short radius 
requirements cause the aircraft weight 
empty to become far more important in 
design than aerodynamic and engine effi- 
ciency. In designing VTO aircraft, 
serious consideration should therefore 
be given to the design mission rules at 
the outset. A definite set of mission 


rules solely for VTOLt aircraft is needed 
for competitive designs and these rules 
should closely approximate the expected 
operational requirements. 2 

The lift fan-turbojet propulsion sys- 
tem achieves vertical thrust with lower 
downwash velocity and at considerably 
lower exhaust gas temperature than 
would result from direct action turbo- 
jets. The installation of this system 
yields a clean aerodynamic configura- 
tion capable of high speed cruise as well 
as high maximum speed for assault 
type operations. Perhaps the greatest 
attribute of the lifting fan system is its 
ability to augment the basic gas genera- 
tor thrust by a factor of almost three and 
still maintain a high specific system 
thrust of almost nine. Continued 
research and engineering with this sys- 
tem envisions considerable growth po- 
tential with a capability of achieving 
specific thrust ratios of up to fifteen. 
The use of the fan system is an attrac- 
tive solution for high performance VTOL 
aircraft. 


++ + 


Biaxial Stress Problems (Continued from page 19) 


inflation, expansion and repackaging can 
be accomplished repeatedly. This is an 
old story to balloonists; Alberto Santos- 
Dumont once built a man-carrying bal- 
loon which he could and did pack into a 
valise* to carry home after a flight. 

In recent years, the development of 
plastic films such as Mylar and poly- 
ethylene has brought into the field a 
competitor of fabrics, and in some ap- 
plications, notably stratosphere _bal- 
loons, films have completely displaced 
fabrics. The decisive factors in bal- 
loons were the low cost of films and 
secondly, their lightness in the sense 
that they can be made thinner than fab- 
rics; the stresses in open-necked bal- 
loons are low. But fabrics have one 
more structural advantage: better tear 
resistance in pressurized structures. 


Test Strength 


Tear strength is one of those qualita- 
tive empirical measurements which often 
seem so nebulous to a stress analyst. 
It corresponds to the brittleness-duc- 


‘tility relationship in homogeneous ma- 


terials. Tear strength depends on the 
stress concentration at the root of the 
tear, which is relieved by yielding in a 
ductile material; but homogeneous duc- 
tile material tends to become brittle as 
the octahedral normal (hydrostatic) 
stress increases in the direction of hydro- 
static tension.6~? Pressure vessels in- 
herently require biaxial tensions, so 
that the seriousness of any stress con- 
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centration is aggravated. With a ma- 
terial like fabric, however, each individ- 
ual yarn is stressed uniaxially and so 
does not lose its ductility as rapidly 
although the fabric as a whole is 
stressed biaxially. 

An objection to this argument is that 
high-tenacity fibers, like other high- 
strength materials in general, have little 
ductility anyway. (The phenomenon 
of decreasing ductility with increasing 
strength is, incidentally, exactly what 
one would expect from the hydrostatic 
tension effect.) Let us examine the mat- 
ter further by first considering the 
problem of an elliptical hole in an iso- 
tropic elastic plate subjected to biaxial 
stresses in the ratio 2:1 (Fig. 1). If 
S, = 2S,, the stresses at the ends of the 
axes of the ellipse are’ 


Sz [(2a/b) + (1/2)] oy 
Sz [(b/a) — (1/2)] oz 
Eq. (1) can be written 

S:/oz = 2/[1 + 4(a/b)] (3) 


Consider next a single ply of woven 
fabric having C threads per unit width, 
subjected to the same state of stress 
S, = 2S,. Assume S, along the warp 
and S, along the fill. Let one warp 
thread per inch be cut, then the stress 
in the fabric is increased to 


oz = [C/(C — 1)] Sz (4) 


if all the other threads share the broken 
thread’s load. But suppose only the 


Or 


0 (1) 
0 (2) 


Gy 
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Fig. 3. Clement Bayard airship (circa 1908), nose tilted up by shear. 


two adjacent threads share it. Their 
loads become o,/C = 1.5 S,/C. If two 
warp threads are cut, and the two adja- 
cent threads take the load, they must 
carry o:/C = 2.0 S,/C. On this basis, 
one can obtain a relation analogous to 
Eq. (3) by inscribing the ellipse of Fig. 
1 so that 26 is the distance between fill 
threads and 2a is the length of the slit 
between uncut warp threads. If n 
threads on either side of the slit pick 
up the load, m threads being cut, each 
effectively assisting thread takes 


= [1 + (m/2n)](S2/C) (5) 
Then 
oz/S; = (2n + m)/2n (6) 
in which 
m = (a/b) — 1 (7) 


Fig. 5. Two-place Infilatoplane in flight. 
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Fig. 4. 


so 


Similar equations can be written for 
slits normal to the Y axis. Fig. 2 
shows Eq. (8) plotted for » = 1 and 


The ZPG-3W airship (bottom), with a ZPG-2W and ZSG-4 (top) above. 


S:/oz = 2n/(2n + a/b — 1) (8) 


Note that a/b is always an integer >]. 


n = 2° higher values of » obviously 
yield lower stress concentrations and 
higher expected tear-strength/tensile 
strength ratios. Eq. (3) for elastic iso- 
tropic sheet is also shown in Fig. 2, and 
the uppermost curve is from burst test 
data on slit cylinders of a 3-ply Fortisan- 
neoprene fabric for which C is assumed 
to be 32 threads per inch so that the 
shortest slit (1/16-in.) cut one thread 
per ply only. The shape of this curve 
agrees well with the theory. 

The left legs of the theoretical curves 
take into account the stresses at the 


ends of the other axis of the ellipse; 
but the right leg is the significant one. 
The maximum S,/o, attainable for an 
elastic isotropic material is 0.667, and 
occurs when a/b = 1/2; while ina 
fabric, it can reach 1.0 and is always 
well above that for the isotropic sheet. 
If o, equals the uniaxial strength of the 
material, it is clear that the fabric will 
take much the higher stresses without 
local failure occurring. 


Plied Fabrics 


Structural fabrics are generally of | 
more than one ply, with at least one of 
the added plies at an angle with the 
straight ply. If we continue to ignore 
the effect of the coating and bonding 
material, the multi-ply fabric should 
have as good or better tear charac 
teristics than the single ply, since a slit 
across one ply, say the straight, is ef 
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fectively shorter for the diagonal or bias 
ply. The experimental curve in Fig. 2 
is for a fabric with a bias ply. Ander- 
son, et al., have found, however,? that 
tear strength may be either increased or 
decreased by varying the properties and 
thickness of the coating material. Rub- 
ber, neoprene, Hypalon and, in the case 
of metal fabrics for high-temperature 
applications such as re-entry, silicone 
compounds are the most common coat- 
ing materials. 


Shear 


Fabrics are not plied primarily to 
get additional strength; single-ply fab- 
rics are available that are stronger than 
necessary for most current applications. 
What single-ply fabrics lack is the 
ability to resist forces applied at an 
angle to its natural axes; that is, to 
the warp and fill thread directions. 
When the principal stress directions are 
known and fixed in a particular struc- 
tural design, this need not be trouble- 
some; but generally, in practice, “it 
manifests itself as low shearing stiffness. 
The addition of a bias ply corrects this. 
Fig. 3 shows an early French pressure 
airship, the Clement Bayard I, whose 
envelope was built without a bias ply. 
Notice the pitching up of the nose of the 
ship due to shear. According to Haas? 
such effects had the serious disadvan- 
tage of making the airship hard to con- 
trol. Fig. 4 shows the largest pressure 
airship ever built, the ZPG-3W, four of 
which were delivered to the Navy 
within the past year. Its envelope is 
constructed with a bias ply, and there 
is no noticeable distortion of the enve- 
lope although the stresses to which it is 
subjected are much greater than in the 
Clement Bayard I. 

In inflatable airplanes, Fig. 5, there 
are also important shear resistance prob- 
lems, as Fig. 6, showing an early in- 
flatoplane in the wind tunnel at Langley 
Research Center, indicates. Here the 
shear problem has to do with the drop 
threads which preserve the Airmat wing 
cross-section. 

The necessity for a bias ply to resist 
shear immensely complicates the prob- 
lem of estimating the strength of a 
fabric. The efficiency of bias plies in 
contributing to the strength of a cylin- 
der being burst by internal pressure was 
teported in various papers published be- 
tween 1909 and 1930 to be 100 percent," 
50 percent,!! 0 percent,!* 2 and 30 per- 
cent.5 A more recent analysis of 
cylinder burst test data at Goodyear 
Aircraft found that it depends mark- 
edly on the properties of the basic 
fibers, and that in two-ply fabrics it 
averaged 105 percent for cotton, 65 
percent for nylon, 50 percent for Forti- 
san, and 80 percent for Dacron.!4 

The. bias ply further befogs the ques- 
tion of the: stiffness of the composite 
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fabric. In order to understand the be- 
havior of multi-ply fabrics, one must 
first understand single-ply fabrics, how- 
ever; and the rest of this discussion will 
be devoted to that subject. 


Analysis of Single-Ply Fabrics in 
Shear 


To describe the behavior of cloth, 
Weissenberg™ suggested a ‘“‘trellis 
model,” which assumes the threads to 
be pinned at their intersections and that 
there is no other interaction. Chad- 
wick and Shorter" further discussed 
this model, but did not extend its use 
beyond simple pulls. All three of these 
investigators seem to have been un- 
aware of the earlier work of Haas,? who 
considered three mechanisms of defor- 
mation. These he called thread shear, 
thread straightening, and thread exten- 
sion; and the first of these may be re- 
garded as a lattice or trellis model. 

Peirce lists four deformation mech- 
anisms for cloth in the order of their 
importance as (1) angle shear and 
thread bending in the plane of the cloth 
(with a homogeneous state of stress, 
the latter cannot occur); (2) crimp re- 
distribution; (3) thread compression, 
with flattening and extension; (4) fiber 
extension. Working from considerations 
on the structure of woven cloths ex- 
pounded by Peirce in another paper,! 
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Womersley™ developed a theory which 
takes into account the first two of these 
four mechanisms and neglects the last 
two. Its application to the problem of 
the determination of the shape of and 
tensions in an originally cylindrical 
seamless piece of cloth, given a longi- 
tudinal extension and fixed-diameter~ 
ends, required numerical methods. 
More recently, Baker and Mikola- 
jewski?! have studied fabric stress-strain 
relationships assuming zero shear strains 
and Owada?? has investigated thread 
straightening in woven wire fabrics. 

Another approach is to assume the 
applicability of the equations of ortho- 
tropic elasticity; this is commonly done 
for glass cloth laminated with rigid 
resin (see, for example, Dietz**). The 
assumption is dubious because of the 
nonhomogeneity of fabric, but es- 
pecially so when a shear resistant resin 
is not present. 


Shear Resistance of a Single Ply 


If we now confine our attention to 
homogeneous states of stress and the 
problem of shear, it will be convenient 
to return to two equations first derived 
by Haas—one from statics and geometry, 
and the other from strain compatability 
requirements. They can be obtained 
from Figs. 7 and 8. 
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Fig. 6. One-place Inflatoplane in wind tunnel under load. 


tan a tan B = o;/o2 (9) 


cos a/cos B = cot 6 (10) 


Eq. (9) is quite general, being simply an 
equation of equilibrium; Eq. (10), how- 
ever, is not, though Haas does not say 
so. 

In any body in equilibrium subjected 
to plane stress, there is a line in the 
plane of the stress that remains parallel 
to its position before loading. The 
orientation of this axis depends on the 
boundary conditions. In Eq. (10), the 
angle 6 is defined as the angle between 
the normal to this axis and the nearest 
thread set, say x. Eq. (10) is valid only 
when one of the principal stress axes is 
parallel to the irrotational axis, labeled 
DG in Fig. 8. Since this excludes shear, 
the value of the results will be limited. 
Let us nevertheless see where Eqs. (9) 
and (10) lead. 


Noting that Eq. (10) also requires 
that the warp and fill be 90° apart to 
begin with, we may use the well-known 
relations 


o1, 02 = [(or + oy)/2] + 
V — oy)/2\? + try? (11) 


(12) 


tan 20 = 2 rry/(o2 — oy) 


Further, by defining the shearing strain 
Yry as the change in the angle between 
warp and fill (see Fig. 8) 


Yry = @ B (1/2) (13) 
Using the trigonometric identity 
cos a = (1 + tan? a)7!/2 (14) 


with Eq. (10) gives, after squaring both 
sides, 


1+ tan?s = (1 + tan*a)cot?@ (15) 


o,dsinB 


6, dcosG 


Fig. 7. Statics of two-thread sets. 
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Solving Eqs. (15) and (9) simultaneously 
gives tan @ and tan in terms of 
and 6—i.e., the given state of stress, 
Then a, 8, and, from Eq. (13), yz, can 
be found. Eqs. (11) and (12) yield r,,, 
oz, and oy. 

Fig. 9 shows 72,/o, as a function of 
for various o:z/o,. The right-hand 
boundary corresponds to the limit where 
o2 = 0; in a fabric this is usually as- 
sumed to be where the fabric wrinkles. 
Since both thread sets are in the o, 
direction according to the theory, a = 
B = 2/2 = yey at this limit. If one 
examines the region of moderately small 
strains, say yz, < 0.2, he finds that the 
slope of each curve in Fig. 9 is very 
nearly equal to o,/cy, or that 


= (16) 
If we define a shear modulus 
Ge=e/y (17) 
we at once get 
Gry = oz (18) 


o, being the stress in the x thread set 
based on the original configuration. 

Let us consider one more special case, 
that with one thread set lying along the 
irrotational axis A—A, the other normal 
to A-A, and subjected to normal ten- 
sions o, and oy. Now let a shearing 
stress parallel to A-A be added. The 
threads parallel to A—A cannot rotate 
if the state of stress is homogeneous, and 
if they do not rotate they can resist no 
shear since they have been assumed 
capable of resisting only tension. Equi- 
librium of the element shown in Fig. 10 
requires that 


tan yey) = Try(u) (19) 
from which, for small yz,, 


Try/Y2y = (20) 
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Fig. 8. Two-thread sets—one principal stress 
normal to irrotational axis. 
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stress 


This is the same as Eq. (16). Eq. (18) 
therefore applies also to torsion, pro- 
vided the thread set x is normal to an 
jrrotational axis before loading. 


Comparison With Experiment 


To experimentally verify this equa- 
tion an inflated single-ply cylinder with 
the warp and fill threads at an angle of 
45° with the cylinder axis was tested 
by pulling on it axially. But even the 
slightest initial pressure will, at least 
theoretically, result in rotating or ‘‘rack- 
ing’ the threads to the equilibrium 
angle for a 2:1 stress ratio (about 
54.8°). The original measurements 
should be taken before inflation, but 
this is not really feasible either. The 
presence of the neoprene coating may be 
expected to hinder thread rotation so 
that (a) this angle will not be 54.8°, and 
(b) G will be higher than the theoretical 


py 


Al \ 


_Trrotationat Axis 


Fig. 10. Two-thread sets in shear—thread 
sets parallel and normal to irrotational axis. 


value; but a coating is necessary to hold 
the inflation pressure. The actual test 
results are shown in Fig. 11. The origin 
used assumes that the threads were 
racked to 54.8° by pressure stresses be- 
fore any additional axial load was ap- 
plied. In this test, ¢z/o, = 1, and the 
relationship found between and 
Yz is in agreement with Fig. 9 to the 
extent of being reasonably linear, but G 
is larger than predicted by Eq. (16). 
Fig. 12 shows the results of torsion tests 
on a single-ply nylon-neoprene cylinder 
at various pressures. The experimental 
G’s are average values from two com- 
pletely reversed cycles of loading. The 
effect of the neoprene would be ex- 
pected to be to increase the stiffness by 
a constant amount in both kinds of test. 
Considering this, the theory is in ex- 
cellent agreement with experiment when 
the stretched diameter is taken into ac- 
count, although pressure affects G some- 
what more than was expected. 
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Fabrics having bias plies, incidentally, 
have much higher shearing stiffness, but 
that question is beyond the scope of this 
paper. 


Conclusions 


Fabrics may be expected to have 
better tear resistance than films or sheet 
material having the same stress-strain 
characteristics. The shearing stiffness 
of single-ply woven fabric, neglecting 


| shearing stress-strain diagrams. 


the effect of coatings, provided either 
one principal stress or one set of threads 
remains parallel to the irrotational axis 
of the body, can be well approximated 
by G = oz, the X axis being along the 
thread set most nearly normal to the 
irrotational axis before loading. If Eq. 
(18) is extended to an Airmat structure 
like the inflatoplane wing, assuming 
that the drop threads remain irrotational 
when pure transverse shear is applied, 
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Fig. 11. Shearing stiffness of single-ply fabric—principal stress axes fixed. 
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edge stiffness is disregarded, and o; is 
allowed to be a stress applied to the 
entire enclosed cross-sectional area, it 
becomes 

G=p (21) 
where ? is the inflation pressure. This 
result has been obtained independently 
by Leonard, et al. (references 24, 25) 
using a more rigorous approach. Shear 
deflections of Airmat beams can then be 
found from 


A= J Vdx/Ap (22) 
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Space Man vs. Space Machine (Continued from page 21) 


doubtless be many instances in which 
this is just several minutes too long 
and the whole experiment will be lost. 
Even if men went along only to main- 
tain the automatic equipment without 
making a single observation—just run- 
ning a maintenance shop—they would 
be worth their weight allocation and 
power consumption. 

Of course, to both of these arguments 
the opponents have ample rebuttals. 
For example, the machine proponent 
claims that a maintenance shop ade- 
quate for all of the automatic equipment 
likely to be found in a spacecraft big 
enough to carry man plus machinery 
would probably be larger than the 
launching rocket. Therefore, any hu- 
man maintenance man, with tools 
limited by the nature of rocketry, would 
very quickly be frustrated, unless the 
machinery were built with excellent 
reliability and numerous self-maintain- 
ing features. But in this case, of 
course, the human maintenance man 
would have nothing to do anyway. 
In fact, the really ardent machine pro- 
ponents claim that if a man is sent along 
he should be put to sleep for the whole 
trip and should never under any cir- 
cumstances be permitted to touch the 
machinery, since the probability is 
that he will hurt it rather than help it. 

There are numerous rebuttals offered 
by the proponents of manned flight. 
Quite naturally, they point to the 
record of the missile builders at Cape 
Canaveral, where automatic machinery 
has yet to prove its infallibility. Fur- 
thermore, all of this automaic machinery 
that is so blithely discussed is apparently 
no more than talk. These things may 
be possible in principle but apparently 
practice and principle are still a long 
way apart. For all of our vaunted 
capacities in the field of automation, 
there is still no complex automatic piece 
of equipment which does not at least 
have a “red light” which turns on to 
call in the human repairman. Most 
automatic equipment requires at least 
one human operator a large fraction of 
the time. Such automatic equipment 
as we do have is seldom found in a 
scientific laboratory. There is. still 
plenty left for the man to do in the 
experimental process, and, from all 
appearances, there will be for many 
years to come. 

You who have listened in to such argu- 
ments or perhaps taken part in them 
(and this is one sport where it is very 
hard to be an observer without becom- 
ing a participant) can undoubtedly add 
much more to the few points and coun- 
terpoints which I have repeated here. 
The important thing is this: Underneath 
both sides of the argument there lies a 
hidden assumption. That is, the pur- 


pose of flight to the planets is to make 
scientific observations. Of course, we 
have by assumption canceled out both 
commercial and military applications, 
so what else is left besides science? 

The answer to this question brings us 
to still another phase of the man vs. ma- 
chine argument and returns us to our 
fundamental question: ‘“‘Why are men 
interested in space at all?” 


Commercial and Military Potentials 


Both the Space Act of 1958, which 
established the National Aeronautics 
and Space Administration, and the 
pamphlet which had been issued pre- 
viously from the Office of the President, 
describing the Government’s philosophy 
on a space program, cited several ob- 
jectives for man’s interest in space. 
The Government recognized the com- 
mercial and military potentialties of 
space development as well as the neces- 
sity for scientific exploration of space 
which might lead to further commercial 
and military possibilities. Further- 
more, both these documents pointed 
out that national prestige would benefit 
from a successful space program. Be- 
yond this, the Government recognized 
the fundamental urge of mankind to 
explore, which includes but goes far 
beyond purely scientific exploration. 

Numerous space technologists who 
have made public lectures on man’s 
mission in space have referred to this 
urge for adventure and for exploration 
as one justification of their own in- 
terest in space. The statement of 
George Mallory in response to the 
question, ‘‘Why do you want to climb 
Mount Everest?”—which was, ‘“Be- 
cause it is there’—has certainly re- 
ceived ample publicity by the space 
enthusiasts. 


Pure Adventure? 


The continual reference to Mallory, 
the mountain climber, suggests that 
there is something in the lore of moun- 
tain climbing which may be comparable 
to a future lore of space travel. I would 
like to pursue this possibility by quoting 
to you some statements made by Sir 
Edmund Hillary, writing in ‘The 
Conquest of Everest.” 


“Ultimately the justification for 
climbing Everest, if any justification 
is needed, will lie in the seeking of 
their ‘Everests’ by others stimu- 
lated by this event as we were in- 
spired by others before us. From 
the response to the news of our suc- 
cess, not only in our own country and 
the Commonwealth but also in many 
other lands, it seems clear that the 
zest for adventure is still alive every- 
where. Before, during, and _ es- 
pecially after the expedition, we re- 
ceived countless gifts and messages 
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of good will and delight, in both 
prose and verse, from all over the 
world, from heads of governments 
and humble folk alike. Very many 
of these messages were sent by. ¢hil- 
dren and young people. The ascent of 
Everest seems to have stirred the 
spirit of adventure latent in every 
human breast.” 


In the final paragraph of this same 
book Hillary writes: 


. and there are many other op- 
portunities for adventure, whether 
they be sought among the hills, in the 
air, upon the sea, in the bowels of 
the earth, or upon the ocean bed; 
and there is always the moon to 
reach.” 


But perhaps the statement by Hillary 
which comes closest to describing the 
motivations which will prompt man to 
travel into space is made when he 
describes not the response of others to. 
his success or his prediction of the role 
of adventure. in the life ahead, but 
rather his own reaction to his success- 
ful conquest of earth’s highest moun- 
tain. 


“Was it worthwhile? For us who 
took part in the adventure, it was so 
beyond doubt. We have shared a 
high endeavor; we have witnessed 
scenes of beauty and grandeur; we 
have built up a lasting comradeship 
among ourselves, and we have seen 
the fruits of that comradeship ripen 
into achievement. We shall not for- 
get those moments of great living 
upon that mountain.” 


Now, although it is convenient to 
appeal to the expressions of mountain 
climbers to describe the adventurous 
spirit of man which may also lead him 
into space exploration, we may be 
indulging in too much romanticism. 
How comparable are these two activi- 
ties? 

A review of the acknowledgment list 
which Sir John Hunt includes in his 
book, ‘‘The Conquest of Everest,” 
does not show the British government 
or British governmental agencies heavily 
involved. There are a few instances, 
such as references to oxygen provided 
by the Air Ministry Gas Factory and 
help with oxygen apparatus from the 
Chemical Defence Experimental Es- 
tablishment, and so on. However, the 
participation of the government in this 
expedition was very light. Almost all 
of the support was given by private 
individuals and organizations. The 
space program is quite different. 


Justified Venture? 


One might argue that private organi- 
zations certainly have a right to finance 
adventurous mountain climbers to trek 
up Everest simply ‘‘because it is there,” 
but does a federal government have the 
right to spend tax funds to send ad- 
venturous spacemen to Mars simply 
“because it is there?” Can we accept 
the spirit of adventure as a justification 
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for manned flight to the planets? 
And can we accept the urge to explore— 
beyond the scientific urge to under- 
stand—as a justification for sending 
equipment to the planets? Or must 
we discount such motives and return to 
clearly defined scientific objectives as 
the rationale of the program? 

If we must rely simply on the scien- 
tific aspects of the space program to 
justify planetary exploration, then we 
are perhaps in an untenable position. 
For example, Robert Buckheim of The 
RAND Corporation has pointed out* 
that all of the funds expended on as- 
tronomy in the United States amount to 
approximately $10 million a year. This 
amount of money is required to launch 
small-to-medium-size artificial 
earth satellites. 

The cost of a single attempt at plane- 
tary exploration (not just going toward 
the planet, but actually getting there 
and making measurements directly of 
planetary characteristics) will certainly 
cost more than this. 

There are many planetary astrono- 
mers who point out that the capabilities 
of ground-based observatories in making 
planetary discoveries have yet to be fully 
exploited, and certainly the use of 
balloons for planetary observation is a 
very fruitful field, still in its infancy. 
Both of these techniques are far cheaper 
than space probes and can contribute 
a great deal to the science of planetology. 

Of course, nothing but on-the-spot 
observation and exploration can give 
us an adequate description of the life 
forms which probably exist on Mars. 
Intriguing though this may be, there are 
biologists who say that we shall have 
achieved the creation of life in the 
laboratory before we land on the planet 
Mars. Here, again, earth-based ex- 
periments, considerably cheaper than 
rocket-borne experiments, will offer a 
tremendous payoff in the field of biology 
even when compared with so exciting a 
prospect as the analysis of life on 
another planet. 

In general, I think it is a fair state- 
ment to say that even the scientists 
intimately concerned with the space 
program do not feel that science alone 
can justify all of the effort. and all of 
the expenses which are going into that 
program. The attitude of many is 
this: Since we are going to the planets 
anyway, we should certainly do valid 
scientific experiments in the process. 
But, why are we “going to the planets 
anyway?” 

We may have reached an impasse 
which, if it were applied to an individual 
instead of to a society, might form the 


* Buckheim, R. W., Space Exploration: 
A National Challenge, presented at 14th 
National Conference on the Administra- 
tion of Research, Ann Arbor, Mich., 
Sept. 20, 1960. 
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basis of a neurosis. Perhaps the funda- 
mental reason for going to the planets 
is that we want to. And yet it is 
difficult to accept ‘‘want to” as a justi- 
fication for the expenditure of federal 
tax money. 

It is not only in the field of space 
exploration that this dilemma presents 
itself. Our economy has for the past 
many years been operating at such a 
high level that it can produce much more 
than the necessities of our society. 
This applies not only to the individual 
consumer, but also to our whole nation 
as a unit. 

As a nation, we are in a position where 
we can afford to do some of the 
things which we want to do instead of 
having to do only those things which 
we need to do. Perhaps it is our puri- 
tanical background which makes this 
extremely difficult for us to accept. 

Occasionally it appears that we are 
stretching military necessity just a bit 
far when we use it to justify manned 
lunar bases. Is it not enough to say 
that we are going to put man on the 
moon simply because we want to? 
We may be in the midst of such a 
problem with our present manned- 
satellite plans—the Mercury program. 
It is difficult to justify this program on 
either the basis of military or scientific 
need. Therefore, we are left with the 
only other need to which we have re- 
course in the space program, namely, 
national prestige. Now, the unfortu- 
nate fact is that the Soviets will very 
likely beat us in putting a man into 
orbit. Thus, all of the publicity and 
effort we are putting into the Mercury 
program on the basis of national pres- 
tige is very likely to destroy the prestige 
rather than to build it. The Russians 
will cash in, and we will be left with our 
press releases hanging out. 


Our “Want To” Indulged 


Would it not be more healthy (and 
more honest) if we were to say that we 
are going to put a man into orbit around 
the earth, because that is what we want 
todo? Might this not also contribute 
to the national prestige? Might it not 
be worth while to show the world that 
our economy, our political system, our 
national strength, is so vital that we 
can afford the indulgence of our national 
urge toward human adventure and ex- 
citement with a program for manned 
space flight? 

It is possible that this dichotomy 
between the national need and the 
national want is the basis for the ap- 
parently endless argument about the 
role of man vs. machine in space. 
Of course, it is difficult to justify the 
position of man in space if we look only 
toward the needs of a space program, 
but then it is difficult to justify the 
program itself. 


If we turn, instead, toward the wants 
of our society which find their outlet in 
the space program, then the role of man 
in space needs no justification. It is 
an objective in itself. It then becomes 
meaningless to argue that a machine 
is better equipped to satisfy the mission, 
since the mission is to put man on 
another planet. 

If man into space is an object in it- 
self, then it is so basic that it needs no 
justification on the basis of any other 
objective. By this same token, it 
becomes an objective which is difficult to 
analyze or to argue—either pro or con. 
It is possible to try to understand such 
an objective in terms of human be- 
havior. It is possible to make a direct 
analogy between the urge to explore 
the planets and the drives which have 
motivated all of our explorers, starting 
with the first man that stood erect to 
become a man and including Marco 
Polo and Adm. Richard E. Byrd. 
It is possible to understand the urge 
to adventure which will carry men to 
the planets as an analogy to the urge 
which carried men to the top of Everest 
and on a traverse across the Antarctic 
Continent via the South Pole. (It is 
intriguing to note that these two ac- 
complishments were actually the work 
of the same man, Sir Edmund Hillary, 
who ascribes to himself the occupation 
of ‘‘beekeeper,”” which proves that the 
sun never sets on British whimsy.) 

Although we may have some success 
in relating the exploratory urge of space 
travel to the exploratory urge of more 
earthly adventures, we cannot prove 
with logic or justify on the basis of need 
the role of man’s adventurous spirit 
in the space program. Either we accept 
it as a matter of faith or we reject it 
on similar grounds. 


No Rationalization 


From this point of view, the role of 
man in space as a part of the national 
space program rests on the beliefs and 
the desires of the individual citizens of 
the country. In that case, each indi- 
vidual must search his own soul to 
find the answer. There is no ration- 
alization which will lead him to it. 
We who are privileged to be connected 
with the space program certainly have 
the primary obligation for such soul- 
searching, and undoubtedly all of us 
have indulged in it. 

Let me then recall to your mind a 
picture which I am sure has been there 
many times before. One day—a day 
that is probably farther away than 
many of us would like to believe, but a 
day that may be closer at hand than 
many of our conservative cohorts are 
ready to admit—one day, a mammoth 
rocket with a man in its nose will 
pursue its thundering logarithmic his- 
tory beyond the speed of escape; 
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and, until the successful burnout of the 
last stage, the whole world will hold its 
breath. On that day, the first human 
being to head for Mars will begin the 
long, silent month of weightless coast- 
ing, and, regardless of the uniform he 
wears, his payload will carry neither 
gun nor bomb. In the months that 
follow, he will describe to his earth- 
bound neighbors the panorama of the 
universe; and it will not matter whether 
or not he reports a single scientific ob- 
servation. One day a capsule carrying 
a man will streak down through the 
atmosphere of Mars, and the incan- 
descent trail which follows that capsule 
across the Martian sky will illuminate 
the hopes and dreams of two billion 
human beings. There will come a day 
when the first man stands erect on the 
red sands of Mars—and you and I can 
help him do it. 

Perhaps we will help him because of 
the scientific enlightenment which will 
result from his exploration. Perhaps 
we will be trying to strengthen the 
national military posture or add to the 
national prestige. But is it not possible 
that we will help simply because we 
want a man to stand on Mars? 


The artist’s rendering of the 
space machine is based on the 
JPL-NASA-Hughes Aircraft Co. 
moon soft-landing “‘Surveyor.”’ 
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Boundary Layer and Geometry 


(Continued from page 23) 


sons, to have the transonic section located down- 
stream of a high Mach number, flexible, supersonic 
nozzle. With this arrangement, however, a thick 
boundary layer develops in the long, supersonic 
nozzle section and enters the transonic section. 

The analysis here presented is based on experi- 
mental data published by ARO, Inc., Arnold Engi- 
neering Development Center, USAF, Tullahoma, 
Tenn. 


Boundary-Layer Effects on Perforated-Test Section 
Performance 


Consider a rectangular transonic working section 
consisting of four perforated walls, the slope of two 
of which can be adjusted, surrounded by a plenum 
chamber in which pressure can be independently 
controlled. With this arrangement, any given 
transonic Mach number can be maintained for dif- 
ferent angular settings of wall divergence (assuming 
that adequate pressure ratio is available). 

The starting point of our investigation is the 
realization of the fact that the geometrical con- 
vergence or otherwise of the perforated walls strictly 
determines, at any Mach number, the changes in 
the boundary-layer flow on the walls, following entry 
to the perforated section. In effect, the boundary- 
layer flow, the streamwise cross-sectional area vari- 
ation, and the flow from or into the plenum chamber 
compensate each other to maintain an essentially 
constant Mach number along the perforated section. 

This interdependence is shown by experimental 
results plotted in Fig. 2, in which the ratio Am/m of 
flow removed from the plenum chamber to total 
flow entering the perforated test section is plotted 
in terms of Mach number for constant wall settings. 
It will be noticed that a unique value of Am/m is 
required* to produce a given Mach number with a 
given wall inclination. 

Thus small changes in the perforated wall inclina- 
tion are synonymous with changes in the boundary- 
layer flow over the perforated walls. This is signifi- 
cant from the point of view of model interference 
properties of the walls, since their characteristics 
are highly dependent on the flow in the immediate 
vicinity of the wall surface. For example, the 
thicker the boundary layer, the more a wall would 


* The Am/m values plotted in Fig. 2 are the minimum ones 
required at any M and 6, and are independent of the tunnel 
pressure ratio, provided the latter is sufficiently large so as 
not to affect the pressure distribution in the perforated section 
(as would be, for example, the case in blowdown tunnels 
operating at high stagnation pressures). At smaller pressure 
ratios, larger Am/m values are required at any Mach number, 
but only ‘the: flow at the downstream end of the perforated 
section is affected without influencing the boundary-layer 
flow over most of the pérforated walls. 
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behave as a free boundary, and the less linear would 
be its flow characteristics. 

The overall effects of boundary layer on the per- 
formance of a perforated working section are indi- 
cated in Fig. 3 (reference 1, Figs. 2, 3, 4), in which 
pressure distributions as measured on the previously 
mentioned, 2-percent blockage, cone-cylinder model 
are shown for three angular settings of the perforated 
walls, at M = 1.2. The boundary-layer growth 
was estimated (Table 1) and the values are indicated 
in Fig. 3, together with the initial boundary-layer 
displacement thickness (Fig. 11). A pressure dis- 
tribution closest to the theoretical was obtained 
with 6, = —30 min—i.e., when the boundary-layer 
growth was almost zero and the ratio é+t4y/¢ = 
0.7 approximately. 

In Fig. 4, the effects of wall porosity and boundary- 
layer thickness are compared using the same model 
(reference 6, Figs. 5and 6). Increasing the porosity 
has the same effect as increasing the boundary- 
layer thickness (i.e., increasing wall divergence), as 
would be expected. 


Determining Perforated-Wall Characteristics 
General Method 


In determining cross-flow characteristics of per- 
forated walls, it is usual to vary the wall inclination 
6. (positive when divergent) while maintaining con- 
stant test-section Mach number by adjusting mass 
flow from or into the plenum chamber. The pressure 
coefficient C, = (ps — Poltenum)/Qo iS measured di- 
rectly and corresponding flow deviation 6 is com- 
puted (for small values of #) from the measured 
plenum mass flow ratio Am/m as 


6 = pv/pu = (A,/A»)(Am/m) (1) 
p= t= 1/16 IN. 22.5% 
m 
ence TAPER STRIPS, SONIC NOZZLE 
Bw 
-30' 


4 4 


08 10 it 12M 


Fig. 2. Effects of minimum suction and wall setting on Mach number 
(references 4, 5). 
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Fig. 3. Effects of boundary-layer thickness on performance of per- 
forated walls (reference 1, Figs. 2, 3, 4). 
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Fig. 4. Comparison of the effects of porosity and boundary-layer 
thickness on performance of perforated walls (M = 1.2, ¢6 = t = 
1/16 in., taper strips, sonic nozzle, reference 6, Figs. 5, 6). 
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Fig. 7. Effect of Mach number on cross-flow characteristics. 
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Fig. 8. Comparison of perforated wall characteristics determined un- 
der constant and variable boundary-layer conditions. 


The foregoing expression applies, for small values of ~ 


6, at M < 1 when a sonic nozzle is used and at / > 1 
when a supersonic nozzle is used upstream of the 
perforated section. At M > 1 with a sonic nozzle, 
part of the flow is removed for expansion from J = 
1 to > 1, so that approximately 


= A,/Ay\[1 + (4/A*)][(Am/m) — 1]} (2) 
where 


A/A* = area ratio corresponding to test-section 
Mach number 


Controlled Boundary-Layer Condition Test Results 


The aforementioned examples show conclusively 
that the boundary layer may vitally affect the per- 
formance of perforated walls. It follows that in 
order to determine perforated-wall characteristics 
from the point of view of model testing, it is neces- 
sary either to control independently the boundary 
layer on the walls and suction quantity Am (i.e., 
flow deflection 6), or to measure both simultaneously 
so that wall characteristics can be obtained by cross- 
plotting. 

These requirements, in general, have not been 
satisfied in experiments, the results of which have 
been published, except in those described in reference 
3, where a small wall sample was tested in a per- 
forated transonic section. The boundary layer at 
the upstream end of the sample was measured and 
independently controlled by varying the divergence 
of the perforated walls. The sample under test was 
rotated through small angles while the main walls 
remained in fixed position and the flow through the 
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sample was measured. The sample wall was only 
3.5-in. long (in the direction of flow), so that the 
changes in the boundary-layer thickness were small 
over its length. 

The results of these tests are shown in Fig. 5 
(reference 3, Fig. 22). The wall geometry is indi- 
cated at the top of the figure: 1/8-in. diameter 
holes in a 1/8-in. thick plate are inclined at 60) 
degrees to the vertical in order to facilitate outflow, 
the porosity (based on the hole radius) being 6 per- 
cent. The values of the boundary-layer displace- 
ment thickness 6+ and of the ratio 6+/@ are indi- 
cated for each set of results. It will be noted that 
on the one hand, provided 6+/¢ < 1, the wall char- 
acteristics are essentially linear for positive and 
negative 6. On the other hand, the measurements 
at 6+/@ = 2.39 show a pronounced nonlinearity. 

Fig. 5 also indicates that there is an appreciable 
variation in the slope Cp, with 6+/¢. The slopes 
near 6 = 0 are plotted in Fig. 6 in terms of 6*/¢. 
Cop is extremely sensitive to 6+/¢ at 6+/@ < 1 and 
relatively insensitive at 6+/¢ > 1. In effect, cross- 
flow characteristics of perforated walls appear to 
be as sensitive to boundary-layer thickness as to 
the porosity and therefore, in general, specification 
of wall geometry has to be related to boundary-layer 
characteristics. 

The variation of cross-flow characteristics with 
Mach number was obtained in the same series of 
tests and is plotted for 6+/¢@ ~ 0.5 and 1.0 in Fig. 7. 
The Cog is seen to increase appreciably with / in the 
transonic speed range for the smaller 6+/¢ value. 
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Fig. 9. Effect of hole size on perforated wall characteristics in presence 
of boundary layer. 
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Test Results Without Boundary-Layer Control 


As already noted, the majority of data on per- 
forated-wall characteristics (references 2 and 7) were 
obtained from experiments in which the boundary 
layer on the test walls was not controlled. 

In the earlier tests (reference 7), perforated plates 
were installed in the movable top and bottom walls, 
the side walls being solid, whereas the later tests (ref- 
erence 2) employed only one movable perforated 
wall. In both series of tests, all Mach numbers 
were established with a sonic nozzle. 


The experimental procedure was to change values 
of C, and 6 by changing wall inclination and auxiliary 
suction while maintaining constant M. 

This procedure does not give the C,, @ curve for a 
perforated wall with a given boundary layer on it 
—i.e., under conditions relevant to absorption of the 
model flow field. As obtained, every point on a C,, 6 
curve for a constant M corresponds to a different wall 
inclination and hence different boundary layer, and 
therefore the C,, 0 curve represents the variation of 
6 with C, and with the boundary-layer thickness. 
This is the reason for which the results of such ex- 
periments indicate, in general, nonlinear wall char- 
acteristics, the slope Cg increasing with 6 increas- 
ing—i.e., as (1) the walls are converged, (2) Am in- 
creases, and (3) boundary-layer thickness decreases. 
This is illustrated in Figs. 8 and 9. ; 

In Fig. 8, characteristics of walls of similar 
geometry, as determined by the two test methods, 
are compared. The broken line, repeated from Fig. 
5, corresponds to a fixed boundary-layer displace- 
ment thickness of 0.0535 in., whereas the full line, 
indicating nonlinear characteristics, was obtained 
with the above-described test method—i.e., under 
variable boundary-layer flow conditions. 


The effect of plate geometry in relation to bound- 
ary layer is further illustrated in Fig. 9 (reference 
7, Fig. 8e). The same test method was used, the 
hole diameter varying from 0.25 to 1.0 in. with ¢ = 
t. At large values of 0, at which the boundary layer 
was appreciably thinned, all walls exhibited similar 
properties, but at small values of 0, the character- 
istic of 0.25-in. wall departed most from the others, 
indicating largest boundary-layer effects, as would 
be expected on the basis of Fig. 5. 


Boundary-Layer Growth on Perforated Walls 


Provided the amount of flow removed from the 
perforated section and the wall settings are known, 
the average boundary-layer growth on the perforated 
walls can be estimated as follows. 


Indirect Measurements 


As the most general case, consider a perforated 
section downstream of a sonic nozzle. Let the 
boundary-layer displacement thickness at the entry 
to the perforated section be 6+, and at exit 6%. Let 
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two of the perforated walls be parallel and the 
angles 6,, of the other two adjustable. Let the entry 
to the perforated section be square, h X h, and the 
length of the section, /. We can now write the flow 
balance. 

Thus 


Effective flow area at entry (at M = 1) = 
(h — 26+;)? 


Equivalent effective flow area at exit (at M = 1) = 
(h — 26*2)(h — 26+2 + 26,1) (A +/A) 


where 4 +/A corresponds to test-section Mach num- 
ber when M>1 and A+/A = 1for M <1. The 
6. is positive for divergent wall setting. Putting 
Am = total amount of flow removed, and m = flow 
at entry to perforated section, we have, in terms of 
effective flow areas at M = 1, 

Am 


m 
(h — 2644)? — — 25+2) (hk — + At/A 
(h — 28+,)? 


and, neglecting second-order quantities, 


5+, = /s{h + 20,1 — (1 — Am/m) X 
(h — 46+,)A/At} (3) 


For M < 1, this reduces to 
d+, = + 201 — (1 — Am/m)(h — 46+1)} (4) 
Should 6+; be unknown, the boundary-layer growth 
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Ad* is obtained by taking 6*; 
= 1/26.) + — (1 — Am/m)A/At} (5) 
and, for M < 1, 

Ad+ = + 1/shAm/m (6) 


Direct Measurements 


Inasmuch as only relatively few direct measure- 
ments of boundary layer were available, the fore- 
going method was used to estimate the boundary- 
layer growth in perforated test sections. The re- 
sults of these calculations are given in Table 1 and 
Fig. 10. The Aé*+ was calculated using Eqs. (5) 
and (6). Ina few cases, Eqs. (3) and (4) were used 
and 6*, — 6+, estimated, taking the values of 6*; 
from measurements at the exit of the sonic nozzle, 
Fig. 11 (reference 2, Fig. 18); a fair agreement be- 
tween the two methods was obtained (Table 1). 
Finally, the calculated average boundary-layer 
growth was compared with the development of 
boundary layer as measured at three stations in the 
perforated test section, Fig. 12. In general, good 
qualitative agreement was obtained, with good 
numerical agreement in several cases. 


Theoretically Required and Experimentally 
Determined Perforated-Wall Characteristics 


As previously indicated in Fig. 1, the theoretically 
desired perforated-wall characteristics can be deter- 
mined from analysis of the flow field of a specific 
model at a specific Mach number. On the basis of 


studies, such as shown in Fig. 1 and test results, 
Fig. 3, Goethert (reference 1) suggested that an im- 
proved performance could be obtained with a “‘dif- 
ferential-resistance” wall (Fig. 1, dash-dot line, 
C,, 0 diagram) as compared with a “‘linear’’ wall 
(Fig. 1, broken line). By providing a greater re- 
sistance to inflow from the plenum chamber than 
to outflow into it, such a differential-resistance wall 
would closely satisfy the no-interference require- 
ments throughout the flow field—i.e., in both the 
compression and expansion regions. 


The differential resistance characteristics were 
sought by means of inclined-hole wall geometry of 
the type shown in Fig. 5. The results obtained in 
Fig. 13 were, in fact, at least in this one instance, 
excellent and showed a very substantial improve- 
ment when compared with normally perforated 
wall performance (Figs. 3 and 4). 


Difficulties are encountered if an attempt is made 
to correlate the performance of normal- and inclined- 
hole walls, as determined with the 2 percent, cone- 
cylinder model at M = 1.2 with their measured 
characteristics. 


First, it is interesting to note that characteristics 
of the inclined-hole wall, with which the excellent 
results of Fig. 13 were obtained, are, under controlled 
boundary-layer conditions and with small 6*/¢ 
ratio, essentially /imear (Fig. 5). Moreover, the 
inclined-hole walls have the theoretically desirable 
property, indicated by Fig. 1, of maintaining out- 


Table 1. 
Wall Type M Am/m in. — in. 
@ = t = 1/16 in. Ref. 4, Fig. 13 0 0.8 0.009 0.027 
r = 22.5 percent, 0.9 0.0285 0.0855 
or taper strips 1.0 0.0365 0.1095 
| 0.041 0.10 
(O, Fig. 10) 1.2 0.052 0.065 
Ref. 5, Fig. 7 —30 0.9 0.0511 —0.0103 
1.0 0.0593 0.01438 
0.0672 0.015 
Ref. 5, Fig. 8 +30 0.9 0.006 0.182 
1.0 0.0175 0.216 
1.1 0.026 0.221 
1-2 0.036 0.i76 
¢@ = 1/4 in. Ref. 4, Fig. 12a 0 0.9 0.0135 0.0405 
# = 1/16 in. 1.0 0.022 0.066 
= 22.5 percent, a 0.0285 0.063 0.085 
or taper strips 1.2 0.043 0.042 0.0625 
Ref. 4, Fig. 17b —30 0.9 0.049/0.048 —0.0181 
1.0 0.0575/0.0565 0.0074 
Ref. 5, Fig. 7 Hae | 0.0625/0.0608 0 
Ref. 4, Fig. 17a +30 Be 0.0025 0.162 0.169 
1,2 0.0175 0.125 0.149 
@=t=1/4in Ref. 2, Fig. 10 0 0.9 0.013 0.039 
r = 6 percent 1.0 0.021 0.063 
a = 60 degrees De I 0.0286 0.063 
(@, Fig. 10) 1.2 0.0444 0.0465 
—30 0.8 0 0295 —0.0751 
0.9 0.045 —0.0286 
1.0 0.052 —0.0076 
0.0593 —0.0075 
12 0.0745 —0.0225 
+30 0.0045 0.158 
1,2 0.0175 0.119 
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Fig. 10. Boundary-layer growth in perforated test sections. 
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Fig. 11. Boundary-layer thickness at station X = 2.4 in. sonic nozzle 
(reference 2, Fig. 18). 
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Fig. 14. Comparison of normally and obliquely perforated walls at 
M = 1.175. 
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Fig. 12. Comparison of measured and calculated boundary-layer 
growth. 
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Fig. 13. Pressure distribution on a cone-cylinder model (reference 1, 
Fig. 20). 
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flow into the plenum chamber (positive 6) at nega- 
tive C,.* 

With regard to the slope Cog, it amounts to about 
8 (Fig. 6 for 5+/@ = 0.25) compared with the 
theoretically required one of about 2.9 (Fig. 1). 
This large discrepancy may be at least partly due 
to the effect of model flow field on the boundary layer 
on the perforated walls. Similar comparison was 
unfortunately not possible for the normally per- 
forated wall, for which no data was available under 
controlled boundary-layer conditions. 

However, comparison can be made of the proper- 
ties of the two types of wall obtained by the method 
previously described—i.e., without boundary-layer 
control. This is presented in Fig. 14. The char- 
acteristic of the inclined-hole wall is much more 
linear than that of the normally perforated wall, 
the nonlinearity of the latter being probably ac- 
centuated by the small hole size (1/16 in.) in relation 
to the boundary-layer thickness. Also, the desir- 
able region of positive @ at negative C, is prac- 
tically nonexistent. 

It would appear that, under the boundary-layer 
conditions of these particular tests, the so-called 
differential resistance wall behaved much more like 
a linear resistance one, whereas the one supposed to 
have linear characteristics had in effect nonlinear 
properties of the wrong kind due to excessive bound- 
ary layer. That, in fact, perforated-wall properties 
may become highly nonlinear in the undesired sense, 
when the boundary-layer thickness is excessive, has 
been indicated by Figs. 5 and 9. 

The foregoing discussion appears to indicate that 
good interference properties may be obtained with 
perforated walls having effectively linear cross-flow 


characteristics and allowing positive 6 at negative 


Application to Perforated Test-Section Design 


It has been shown that linear characteristics of 
perforated walls are desirable and that their proper- 
ties are sensitive to the ratio 6+/¢ for small 6+/¢ val- 
ues. For these reasons—i.e., to ensure linearity and 
an adequate degree of control of Cr, (by varying the 
slope of the perforated walls), it is suggested that a 
suitable criterion for the minimum hole diameter is 
> 26+.f 

The upper limit of the hole diameter is dictated 
by the maximum permissible amplitude of disturb- 
ances caused by the perforations at the model loca- 
tion. Experimental measurements (reference 8) 
indicate that with a thin boundary layer (6+/¢ < 
1/3) Mach-number fluctuations are reduced to a 


* This property was also observed by Maeder (references 
9-12) for normally perforated walls with negligible boundary 
layers and was attributed to free-boundary mixing which 
induces a normal velocity component. 

t Note that this is based on data for 60-degree inclined 
holes. 
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minimum (of AM = +0.002) at M = 1.2 at a dis- 
tance of about 25 hole diameters from the wall. 
This is true for both normal and inclined holes, 
provided the major diameter of the opening of the 
latter is taken. These experiments also indicate 
that the amplitude of disturbances at a given loca- 
tion increases with the boundary-layer thickness 
decreasing (for normal holes) and with the Mach 
number increasing. 

Taking a distance of at least 25 hole diameters 
from the perforated wall to the model, the maximum 
hole size is 4/50, where h = smaller wind-tunnel 
height. 

We thus have, for 60-degree inclined walls, 


6+ < $/2 < h/200 
or, for the limiting conditions, 


6+ < h/200 


We therefore conclude that a satisfactory 60-degree 
inclined hole perforated transonic test section can be 
built provided the boundary-layer displacement 
thickness does not exceed about 0.5 percent of the 
smaller test-section dimension. 

If the same criterion of ¢ > 26+ were applied to 
the normally perforated walls, we would obtain 
6+ < h/100 for the limiting boundary-layer dis- 
placement thickness. 
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Gimbals (Continued from page 25) 


Derivation of the Equations 


A complete derivation of the deformation equa- 
tions, which is based on Castigliano’s theorem, is too 


’ Jong and complicated to be gone into here. However, 


a condensed version is given, for two reasons: first, 
so that the reader may appreciate how the equa- 
tions came about, and second, to aid the reader who 
wishes to elaborate upon the derivation and its sub- 
sequent results. 

For the derivation, we consider a composite of 
Figs. 1(b) and 3(a) to treat tensile and bending loads 
simultaneously, as in Fig. 6. Thus tensile loads F, 
and bearing reaction P,, preside. 


Deformations of a Ring Using 
Castigliano’s Theorem 

Castigliano’s theorem operates on the strain energy 
associated with the bending of the ring. Thus, be- 
fore commencing with the derivation, let us write the 
expression for the bending moment in the ring. In 


Fig. 6, the bending moment M at any point in the 
ring is given by 


M = M, + HR (1 — cos a) + VRsin a (1) 
for regions 1, 2, and 3 (0 <a < 2/2). 
M = M, + HR (1 — cos a) + VRsin a — 

PR(1 — sina) — FReosa (2) 


for regions 4, 5, and 6 (7/2 < a < =) where R is the 
radius of the ring, and a is the defining angle for the 
region in question. 

Now, by Castigliano’s theorem (see Fig. 7) 


6; 0U/dF; 0; = 


where 6;, 6,, and F;, occur in direction i and are as 
shown in Figs. 6 and 7. U = strain energy in a 
beam = (M?/2EI)ds, where E = modulus of 
elasticity and J = moment of inertia of cross section 
of ring containing the bearing. 

Since ds = Rda, 


ra) 
n= 


This is the deformation equation for a closed loop. 
Placing the partial derivative operator inside the 
integral sign, we obtain 
J, oF, 2EI 
o(M*) oM R 
0 


OF, 2EI (3) 
_ f 
OM R 


Boundary Conditions 


Since the ring is symmetrical in loading as well as 
in geometry about the vertical axis, the integration 
may be conducted on one half of the ring. Note that 
for symmetrical closure 6) = 0 and 6z, = 0. 

All deformations resulting from this computation 
are given with respect to the fixed point A. Note 
that the total horizontal deflection 5y is found by 
doubling the deflection of load F with respect to X. 

We now define the coefficient of nonuniformity K* 
as the ratio of J to Im, the moment of inertia of the 
main portion of the ring. Thus J = /,,/K and 


oM 
4 2M, J, 


= km? Hy 0 (6) 
To find My and Ho, we proceed as follows. We 
take the corresponding derivatives of the moment 
[Eqs. (1) and (2)] and substitute into Eqs. (5) and 
(6). Then we perform the indicated integrations 
and obtain the simultaneous equations in Mp and A. 
If we then substitute M, and H) back into Eq. 
(3) for each subscript, we obtain: 


Horizontal Deflection 
oU OM 


Vertical Deflection 


* K is a variable having several discrete values K,. 


ORIGINAL UNDISTORTED 
RING 


RING DISTORTED 
BY LOAD 


Fig. 2. Deformation of a uniform ae ring in tension. [See Fig. 
1a). 
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Fig. 5. Circular ring with nonuniform cross section. [Loads are as shown 
in Figs. 1(b) and 3(a).] 


Fig. 3. Ring in bending in its plane 
By 
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BEARING 
fe) 
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L | 
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Fig. 4. Deformation of a uniform circular ring in bending. 
(1) A is the midpoint between top and bottom of distorted gimbal 5Ry, 
> BU yy; € = 1/2 (5Ryy — 5Uyy). (2) See Fig. 3(a).] 


Pp 


r—R(i-sina) 


2B 
Fig. 6. Analysis of the nonuniform circular ring. [Note: (1) All 
deflections taken with respect to the fixed point. (2) Loads are shown 
in Fig. 5.] 
op Fig. 7. 
y Analysis of deformations of a 


Fig. 8. Noncircular ring with nonuniform cross section. ring using Castigliano's theorem. 
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Table 1. Geometric Influence Coefficients 


da A B Cc 
10° 0.16 0.12 0.02 
15° 0.22 0.17 0.04 
20° 0.25 0.22 0.08 


Radial Deflection (6,y) 


oU R 


ow = Sp = Ey, 


KM 
oP, 


And finally, $;;, the compliance coefficient, is given 
by 


Si; = (8) 


where 6;; = deflection in direction 7 due to load in 
direction j; F; = load or force in direction j; P,; = 
bearing reaction due to load (2P;) in direction 7. 

S;;, which is also called ‘gimbal compliance’’ or 
“elasticity,” is thus the ratio of deflection to the load 
producing it. Note that this is the reciprocal of 
“spring rate.” 


Deformation Equations 


The deflections and subsequent compliances were 
evaluated for several values of ¢, and K, in order to 
cover a wide variety of gimbal configurations. 

The compliance of any particular gimbal in terms 


of its geometry is given below for nonuniform, circular 
gimbals: 


Sry = —0.1366 [1 + 1.54(Ki+ K3 — 2)]X ) 
R3/EIn 
Syy = 0.1488[1 + 24(K; — 1) + 
— 1)]R/EIp 
Sry = —0.0934[1 + 24 (K, — 1) + 
B(K; — 1)|R3/EIn 
—0.0190[1 + 24(K; — 1) + 
— 1)]R?/EIn 


Values for the constants A, B, and C are given in 
Table I for three values of @,. 

Note: Region 3 in the foregoing equations in de- 
fined as that which includes the tensile load F, and 
the bearing reaction P,,. If these loads are applied 


(9) 


Table 2. Geometric Influence Coefficients 
6 1.0 1.0 1.33 1.33 1.5 
te. 0.5 0.67 9.57 0.86 0.80 
A 0.1861 0.1753 0.1273 0.1166 0.1018 
B 0.22077 0.1730 0.2016 0.1515 0.1801 
C 0.5803 0.50389 0.53842 0 0 
D_ 0.1861 0.1753 0.2413 0.2073 0.2342 
E 0.5083 0.5039 0.5443 0.432 0.4839 
F 0.2077 0.1730 0.1671 0 . 0 
G 0.1414 0.14388 0.13873 0.13865 0.13825 
H 0.4814 0.3570 0.3982 0.2879 0.3346 
J 0.4814 0.3570 0.3982 0 0 
L 0.1419 0.1262 0.09563 0.06829 0.07096 
M 0.1690 0.1347 0.1732 0.2899 0.1470 
N 0.5403 0.4886 1.1980 0 0 
Q 0.1419 0.1262 0.1792 0.1372 0.1588 
T 0.5403 0.4886 0.5183 0.43825 0.3086 
U_ 0.1690 0.1347 0.1228 0 0 


in other directions, signs and subscripts change ac- 
cordingly. 


Nonuniform Noncircular Gimbals 

Now consider a noncircular ring, (refer to Fig. 8) 
and define the geometry of the ring as follows: 
width of ring _ 2(B + p) 7 B+ p 
height of ring 2(4 +p) Atop 

A+B+2 
r= ee (i.e., 27 = average of height 
and width) 


Let B = 


¢ = (p/r)(2p/A + B + 2p) 


Using methods similar to those presented earlier, 
we find that the compliance of such a gimbal is given 


by the following equations with the constant values 
found in Table 2. 


Sw = 
r?/EIm 

Sy = | ag 
r3/EIm | 

Se = | 
r*/EI,, 

Conclusion 


The compliances of nonuniform circular and non- 
circular gimbals are given in terms of gimbal geome- 
try. Derived on the basis of strain energy and 
Castigliano’s theorem, the results are defined in 
Eqs. (9) and (10) and geometric influence coefficients 
are tabulated in Tables 1 and 2. 


++ + 


Stainless Steel (Continued from page 27) 


Test Results 


The results of the tests are given in the form of 
stress-strain curves (Fig. 5). Of the three curves 
for each hardness, one describes the stress-strain 
curve of the strip material, one the formed column 
and the last describes the formed material in the 
rounds of the formed specimen, deduced from the 
two previous curves. It can be seen from these data 
that forming tends to increase the values of the tan- 
gent modulus, and is more effective for softer tempers 
and for longitudinal materials. 

In order to appraise properly the effect of the 
harder round material on overall design, a method 
was developed for considering this effect, plus the 
effect of spot welding, on the capacity of stainless 
steel columns. This development appears in the 
latter part of this paper. 


Effect of Spot Welding 
This study of the effect of spot welding on column 
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strength is a by-product of the investigation of the 
effect of forming. 

Accordingly, some welded and unwelded strips 
of stainless steel were tested in compression in the 
steel guides. The effect was as suspected. The 
overall modulus of elasticity of the welded specimens 
was less than that of the virgin material. The welds 
were smaller than standard and thus the tests re- 
ported herein are of qualitative, rather than quanti- 
tative, significance. In spite of this, the data were 
used in estimating the effect of the welding in pre- 
dicting columns strengths. Since the welding in 
these column cross-sections is mostly on the neutral 
axis, the effect of welding is significant only in the 
short column range. 


Specimen Preparation 


Welded test specimens (Fig. 6) were made by 
spot-welding flat strips of steel similar to those used 
in earlier tests (2'/,in. by 1/2 in.). That is, current 
was passed through the strips to form two rows of 
1/16 in. spots, spaced 7/16 in. apart. The spot size 


Table 1. Test Specimens 


Formed Columns Flat Strips 
No. of Thickness Thickness 
Hardness Specimens (in.) Number (in.) 
Annealed 1 0.025 5 0.025 
1/,-hard, long. 1 0.020 5 0.020 
1/,-hard, transv. 1 0.020 4 0.020 
1/,-hard, long. 1 0.025 4 0.025 
1/o-hard, transv. 1 0.025 5 0.025 
Full-hard long. 2 0.030 4 0.030 
Full-hard trans. 2 0.030 4 0.030 


Table 2. Stainless Steel for Tests on Rounds 


Thick- Tensile Tensile 
ness Yield Point Strength 
Hardness (in. ) (psi) (psi) 
Annealed 0.025 40,800 119,000 
1/,-hard 0.020 
1/,-hard 0.025 
Full-hard 0.030 147,200 196 ,000 
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and spacing was the same as that used in stability 
tests of stainless steel columns.* Any roughness due 


to welding was smoothed with emery paper to 


prevent friction against the guides during the tests. 


Table 3 describes the welded strips. 


The unwelded specimens were similar to the steel 
strips used in the earlier tests on forming. 
Test Equipment and Procedure 


The test equipment (steel guides, strain gages, and 
testing machine) and procedure was the same as 


was used in earlier tests on forming. 


Test Results and Conclusions 


The results of the test are presented in Fig. 7 in 
the form of stress-strain curves. 

These results indicate that spot welding reduces 
the average modulus of elasticity for stainless steel. 
The effect is least apparent for annealed material, 
becoming more pronounced as the hardness in- 
creases. 


Calculation of Column Curves Considering Effect 
of Forming and Welding 


Methods of computing column loads are, in gen- 
eral, based on the tangent modulus or reduced modu- 
lus theories. Some aircraft manufacturers, through 
their experience with aluminum columns, tend to 
accept the tangent modulus theory rather than the 
reduced modulus theory. They contend that when 
stainless steel columns exceed the load capacity 
predicted by the tangent modulus formula, it is 
probably due to work hardening of corner material 
in the columns. 

One paper has been published’ in which stainless 
steel column calculations are based on a tangent 
modulus that is found by testing a very short column 
to determine its average stress-strain curve. Duberg 
and Wilder‘ have, as a result of theoretical calcula- 
tions based on curvature of the stress-strain curve, 
reported that the computed colunin curve should lie 
between the tangent modulus and reduced modulus 
column curves and that for stainless steel it is 
roughly midway. 

The following analysis presents a method for com- 
puting the critical load for a column in which the 
physical properties of the material may vary through- 
out the cross-section, provided that the variation is 
nearly the same for any two cross-sections. A prac- 
tical example is a formed and welded stainless steel 
column; work hardening due to bending changes the 


Table 3. Welded Strip Specimens 


No. of 
Speci- Size Thickness 
Hardness mens (in. ) (in. ) 
Annealed 1 21/4 X 3/2 0.025 
1/,-hard, long. 1 21/4 X Ye 0.020 
1/,-hard, transv. 1 21/4 X 3/2 0.020 
1/.-hard, long. 1 21/4 We 0.025 
1/,-hard, transv. 1 21/4 X Ve 0.025 
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Fig. 3. Steel guide. 
physical properties of the “corner material,’ and 
welding has a local annealing effect. Since the 
amount and disposition of formed and welded mate- 
rial depends on the cross-section, columns of dif- 
ferent cross-section may have different column 
curves even though they are made of the same mate- 
rial. 

Consider the pin end column shown in Fig. 8(a). 
Assume that the column is perfectly straight and is 
compressed by a centrally applied load. If the load 
is less than the critical value the column remains 
straight and undergoes axial compression only. By 


Fig. 4. Optical strain-gage arrangement. 


gradually increasing the load, a condition of unstable 
equilibrium develops. That is, a small lateral force 
applied at the middle of the column may produce a 
lateral deflection which will not disappear when the 
lateral force is removed. We may then define the 
critical load as the axial load sufficient to keep the 
column in a slightly bent form. The critical load 
may be calculated from the differential equation of 
the deflection curve. 

Using coordinates as indicated in Fig. 8(b), the 
bending moment at any point xis Py. Ina column of 
material with a constant cross-section and constant E 
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Fig. 5. Effect of forming on stainless steel. 
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50 % 


welding on stainless steel. 


throughout the cross-section, the resisting moment 
of the material at x is El(d*y/dx*). The solution 
for the critical load in this case results in the tangent 
modulus column formula if E is taken to be the tan- 
gent modulus. 

Now consider the column in which E may vary 
throughout the cross-section. And assume that the 
fundamental theory of bending still holds; that is, 
that plane sections remain plane after bending, or the 
fiber strains caused by bending vary directly with the 
distance from the neutral axis (axis of constant 
strain during bending). The location of the neutral 
axis may be determined from the condition that the 
total increase in stress over any cross-section during 
bending must be zero. Consider a typical cross- 
section as shown in Fig. 9. The location of the 
neutral axis v-v is to be determined. If ¢ is the fiber 


(b) 


(a) 
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strain on the neutral axis, w is measured from the 
neutral axis and w’ is measured from an arbitrary ref- 
erence — vo then 


f E(e, — €) da = 0 


where f indicates integration over the cross-section 
A 


Buda =0 
A 


where £ may vary throughout the cross-section. 
Since « = u’ — h we may write 


f — h) da = 0 = [ew da Eda 
A A A 


area, Or 


Fig. 8. Column sketch. 


Fig. 9. Column cross section. 


+ 


/ 
7 
y 
50% 50 % L. 
4 | 
j 
= 
P 
| 
a 
V 
bh 
Y 
: 
Vo 


140 
7 from which h = E da/ Eda. 
A A 1875 375 
The resisting moment at a point may now be 
written as 
(d*y/dx?) f E u? da 
A 
100 O/t OF CORNERS=6 


AREA = .0538 IN? 
For small deflections we can assume E uw? dais 
A 


independent of x and y. 
The complete differential equation for the deflec- 
tion curve becomes 


COLUMN CURVE CONSIDERING 
THE INFLUENCE OF CORNER 
AND WELD MATERIAL 
> 


(d?y/dx?) f Eu? da = Py 
A 


STRESS (psi x!075) 


60 
Using the notation k? = P/ E u* da the fore- “2 
{ A pes | 
going equation may be written 
n 2 2 — 
(d*y/dx?) — = 0 
The general solution for this equation is y = A cos i a 
the kx + B sin kx in which the constants A and B must jaeant mons Teo ——/ 
ref- be evaluated to satisfy the conditions: | | 
0 
When x=0 y=0 20 40 60 80 100 120 
x= L/2 u = SLENDERNESS RATIO 
ve 
x=L y=0 
Fig. 10. Column curves for annealed material. 
The first condition is satisfied if A = 0. The third 
1on- | condition is satisfied if RL = mn where n is any 
integer (disregarding the trivial solution 6 = 0). 
We are interested in finding the smallest value of P 
which satisfies the differential equation, hence the 
smallest nonzero value of kL will be taken. Let = 
- +1; thenkL = wr. From the definition of k we have 140 
-375—> 
kL = fF =m 
la 120 
from which x - 
t=.020 
E u* da D/t OF CORNERS =6 
A 100 


AREA = .0443 IN.2 


Ze 


In general for an end-fixity constant C, 


COLUMN CURVE CONSIDERING 


80 


P= (wc/L4) f E u? da 
A 


THE INFLUENCE OF CORNER —— 
AND WELD MATERIAL 
a 
It can readily be seen that if E were constant \ 
throughout the cross section, it could be taken out- 8 60 N Ko 
= ° 
. . = 
side the integral sign and f u? da isthe moment of ” SA REDUCED MODULUS THEORY 
A 
inertia of the section. If E were the tangent modu- - New 
lus, then, the foregoing equation would yield the NXg 
tangent modulus column curve. TANGENT MODULUS THEORY RS > 
We can express E as the sum of two terms, (say) “Siz 
E=E'+ AE. In this expression E’ is, for a given 
load, constant throughout the cross-section, but 
AE may vary. This notation is convenient since E’ 
may be taken as the tangent modulus of the mate- 
Vo tial that has not been subjected to operations that L 
change its properties. Now the aforementioned pvc 


equation can be written 


Fig. 11. Column curves for 1/4-hard longitudinal material. 
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Fig. 12. Column curves for 1/2-hard longitudinal material. 


P = (x°C/L*) f (E’ + AE) u? da 
A 
P = (x°C/L’) E' u2da + AE u2 da 
A A 


P = (#C/L’) E'I + (2°C/L?) f AE u? da 
A 


Consider as a particular example the formed and 
welded stainless steel column. By experiment we 
may find an average value of for ‘“‘corner material.” 
Similarly, if we call an arbitrary portion of the cross- 
section “‘weld material’’ we may find an average 
value of E for it. Referring to these average values 
we may speak of AEw as the average change in E 
due to welding and AF, as the average change in E 
due to forming the corners. For a given load AF w 
will then have a constant value for the “weld mate- 


rial’ and will be zero elsewhere. Similarly, AZ will 


76 


have a constant value for the corner material and will 
be zero elsewhere. Then 


P = (#°C/L?)E'I + da + 
A 
(r?C/L?) f AEw u? da 
A 
P = (#C/L*)E'T + (#°C/L*)AEc Ie + 


(2?C/L?) AEw Iw 
where J, is the moment of inertia of the corner mate- 
rial about the neutral axis and Jw is the moment of 
inertia of the weld material about the neutral axis. 
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The formula is now composed of three terms: the 
first term is the load the column would carry dis- 
regarding the forming and welding operations, the 
second term is the correction for the influence of 
forming, and the third term is the correction for the 
influence of welding. Since P, E’, AEc, and AEy all 
depend on strain, it facilitates computation to write 
the foregoing equation in the following form 


L?/C = (#°/P) (E’I + AEcle + AEw!w) 


Using this form, curves (Figs. 10, 11, and 12) for 
critical column loads were computed and plotted 
for annealed, 1/4-hard, and 1/2-hard stainless steel 
material. On the same charts are column curves 
based upon the tangent modulus neglecting the 
effects of forming and of welding, and curves based 
upon the reduced modulus as recommended by 
Watter-Lincoln.® To make an accurate comparison, 
these last curves are based upon the properties of the 
material used in these experiments. 


Discussion 


Inspection of Figs. 10, 11, and 12 shows that the 
curves that are corrected to include the effects of 
welding and forming fall between the column curves 
based on the tangent modulus and those based on the 
reduced modulus, except in the L/p values below 
70 for the annealed material. The curves also indi- 
cate that as the hardness increases, the effect of 
forming becomes less noticeable. The effect of 
welding is to reduce the modulus of elasticity and 
thus, to some extent, to counteract the effect of 
forming which increases the modulus. However, 
for the examples we have presented, the spot welds 
are along the neutral axis of the sections and there- 
fore, except for the very short columns, have very 
little influence upon the column load. Thus, the 
major difference between the corrected curves and 
those based upon the tangent modulus is due to the 
effect of forming. The significance of the curves is 
more apparent when compared with the results of the 
column testing. 


Comparison of Test Data With Column Curves Based on 
the Effect of Forming and Welding 


Data obtained from column testing of a double hat 
section are plotted in Figs. 10, 11, and 12 for compari- 
son with the column curves based on the tangent 
modulus theory, reduced modulus theory, and the 
column curve that contains the influence of corner 
material and weld material. Figure 10 shows the 
comparison for 0.025 thick, annealed material, 
Fig. 11 for 0.020, 1/4-hard, longitudinal grain mate- 
rial, and Fig. 12 for 0.025, 1/2-hard, longitudinal 
grain material. The agreement between theory and 
test points is good for annealed and practically per- 
fect for the 1/4- and 1/2-hard materials. The tests 
on the rounds and welds were confined to one thick- 
ness in each hardness. 
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Conclusions 


The results of these tests demonstrate that if one 
has sufficient information on physical properties of 
material comprising the cross section of a column, 
he can predict quite accurately the load a column 
of stable cross section can sustain. The variation 
of tangent modulus throughout a column cross 
section can be approximated by experimental 
methods. From the compression tests described 
in this paper, we found on the one hand that form- 
ing tends to increase the tangent modulus, and the 
increase is greater for softer tempers and longi- 
tudinal material. On the other hand, spot welding 
tends to reduce the tangent modulus of the mate- 
rial in the neighborhood of the weld. This effect is 
least apparent for annealed material and becomes 
more pronounced for harder materials and in the 
transverse grain direction. 

The strengthening effect of forming on the column 
is more pronounced on softer materials and in the 
longitudinal grain direction. For all material the 
effect of forming becomes less pronounced for longer 
columns where the critical stress is relatively low. 

While the analysis presented may prove difficult 
to apply for design purposes and additional test data 
would be desirable to establish values for AEc, and 
AEw it would be reasonable to suggest tentatively 
the following rules for a semi-empirical approach’: 
(1) allowable column curves for 1/4-, 1/2-, and full- 
hard material should be based on the reduced modu- 
lus for grain in the longitudinal direction, and upon 
the tangent modulus for grain in the transverse 
direction, (2) for annealed material, it is suggested 
that the tangent modulus be used in preparing 
allowable column curves for slenderness ratios over 
70; for lower slenderness ratios, column curves 
based upon the reduced modulus should be used, 
although actual strength may be appreciable greater 
than predicted. These rules should be applied only 
where the column under question is clearly strength- 
ened by the forming process. Since it is difficult to 
control accurately the hardness of cold-rolled strip 
material, a criterion for use of these rules should be 
that the tangent modulus of the actual material 
being used will be greater than or equal to the tan- 
gent modulus upon which the column curves were 
based for all stresses up to the critical stress. 
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Cosmic Rays (Continued from page 29) 


ment that describes the amount of ionization neces- 
sary to produce a given biological effect. It is re- 
lated to the roentgen by a factor known as the Rela- 
tive Biological Effectiveness (RBE), which is an ex- 
perimentally derived or extrapolated quantity. The 
RBE takes into account the geometrical distribution 
of the ionization and depends in a complicated man- 
ner upon many biological properties. With the ex- 
ception of several anomalous neutron experiments, 
the RBE is always of value one or more, and one 
roentgen is always comparable to at least one rem. 
An acute whole-body radiation dose of about 400 rem 
can be expected to kill 50 percent of the exposed per- 
sonnel within one month. The minimal lethal dose 
is close to 225 rem, while illness without serious 
disability occurs for a dose of 100 rem. 

Maximum permissible doses,*~* as set forth by the 
NCRP, are those not expected to cause appreciable 
bodily injury to a person at any time during his life. 
Whole-body limits before 1958 were 0.3 rem per 
week, or 15 rem per year. For radiation workers, 
these are the revised limits: 

(1) The whole-body, penetrating radiation dose in 
any 13 consecutive weeks should not exceed 3 rem. 

(2) Thedose in any year should not exceed 5 rem, 
unless there has been no previous employment in 
radiation work, or records show incurrences of less 
than 5 rem for each year past the age 18. In these 
low-exposure cases, up to 12 rem per year, cor- 
responding to four periods of 13 weeks with a 3-rem 
dose in each period, can be absorbed until the cumu- 
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Fig. 2. Typical diurnal cosmic ray intensity variation (reference 7). 
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3. Typical 27-day cosmic ray intensity variation (reference 7). 
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Fig. 4. Annual cosmic ray intensity variation from normal for 1937 
through 1952 (reference 8). 
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Fig. 5. lonization extremes observed near the north geomagnetic pole 
during solar maximum and minimum (reference 9). 
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Fig. 6. July 1959, Forbush decrease observed at Mt. Washington, N.H. 
(reference 12). 
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Fig. 7. High-altitude cosmic ray intensity during a sudden commence- 
ment associated with the solar flare of July 14, 1959 (reference 20). 
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lated dose amounts to 5 rem for each year beyond 
age 18. 

(8) An emergency, whole-body dose of 25 rem, 
occurring only once during the lifetime of the ex- 
posed person, does not change that person’s exposure 
status. 

Average radiation intensities to produce the afore- 
mentioned doses are given as a function of exposure 
time in Fig. 1. 

The extent to which foregoing recommendations 1, 
2, and 3, regarding radiation workers, should be ap- 
plied to spacecraft occupants is somewhat question- 
able, and a few comments are in order. At a con- 
ference® sponsored by the National Aeronautics and 
Space Administration (NASA) in June 1960, it was 
pointed out that a dose of 5 rem could be safely in- 
curred at any time during one year. This statement 
tends to render recommendation 1 somewhat in- 
applicable. Current opinion, indeed, seems to favor 
limitations of 5 rem per year for routine exposures in 
space, and 25 rem as the maximum nonrecurring dose 
for a single extraterrestrial mission. However, there 
exist circumstances in which a 25-rem dose can be 
received only by chronic exposure to space radia- 
tions over a period of at least several months. The 
problem area which remains unsolved, and which 
should be resolved with space flight in mind, is 
whether or not the allowed chronic nonrepeated dose 
can be greater than 25 rem. Ifa greater dose can be 
permitted, its value should be ascertained. Close 
inspection of recommendation 2 shows that the most 
desirable astronauts, at least for routine flights, are 
those of age 30 to 40 who have had very limited 
prior radiation exposure. They can be allowed a 
12-rem yearly dose on a regular basis over a period 
of several years. Afterward, they would be free to 
incur the normal yearly dose of 5rem. Until authori- 
tative statements to the contrary are made, it would 
be assumed the 12-rem yearly dose would result from 
exposures not exceeding 3 rem each 13 consecutive 
weeks. Finally, statements are occasionally made 
that in extremely urgent cases, perhaps of a military 
nature, the acute dose could be allowed to reach the 
vicinity of 100 rem. Such a dose might be received 
from a solar flare, for example. After an exposure of 
this magnitude, no other missions could be under- 
taken for quite a long time, for even low-level in- 
tensities could not be tolerated. Besides, to attempt 
to shield a space voyager against solar flares in such 
a way as to allow only 100, and not 200, 300, or in- 
deed 400, rem (acute doses which could produce 
death) is an almost impossible task, because of the 
statistical nature of flares. 


Cosmic Radiation 


Time Variations 


Radiation time variations, either periodic or 
during states of abnormal solar activity, are presently 
being more fully observed and understood through 


satellite experiments, deep-space probes, sounding 
rockets, and balloon observations. Variances have 
been measured only in the relatively low-energy (less 
than 10'* ev) cosmic ray spectrum, and are related 
to solar conditions. There is no evidence of a change 
at high energy. The following kinds of variation 
have been noted and are illustrated, as they typically 
occur, in Figs. 2 through 7: 

(1) A solar diurnal modification’ with a resultant 
deviation of a fraction of one percent from the normal 
intensity. 

(2) A 27-day recurrence corresponding to the solar 
rotational period and having an associated deviation 
of several percent. 

(3) A deviation apparently associated with the 
solar sun-spot cycle, which is often taken to recur 
every 11 yr. This variance from normal is about 2 
percent for observatory measurements* at an at- 
mospheric depth greater than 700 gm/cm? (10,000 ft 
altitude). There is also evidence of a flux variation 
by a factor of about two at high altitudes®’ near 
the north geomagnetic pole. 

(4) Forbush decreases.'°~!* These can be as large 
as 17 percent and are often, but not always, asso- 
ciated with magnetic storms and possibly the solar 
emission of ionized corpuscular clouds. Forbush 
decreases can last as long as several days, with 
gradual recovery taking place, after some events, 
over a period of months. The magnitude of the de- 
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Fig. 8. Cosmic ray intensity as a function of geomagnetic latitude 
and altitude for near-earth positions of measurement during a period of 
low solar activity. 
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creases is the same far beyond the influence of the 
geomagnetic field'‘—"* as it is at the earth. 

(5) Ionizing radiations emitted by the sun. 
These radiations are directly related in intensity and 
frequency to solar activity. While a complete treat- 
ment here is not appropriate, the following flare 
qualities are of interest. A typical!” solar flare 
of importance 3 or 3+ produces a particle flux of the 
order 10‘ protons/cm*-sec at the top of the at- 
mosphere above a magnetic latitude of about 65°. 
Latitudes as low as 50° can be reached if there is a 
concurrent geomagnetic storm. The typical dif- 
ferential proton energy spectrum is inversely pro- 
portional to the fifth power of the energy, above 30 or 
40 mev. The flux decreases,® several hours after 
onset, by a factor inversely proportional to the 
square of the time. Approximately 40 events have 
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Fig. 9. Cosmic ray intensity above the atmosphere at an altitude of 
400 km or less as a function of geomagnetic latitude during a period of 
low solar activity. 
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occurred so far during the latest solar cycle, which 
began in 1956. In actuality, the flare spectrum 
varies quite widely from event to event, '8—** as well 
as with time”! for a single event. On July 15, 1959, 
there was an intense flare condition that lasted at 
least a full day. The maximum flux may have been 
as high as 10° protons/cm?-sec, with a corresponding 
ionization intensity as high as 10‘ roentgen,hr. 
The radiation was observed!® at an atmospheric 
depth of 200 gm/cm?, which corresponds to an alti- 
tude of 40,000 ft. After the abrupt commencement 
of another flare on May 12, 1959, the intensity de- 
creased by a factor of 10 every four hr. The flare of 
February 23, 1956, resulted in an increase in 
primaries by a factor of three for an energy range 10 
to 100 Bev and by a factor of 50 for an energy range 
1 to 10 Bev. Most flares, however, do not have an 
energy range extending above 1 Bev. The typical 
flare energy spectrum and flux suggest that a biologi- 
cal shield may weigh 50 gm/cm?(100 lb/ft?) in order 
to provide protection against most outbursts and 
allow in free space an integrated, or total, flare dose 
of the order of one rem. Relatively rare, high- 
energy flares, such as that of February 23, 1956, 
would necessitate even more shielding. It cur- 
rently appears as though flare particles are isotropic 
in direction when they reach the vicinity of the earth. 
Still another adverse characteristic is that a reliable 
flare prediction method which is valid for a period 
greater than a few days has not yet been developed. 
Some flares which even occur on the back side of the 
sun*—*8 produce particle fluxes which are observed 
on the earth. Therefore, unless considerable shield- 
ing is provided, manned space flight should take 
place only when the sun is in a relatively quiet state, 
since flare intensities are measured in tens to possibly 
thousands of rem per hour, and a harmful dose would 
be absorbed in a short time by any unshielded space 
traveler. Periods of low solar activity are ex- 
pected**—*! during 1963-67, and 1975-78. Activity 
maximums** are expected in 1971 and 1982. The 
transition time from solar quiescence to high ac- 
tivity is relatively short, compared with the more 
gradual decline from maximum to minimum ac- 
tivity. 


Steady-State Values Above the Atmosphere 


Steady-state cosmic radiation values that were 
generally accepted**—*4 for a number of years were 
based on the belief that the primary spectrum con- 
tained few particles in the energy region below one 
Bev per nucleon. This meant the ionization at 
geomagnetic latitudes greater than 60° was taken to 
be the same as that at 60°, and this indeed appeared 
to be true during 1950-52. However, in 1954 (a 
time of minimum solar activity) low-energy protons 
caused an increase®* in the ionization levels at lati- 
tudes above 60°, and even above 40°. This back- 
ground ionization during low solar activity on occa- 
sion may be as much as 1.5 to 2 times as high at polar 
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Fig. 10. Cosmic radiation intensity as a function of geomagnetic 
latitude for high altitude during a period of low solar activity. 


regions, or high altitudes, as it would be at solar 
maximum. Nevertheless, for all but low-altitude, 
low-inclination orbits the most favorable periods for 
extended space flight are still these same times of 
low solar activity, because of the great reduction in 
difficult-to-shield flare outbursts. For this reason, 
values of the ionization rate that include the effect of 
the increase during solar minimum are shown in 
Figs. 8 through 10 as functions of altitude and geo- 
magnetic latitude, both for near-earth and high-alti- 
tude positions of measurement. During high solar 
activity the background cosmic intensity at an alti- 
tude 400 km or less would be 2/3 to 1/2 the 90° value 
of Fig. 9, and nearly constant between 60° and 90° 
geomagnetic latitude. The intensity below 40° 
latitude would remain unchanged from Fig. 9. The 
intensity along the 90° latitude line in Fig. 10 would 
be reduced by 2/3 to 1/2, and all other lines of con- 
stant latitude would approach this 90° latitude line 
with increasing altitude. 

Fig. 8 shows that as the surface of the earth is ap- 
proached, there is an ionization increase, followed by 
a decrease. The increase begins at 130,000 ft, con- 
tinues down to heights of 80,000 ft (at 90° latitude) 
or 50,000 ft (at 0° latitude), and has its source in the 
shower, or cascade, formation of mesons, nucleons, 
electrons, and high-energy photons, all of which are 
created by interaction of high-energy cosmic par- 
ticles with atmospheric constituents. The decrease 
in ionization with decreasing altitude below 80,000 


to 50,000 ft comes about through atmospheric radia- 
tion absorption, while the decrease with decreasing 
magnetic latitude results from the increased shielding 
offered by the earth’s magnetic field against rel- 
atively low-energy cosmic particles. Fig. 10 shows 
that the increase in cosmic detector ionization at in- 
creasingly great distances from the earth arises from 
a combination of the decrease in the solid angle sub- 
tended by the earth and the decrease in geomagnetic 
field strength, with a corresponding decrease in cos- 
mic particle deflection. Atmospheric depths®* cor- 
responding to 130,000, 80,000, and 50,000 ft are 3, 29, 
and 120 gm/cm?, respectively, as shown in Fig. 11. 
The way the interplanetary galactic cosmic radia- 
tion varies with increasing distance from the sun de- 
pends on the interplanetary magnetic field model 
considered to be valid. A disordered field®’ of 
strength one gamma (10~> gauss), beginning at 
about 1.5 astronomical units and extending to 5.5 
astronomical units from the sun, in the presence of 
a solar wind of 500 to 1,500 km/sec velocity, 500 
particles/em* density, would lead to only a slight 
change in cosmic intensity from the earth to Mars 
(1.5 astronomical units from the sun). A solar field 
that is predominantly a dipole®* and contains mag- 
netic scattering centers leads to a significant galactic 
intensity increase with increasing distance from the 
sun. Pioneer V measurements*® indicated the 
presence of a small negative cosmic intensity gradient 
in the direction of the sun of value —(15 + 20) per- 
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Fig. 11. Atmospheric depth, gm/cm?, as a function of altitude above 
sea level (reference 36). 
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Fig. 12. Relative biological effectiveness for cosmic rays as 
a function of altitude and geomagnetic latitude during a time of low 
solar activity. 


cent per astronomical unit. This gradient was 
measured in the vicinity of the earth during the 
active portion of the solar cycle and showed, at that 
time anyway, the cosmic ray change was slight at 
different distances from the sun where flight is 
presently considered reasonable. In the following 
analysis, it will be assumed the change also is slight 
during low solar activity. 


Steady-State Biological Intensities 


An estimate of the biological whole-body radiation 
intensity as a function of altitude and geomagnetic 
latitude can be obtained from Figs. 8 through 10 only 
if the conversion can be made from the ionization 
itself, in units of roentgens, to rem, the unit that gives 
an idea of the biological effectiveness*~*? of the 
ionization. The factor of conversion, the Relative 
Biological Effectiveness (RBE), yields a measure of 
the degree of localization, or nonuniformity, of tissue 
ionization. Ionization localization along the paih 
of penetration is singularly important for heavy 
particles (atomic number 6 or greater). Although 
all atomic species through iron have regularly been 
observed** in the primary radiation, the biologically 
noteworthy heavy constituents are carbon, nitrogen, 
oxygen, the magnesium and calcium groups, and 
iron. When these medium- and high-energy par- 
ticles enter tissue, they first produce an ionization 
trail of great density. The high-energy particles, in 
general, undergo nuclear disintegration during the 
penetration process. This results in a large reduc- 
tion in specific ionization, since afterward the ioniza- 
tion is caused by several particles of reduced charge 
traveling in different directions. However, those 
primaries that have a reduced impinging energy 
have a significant probability of being stopped com- 
pletely through ionization only. This leads to ex- 
tremely large specific ionizations near the ends of 
their paths, since the rates of energy loss increase as 
the particle energies decrease, down to very low 
energies. These thin-down hits are capable of caus- 
ing cell destruction. Their effects in nonreparable 
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regions of the body, such as certain brain areas, have 
not yet been demonstrated. An RBE conversion 
from roentgen to rem is presented in Fig. 12. It isa 
whole-body RBE, and was obtained from a weighted 
analysis** of particle type and tissue ionization 
characteristics between 30° and 55° latitude at the 
top of the atmosphere, and from extrapolation else- 
where. The increase with increasing altitude and 
geomagnetic latitude is explained by noting that at 
a position requiring decreased particle penetration 
of the magnetic field, there is a slight increase in the 
relative number of heavy constituents, compared 
with hydrogen and helium. At the same time, the 
heavy component energy range extends to lower 
values. It must be emphasized, however, that little 
actual biological experimentation has been per- 
formed to test the relation between ionization track 
density and RBE for particles of large atomic num- 
ber, which produce the greater fraction of the un- 
shielded cosmic ray biological intensity. In reality, 
any whole-body RBE is a quantity averaged over all 
effects produced in all parts of the body, and, in par- 
ticular for cosmic rays, must be recognized as a 
tentative quantity at present. Where the maximum 
RBE value of Fig. 12 is 5.1, a value of 7 has also been 
suggested. The RBE during high solar activity 
can be expected to be somewhat lower than shown in 
Fig. 12, above 40° latitude. It will also have nearly 
constant value between 60° and 90° because of the 
absence of particles of energy much below one Bev 
per nucleon. 

A cosmic particle biological shield is not ordinarily 
advisable, since the radiation is isotropic and a thick- 
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Fig. 13. Approximate geomagnetic latitude as a function of geo- 
graphic latitude. 
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ness greater than 25 gm/cm? (50 lb/ft?) is required 
just for heavy particles. In fact, the biological dose 
is increased unless the shield thickness is greater 
than 15 gm/cm? for carbon, nitrogen, and oxygen; 
greater than 10 gm/cm? for the magnesium group; 
greater than 6 gm/cm? for the calcium group; and 
greater than 5 gm/cm? for iron. At least several 
hundred Ib/ft? shielded area are necessary for hydro- 
gen, helium, and the secondary radiations formed 
through high-energy nuclear interaction processes. 
Hydrogen and helium have far higher penetrating 
power than do heavier particles, and constitute 99 
percent of the primary cosmic flux, but only 15 per- 
cent of the unshielded biological dose. 


Orbital Dose Rates 


The information in Figs. 9, 10, and 12 is sufficient 
to determine the biological dose received from the 
steady-state cosmic radiation at times of low solar 
activity. However, an approximate transformation 
from geomagnetic to geographic coordinates should 
be made. This is accomplished by observing that 
the magnetic field, with the exception of several 
anomalies, is well represented above an altitude of 
several hundred kilometers by a dipole whose center 
is displaced‘**~* about 500 km along a radius vector 
intercepting the earth’s surface at close to 6° N lati- 
tude, 162° E longitude,’ and whose axis is tilted 
11.5° from the rotational axis of the earth. Dis- 
placement of the magnetic center causes a longitu- 
dinal effect in the cosmic-ray intensity. This 
amounts to about 4 percent at 0° magnetic latitude, 
and is neglected in the geomagnetic-geographic 
transformation. Approximate maximums, mini- 
mums, and averaged geomagnetic latitudes are given 
in Fig. 13 for various geographic latitudes. The 
absolute geomagnetic values averaged over the 
geographic longitudinal angle are those that are 
pertinent for discussion of a long-lifetime satellite 
whose orbital period is different from the length of a 
mean sidereal day (23 hr, 56 min in solar time). 

It has been demonstrated that long-duration flights 
should not pass through the region where high-energy, 
magnetically trapped protons are to be found, since 
effective shielding*- of the proton spectrum*® 
measured even at 1,200-km altitude is very difficult 
and requires a thickness of many gm/cm?. This 
near-earth, high-intensity penetrating component 
commences at a longitudinally dependent altitude be- 
tween 400 and 1,400 km, because of the displacement 
of the center of the magnetic field from the geo- 
graphic center of the earth, and extends*® to at least 
3,500 km. Relatively soft radiation zones also sur- 
round the earth. They are highly time-dependent in 
structure, have maximum extent 40,000 to 80,000 km 
from the earth, contain high electron fluxes, and re- 
quire shielding thicknesses on the order of 2 to 3 
gm/cm? for adequate biological protection. This 
shielding is necessary to attenuate both the brems- 
strahlung radiation arising from electron stoppage 
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Fig. 14. Average cosmic ray intensity as a function of satellite orbital 
inclination for several near-circular orbital altitudes. 
and the electrons themselves, the flux of which can 
be as great as 10''!/cm?-sec and whose energy range 
lies in the region of several tens to several hundreds 
of kev. Since locations at or below 400-km and 
above 3,500-km altitude may be safe from the 
hazards of geomagnetically trapped particles when 
only a relatively light shield is used, near-circular 
orbital dose rates are next calculated for sample alti- 
tudes of 400 km or less, 3,500, 6,000, and greater 
than 30,000 km. It is noted that the cosmic ray in- 
tensity is constant® between 50 and 400 km. 
The average unshielded, long-term orbital cosmic 
radiation intensity, I(z), is given in rem/hr by: 


1 lotr 
= T[A,(t), h(t) Jat 


where 
t = time in hours 
T = period in hours for one complete 
satellite revolution 
h(t) = satellite altitude as a function of 


time 
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Fig. 15. Allowable exposure time to cosmic rays for accumulated 
doses of 3, 5, and 25 rem as a function of orbital inclination for a near- 
circular orbit of altitude 400 km or less. 


a = inclination of the satellite’s orbital 
plane to the plane of the equator of 
the earth 

(2) = geographic latitude as a function of 


time for a constant satellite orbital 
inclination 7, and is obtained from 
the satellite equations of motion 
A,(t) = 
f[A(t)] = average geomagnetic latitude as a 
function of time as determined 
from Fig. 13, on specification of 


h(t)] = steady-state cosmic ray intensity, 
r/hr, for latitude A,(t) and altitude 
h(t), as obtained from Figs. 9 and 
10 
h(t)} = RBE, rem/r, for the steady-state 


radiation, as obtained from Fig. 12 


I(i) is plotted in Fig. 14 as a function of orbital 
inclination for near-circular orbits of altitude 400 
km or less, 3,500, 6,000, and greater than 30,000 km. 
All orbits of 400-km perigee or more will experience 
average intensities lying between those at 400 km 
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and 30,000 km. The maximum dose rate of 7.9 
mrem/hr is received at an altitude substantially 
outside the earth’s environment, taken to be greater 
than 30,000 km, and is independent of orbital in- 
clination. The safest orbit lies just above the at- 
mosphere, has an inclination of 0°, and subjects a 
man to 0.15 mrem/hr. 

Allowable orbital exposure times to the cosmic ray 
intensity alone are plotted in Figs. 15 through 17 for 
received doses of 3, 5, and 25 rem. It can be seen 
from Fig. 15 that the near-circular orbit which lies 
below the inner magnetically trapped radiation belt 
does not expose a man to a dose greater than the 5 
rem yearly allowable dosage for any inclination less 
than 48°, or greater than the allowable 13-week dose 
of 3 rem for any inclination less than 71°. The mini- 
mum allowable time in orbit at an altitude of 400 km 
or less is about 80 days for an inclination of 90°. 
Figs. 16 and 17 show that the allowed exposure times 
at high altitudes are considerably reduced. At 3,500 
km, only orbits of less than 25 degrees inclination do 
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Fig. 16. Allowable exposure time to cosmic rays for accumulated 
doses of 3, 5, and 25 rem as a function of orbital inclination for a near- 
circular orbit of altitude 3,500 km. 
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Fig. 17. Allowable exposure time to cosmic rays for accumulated 
doses of 3, 5, and 25 rem as a function of orbital inclination for a near- 
circular orbit of altitude 6,000 km. 


not expose a human being to a dose greater than 3 
rem during a 13-week period, while at 6,000 km all 
orbits have an exposure time limitation. 

Fig. 18 shows that for flights at an altitude greater 
than 30,000 km during a time of low solar activity, a 
3-rem dose is incurred in only 16 days. In just 26 
days the dose is 5 rem. Unshielded lunar trips can 
be undertaken, but excursions to other planets, where 
distances are millions of miles instead of several 
hundreds of thousands, in 130 days will subject men 
to a whole-body dose from cosmic radiation alone 
equal to 25rem. This is the recommended maximum 
allowed for a single emergency exposure. It also 
must be remembered that some additional radiation 
probably will be received in traversal of the radiation 
belts that surround the earth. For long-distance 
trips, several courses of action can be taken: 


(1) For unshielded flight, the decision must be 
made as to whether or not the first space travelers 
can take unusual risks, or at least receive more than 
25 rem. If they can, the flight time increases pro- 
portionally to the incurred dose. If they are allowed 
to receive only 25 rem, however, their space vehicle 
must be propelled over its path in about 130 days. 
It is seen that time is indeed a precious commodity 
that is not to be squandered. While voyages of 130 
days or less during minimum solar flare activity 
should be relatively safe, to keep the interplanetary 
flight time within this limit will require an extremely 
large chemical propulsion system, a high-thrust and 
high specific impulse unit such as a nuclear reactor, 
or an initial high-thrust system of relatively short 
duration backed up by a long-duration low-thrust 
device. Long-term spiral orbits created by low- 
thrust propulsion units when escaping from a plane- 
tary gravitational field are not advisable, because 
too large a radiation dose is received. Calculations 
of the minimum mass loss for interplanetary trajec- 
tories should also include the 130-day time require- 
ment if it is desired to keep within recommended 
radiation exposure limits. It appears that this time 
restriction can be met successfully by propulsive 
systems envisioned for 1975-78, at least for a Venus 
trip. The energy requirement for a short-term Mars 
voyage may still be prohibitive. Asstated previously, 
unshielded travel should not be undertaken during 
times of high solar activity (next low-activity times 
are expected in 1963-67 and 1975-78). 

(2) For shielded flight during low solar activity, 
a thickness 200 gm/cm?, roughly equivalent to the 
protection afforded by the overhead atmosphere at 
40,000 ft altitude, would attenuate the interplanetary 
cosmic radiation to a biological level nearly 1/10 the 
unshielded value given in Fig. 18. Essentially all 
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Fig. 18. Received dosage from cosmic radiation as a function of the 
time of flight at altitudes greater than 30,000 km. 
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Fig. 19. Allowable reactor radiation intensity at the position of in- 
terest as a function of exposure time for various total doses, including 
the cosmic ray background at an altitude greater than 30,000 km. 


heavy primary particles would be eliminated by ab- 
sorption or fragmentation in such a shield, and 
chronic exposure to the resultant radiation would be 
safe. 

(3) For shielded flight during high solar activity, 
the unshielded steady-state radiation would be ex- 
pected to be somewhat less than 2/3 to 1/2 the 
values of Fig. 18 because of solar modulation of 
cosmic-ray particles with energy below one Bev, and 
a 200 gm/cm? shield would further reduce the biologi- 
cal intensity to nearly 1/5 this unshielded value. If 
protection against just solar flare radiation were de- 
sired, a shield of thickness 50 gm/cm? would afford 
safety in most events. The background cosmic ray 
level, in roentgens, would be greater than the un- 
shielded level, but at the same time the proportion 
of highly ionizing heavy particles would be decreased. 
The whole-body biological intensity would be an 
estimated 1/2 to 1/4 the unshielded value. 

It is vital that more complete knowledge of the 
spatial variations in interplanetary cosmic rays, and 
their low-energy cut-off changes with changing solar 
activity, be determined. Values of the relative bio- 
logical effectiveness for cosmic rays, as well as the 
occurrence probability and physiological effect of 


86 Aerospace Engineering + April 1961 


thin-down hits, will undoubtedly become better 
known when the results of prolonged space environ- 
ment experiments are established. 


Allowable Radiations From 
Nuclear Power Sources 


It has been clearly shown that when discussing 
manned, long-term, extra-atmospheric missions, the 
cosmic radiation is a significant constituent of the 
space environment. The use of nuclear power 
sources in the same environment means that the cos- 
mic-ray level will often be comparable to the allow- 
able reactor radiations when averaged over the total 
time of flight, and it will be the sum of all radiation 
sources that determines the received dose. When 
this is taken into account, justification for the use of 
a reactor comes about through the need in the near 
future of either a high output power source to act 
over a long period, or a high-thrust, high specific 
impulse device to act over a limited time. 


Propulsive Power 


Determination of the allowed dose from a nuclear 
source of propulsive power is accomplished by first 
establishing the maximum total mission dose which 
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Fig. 20. Allowable reactor radiation intensity at the position of in- 
terest as a function of time in orbit for various orbital inclinations, where 
the total dose is 3 rem and includes the cosmic ray background at an alti- 
tude 400 km or less. 
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can be tolerated. The background cosmic ray in- 
tensity is next calculated (as in Fig. 14, for example, 
during low solar activity). The product of the cosmic 
ray intensity and the flight time gives the incurred 
cosmic ray dose, and when this dose is subtracted 
from the maximum mission dose the allowable reac- 
tor dose at the crew position is obtained. The reac- 
tor dose can, therefore, be increased as the inter- 
planetary or orbital time of travel is decreased. 


Auxiliary Power 


For auxiliary power, allowable reactor radiation 
levels also are those that, when added to the cosmic- 
ray background, yield given radiation doses at the 
crew location over a specified time. These levels 
during low solar activity are shown in Fig. 19 for the 
important interplanetary case, where the altitude 
above the surface of the earth is taken to be greater 
than 30,000 km. It is seen that essentially no in- 
crease in the allowable time of flight occurs for a reac- 
tor radiation intensity less than 10~‘ rem/hr, while 
above 10~? rem/hr the flight time decreases almost 
proportionally to the intensity increase. There is 
indeed little advantage in designing a reactor shield 
to reduce the reactor radiation level below 10~‘ 
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est as a function of the time in orbit for various orbital inclinations, where 
the total dose is 5 rem and includes the cosmic ray background at an alti- 
tude 400 km or less. 
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Fig. 22. Allowable reactor radiation intensity at the position of in- 
terest as a function of the time in orbit for various orbital inclinations, 
where the total dose is 25 rem and includes the cosmic ray background 
at an altitude 400 km or less. 


rem/hr at the crew location. All that is gained in 
such a design is additional weight with no increase 
in crew protection for the allowed flight time. 

Allowable reactor radiations, again during low 
solar activity, are shown in Figs. 20 through 22 asa 
function of time for various inclinations of near- 
circular orbits at an altitude of 400 km or less and 
for doses of 3, 5, and 25rem. The 3,500- and 6,000- 
km cases are similar, do not give any new concepts, 
and are not considered further. Fig. 20 shows that 
for a dose of 3 rem, there are minimum reactor radia- 
tion intensities that allow an exposure time of 91 
days (13 weeks) for orbital inclinations less than 71°. 
Intensities below these minimums give less than 3 
rem in 13 weeks when added to the cosmic-ray back- 
ground and are not of interest. The same is true 
(Fig. 21) for inclinations of less than 48° when con- 
sideration is given to a 5-rem yearly dose. Fig. 22 
shows that the allowed reactor radiations can reach 
a comparatively high level and still allow consider- 
able time in orbit before 25 rem are received. It is 
noted that in all cases the allowed exposure times 
increase only slightly for reactor intensities below 
10~4 rem/hr. 
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From Figs. 19 through 22 it appears that allowable 
reactor intensities will lie between 10! and 10> 
rem/hr for long-duration manned flights above the 
atmosphere, with reactor levels above 10~? rem/hr 
at the position of interest very nearly constituting the 
total radiation intensity, and levels below 10~4 
rem/hr allowing the cosmic radiation to make up 
most of the total intensity. 


Minimum Weight Reactor Shield Hand Calculation 


General Formulation 


Operation of a nuclear reactor above the at- 
mosphere will generally require that constraints be 
placed upon its radiation intensity at specified ex- 
ternal locations. These constraints may be neces- 
sary to prevent radiation damage to material com- 
ponents, to keep the biological dose within a recom- 
mended value, or to keep heat generation at various 
positions within a specified rate. The constraints 
are accomplished with shielding, and for space ap- 
plications, since mass is an important cost considera- 
tion, the shielding should be of minimum mass. 

From the theory of Lagrange multipliers, if it is 
desired that a function of variables, 9(x1, . . ., Xj, 
..., X,), be a maximum or minimum where the same 
variables satisfy 7 compatible constraint equations of 
the form: 


nit = (1) 


with k = 1, 2,...,7< m, then the proper x; must also 
satisfy : 


with 7 = 1, 2,...,. Equations 1 and 2 represent 
n + j equations in the variables x, and the j con- 
stants A,,, where the A,, are Lagrange multipliers. 

For a multilayer shield, components of which have 
thicknesses 1), to, . . ., tn, the external reactor radiation 
intensity at time 7 and location p can be expressed 
in the form: 


I,(r) to, ta, T) (3) 


Typical constraints upon the ?’s are that the total 
reactor radiation intensity at position p be a specified 
value, or the total dose be a certain quantity. The 
first type of constraint can be written: 


I,(r) = I, = K, = Constant (4) 


and is generally applicable to long-term nuclear 
auxiliary power sources. The second type of con- 
straint can be written: 


{ ‘Ip(1)dt = C, = Constant (5) 


and, in general, applies to a nuclear source of propul- 
sive power, which is a time-dependent radiation 
source. The initial reactor start-up time is taken to 
be 79, and the time of final shutdown 7;. It has been 
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Fig. 28. Dose attenuation factor for 4.0-mev gamma rays in lead: 
Plane monodirectional source. 


shown previously that for long-duration flights above 
the atmosphere, the allowable radiation intensity of 
an auxiliary nuclear reactor for a manned system 
will have a value between 10~! and 10° rem/hr in 
nearly all cases. The necessary reactor shield weight 
increases with no corresponding increase in the al- 
lowed flight time for reactor intensities below 10-5 
rem/hbr, and a 25-rem dose is received in just 10 days 
for a reactor intensity of 10-! rem/hr. A constraint 
of the type in Eq. 4 should therefore have K, such 
that it lies between 10~—! and10~ rem/hr for manned 
extra-atmospheric flights of duration longer than 10 
days, with position p denoting the region where the 
crew is located. Specific shielded reactor intensity 
requirements are determined for low-altitude orbits 
by choosing the desired total dose, orbital inclina- 
tion, and time in orbit. These variables are not in- 
dependent, since for each dose and orbital inclina- 
tion there is a corresponding definite flight time. 
For distances beyond 30,000 km from the surface of 
the earth, the desired dose and flight time alone de- 
termine the allowed reactor radiation intensity at the 
position of interest. 

The function that is to be minimized is the expres- 
sion for the shield mass, which can be expressed as: 


W = Wh, te, . . ., tn) (6) 


and the equations which must be solved are those of 
Egs. 4 or 5 and: 


W/dt, = d(OIp/dt,) (7) 
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where 


Radiation Attenuation by a Two-Region Shield With 
Axial Symmetry 


The method of radiation attenuation utilized is 
one where the spatial change in radiation intensity in 
a one-dimensional absorbing medium is considered 
to be proportional to the intensity, J, itself. The 
factor of proportionality is an arbitrary function of 
the spatial coordinate, x. In other words: 


dI = —IX(x)dx (9) 
and a particular solution is: 
t 
I = I, exp (10) 
0 
1.0 
8 
6 


> 


Fast Neutron Attenuation Factor 
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Fig. 29. Dose attenuation factor for fast neutrons, measured after 


traversal of a thick lithium hydride shield: Lead shielding thickness pre- 
cedes lithium hydride. 
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where ¢ is the distance in the one-dimensional absorb- 
ing medium where / is measured, and J = J) at ¢ = \). 
A special case might be one in which the absorber 
is distributed along the axis of a cylindrical reactor. 
Thus 2(«) depends on the energy of the incident radi- 
ation, and where one shielding substance follows an- 
other, the incident energy spectrum entering the sec- 
ond absorber is a function of the total preceding shield 
thickness. Near all material interfaces, the radiation 
intensity is also a function of the reflection coefficients 
of both bounding materials. Complete accuracy 
for a finite multilayer shield therefore requires = to be 
a function of all involved materials, and to be repre- 
sented as D(x, X2,..., X,). However, in most prac- 
tical considerations, if the intensity outside a thick 
reactor shield is formulated as a nonseparable func- 
tion x2, . . ., X,), it differs insignificantly from 
the product of separable functions, such as 
F\(x1)Fo(x2) . . Fr(x,), with each F, dependent 
only on a single shield material. The nomenclature 
to be used here is that of the separable function 
product. With this condition, any attenuation fac- 
tors by definition can be cast into the desired form 
of Eq. 10. The attenuation of reactor primary 
gamma radiation by the gamma shield of thickness 
7) can be represented by the factor: 


Ti 
Any = exp — (11) 


the attenuation of primary gamma radiation by the 
neutron shield of thickness 72 can be represented by 
the factor: 


T2 
yw Li(te)dty (12) 


the attenuation of fast neutrons by the gamma shield 
can be represented by the factor: 


Ti 
Ay = exp: — (13) 


and the attenuation of fast neutrons by the neutron 
shield can be represented by the factor: 


Te 
Ax Lo(te)dte (14) 
0 


When secondary gammas produced upon neutron 
absorption and inelastic collision in the shield are 
neglected, the total reactor radiation intensity, /,, at 
the desired position from the reactor is given by: 


IT, = + = K, (15) 
= + (16) 
= + (17) 


I, is the unshielded primary gamma ray intensity at 
the desired position of measurement, J; is the un- 
shielded fast neutron intensity at the desired posi- 
tion of measurement, and K, is the constant of con- 
straint, as in Eq. 4, of value between 107! and 10~° 
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rem/hr. Only direct-beam radiations are con- 
sidered, and no structure is assumed to be located in 
a position to make scattering a problem. 

The choice of basic shield geometry for an actual 
design is somewhat arbitrary and depends on items 
such as the number of coolant ducts necessary to re- 
move heat which is generated in the shield through 
the radiation absorption process. For idealized 
geometries with reactor, shielded region, and shield 
all having axial symmetry, weight functions for 
configurations like those of Fig. 23 can be shown to 
be: 


Spherical Segment Shield 


W = (2486/3) {Tip [3C? + + 
+ Tope[8(C + Ty)? + 
3T(C + T:) + (18) 
= + Ti)? + 
pol + 27, + T2)] (19) 
OW/OT2 = + 7, + Tr)? (20) 


Truncated Conical Segment Shield 


W = {Tip [(C’ + + 
CC"(2C’ + C"T:)] + 
Top2[(C’ + + C"T2)? + + Ti) X 
(2C’ + 2C"T; + C"T»)]} (21) 


OW/0T, = (#/3) { m[(C’2+ 2CC'C") + 
+ 2C’) + 
3C"*7 + 2[(2CC”? + 
4C'C") + 3C°(2T; + T:)]} (22) 


OW/OT> = (mp2/3) + 2CC'C” + 
2C"(CC” + 2C’)(T; + T2) + 
3C"(T, + Tr)?} (23) 
Disc Shield 
W = + poT2) (24) 
OW/d7T, = (25) 
OW/OT2 = (26) 


In the foregoing, 


C =rlat+h)/(rp thtk 

C= (ht —7)/ath+n 

C” = (rp — + h) 

ele h+ w/e + 

’- = core radius of the axially symmetric reactor 

ry = radius of the axially symmetric crew com- 
partment, which is located at the desired 
position where the reactor intensity is 


measured 

a = distance along the reactor axis from the 
near reactor assembly face to the crew 
compartment 


h = length of the reactor assembly 

k = distance from the reactor assembly to the 
shield 

pi = gamma shield density 

p2 = neutron shield density 


Lead-Lithium Hydride Shield for Specific 300-Kw Electrical 
Reactor 


Choice of shield materials is somewhat arbitrary, 
but lead and lithium hydride are nearly as efficient 
shielding substances for gamma rays and fast neu- 
trons as any that have been experimentally tested. 
However, various heavy metals have been tried in 
the nuclear-propelled aircraft (ANP) program in the 
search for a gamma shield, and a material which may 
prove to be an even better neutron shield than 
lithium hydride, which has a density of 0.87 gm/cm', 
is lithium borohydride,*' of absolute density 0.66 
gm/cm', and bulk density 0.2 to 0.25 gm/cm*. The 
shielding properties of both lead and lithium hydride 
are well known, and so final selection here is of these 
substances. 

The percent abundances® of the isotopes of lith- 
ium, as it occurs in nature, are 7.5 for lithium-6 and 
92.5 for lithium-7. The neutron activation cross 
sections for gamma emission are 28 and 33 milli- 
fermis* for lithium-6 and lithium-7, respectively, and 
the gamma emission half-life upon neutron absorp- 
tion in lithium-7 is 0.85 sec. The thermal neutron 
absorption cross section®*—** for lithium-6 is about 
945 fermis with no gamma emission, while the ab- 
sorption cross section of hydrogen is 0.33 fermis, and 
there is gamma radiation produced in the absorption 
process. It is therefore seen that the secondary 
gamma source in lithium hydride arises mainly from 
neutron captures in the hydrogen, with the subse- 
quent average emission per capture of one photon of 
energy between 2 and 3 mev. This is a relatively 
minor reaction compared with neutron absorption 
in the compound in regions where the primary reactor 
gamma radiation is not extremely small compared to 
the neutron flux. The compound itself has a fast 
neutron removal cross section of 0.338 in.—!, while 
the cross section of water is 0.256 in.~!. The neutron 
attenuation of lithium hydride can therefore be esti- 
mated on the basis that 3 in. of it are as effective as 
are 4 in. of water, whose attenuation characteristics 
are directly measurable. This being so, the esti- 
mated fast neutron dose attenuation factor for an 
electrical power reactor is presented in Fig. 24 as a 
function of lithium hydride thickness along a 28-in. 
diameter disc fission source® axis. Gamma-ray at- 
tenuation in lithium hydride is found by again first 
calculating the attenuation in water. The gamma 
leakage energy spectrum emitted by the desired reac- 
tor can be decomposed into energy groups. The 
contribution of each group to the dose at any posi- 


* One fermi equals 10724 cm?. 
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Fig. 30. Macroscopic absorption coefficients for fast neutrons and 
4.0-mev gamma rays in lead. 


tion in the shield can be found by using both the dose 
buildup factor for the appropriate energy, and the 
macroscopic cross section for the same energy. Use 
of the plane monodirectional buildup factor is justi- 
fied if the source is not too close to the shield. Other- 
wise, a surface source integration must be performed 
to obtain good accuracy. However, an average pri- 
mary gamma shield energy of 4 mev is often in good 
agreement with experimental data. For water, 
then, the macroscopic total cross section® for this 
energy is 0.0863 in.—!, while the plane monodirec- 
tional dose buildup factor, B,(t), is plotted in Fig. 25 
as a function of water thickness, ¢. The product of 
B,(t) and exp-0.0863¢ gives the gamma dose attenua- 
tion factor for water, and since gamma attenuation 
is proportional to the electron density of the material 
of traversal, 1 in. of water is as effective as are 1.22 
in. of lithium hydride. The dose attenuation factor 
for lithium hydride is given in Fig. 26. 

The percent abundances of the isotopes of lead, as 
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it occurs in nature, are 1.5, 23.6, 22.6, and 52.3 for 
lead-204, 206, 207, and 208, respectively. The neu- 
tron activation cross sections of lead-204 and 208 are 
0.7 fermis and 0.6 millifermis, respectively, with half- 
lives of 5 X 10’ yr and 3.2 hr. The thermal neutron 
absorption cross sections are 0.8, 0.025, 0.70, and less 
than 0.030 fermis for isotopes 204, 206, 207, and 208, 
respectively, which means that the cross section of 
the naturally occurring element is 0.17 fermis. Over 
90 percent of the gamma radiation emitted in the ab- 
sorption process has energy 7.4 mev. Neutron in- 
elastic scattering in lead is an additional source of 
gamma radiation. In this process, at least part of 
the lost neutron energy goes to excite the lead nu- 
cleus, which in turn loses its excitation energy 
through the emission” of one or more gamma rays. 
The neutron threshold energies for this interaction 
are ().90, 0.80, and 2.62 mev for the isotopes 204, 206, 
and 208, respectively. The emitted gamma energies 
extend from 3.0 to 0.35 mev in discrete levels for in- 
cident neutron energies less than 4 mev, while they 
are essentially continuous for neutron energies above 
4 mev. The total inelastic cross section for lead is 
quite low for neutron energies less than 1 mev, is 
about 1 fermi at 3 mev, and increases almost pro- 
portionally with energy up to 7 mev. This shows 
that neutron absorption in lead is quite small, and 
production of secondary gamma radiation can be 
significant only if it occurs in a region where the 
primary gamma radiation is low in intensity com- 
pared to the neutron flux. The primary gamma ray 
attenuation in lead is found by once again noting 
that an average energy of 4 mev in the radiation 
shield is in good agreement with experimental data. 
The macroscopic total cross section then, for this 
energy, is 1.19 in.~'. The plane monodirectional 
dose buildup factor, B2(t;), is plotted in Fig. 27 
as a function of lead thickness, ¢,, and the product 
B,(t;) exp-1.194,, which is the gamma dose attenua- 
tion factor for lead, is given in Fig. 28. The fast 
neutron removal cross section concept is of value in 
dealing with the shielding effect of lead when it is 
followed by a substantial thickness of lithium hy- 
dride, as it indeed is for a neutron radiation in- 
tensity reduction by a factor of at least 10’. It is 
not currently known at what material thickness the 
removal cross section ceases to be a constant, but it 
is assumed here to be unchanging for a thickness up 
to 18in. The commonly accepted macroscopic value 
is 0.294 in.—'. Exp-0.294t, is plotted in Fig. 29, and 
represents the change in neutron intensity outside a 
thick lithium hydride shield upon adding a lead 
thickness f; to the reactor side of the lithium hydride. 
The gamma and neutron macroscopic absorption 
cross sections are obtained from Figs. 24, 26, 28, 
and 29, and are plotted in Figs. 30 and 31. 

The specific reactor to be shielded, as an example, 
is a nuclear thermionic unit® which operates at a 
thermal power of 2 Mw, produces 300-kw electric 
power through the use of cesium vapor thermionic 
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converters, and has a lithium-7 cooled core which is a 
cylinder of length 11 in. and radius 5.25 in., and is 
surrounded on all sides by 3-in. thick beryllium re- 
flectors. The total weight of the reactor system, in- 
cluding the lithium radiator, is 2,300 lb. The un- 
shielded radiations 30 ft from the core center along 
its axis, where, in general, the greatest intensities 
occur, are 1.6 X 10° rem/hr for fast neutrons and 
2.4 X 104 rem/hr for gamma rays.*® The extent to 
which the previously given attenuation factors and 
macroscopic cross sections are now applicable is not 
known, unless the exact spectra of neutrons and 
gamma rays are specified for the unshielded reactor 
system. It is important that the attenuation factors 
utilized be indeed reasonable for the particular reac- 
tor to be shielded, especially when only thin shields 
are needed. The macroscopic cross sections in Fig. 
31, for instance, have very different curvatures. 
The neutron cross section is higher at small lithium 
hydride thicknesses than it is at large ones, while the 
opposite is true for the gamma cross section. How- 
ever, the neutron source is a fission disc, and many 
neutrons are of low-energy and easily absorbed in the 
lithium hydride next to this source. The gamma 
rays, on the other hand, are from an idealized mono- 
energetic source. Buildup occurs because there is 
scattering as well as absorption within the shield. 
As will be seen, though, the shields in the example are 
thick, so it is assumed the attenuation factors and 
cross sections already given are applicable. It has 
also been noted that regions where secondary gamma 
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Fig. 31. Macroscopic absorption coefficients for fast neutrons and 
4.0-mev gamma rays in lithium hydride. 
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Fig. 32. Optimized reactor shield weights as a function of reactor 
radiation intensity at the position of interest for various separation 
distances. 


production are important are where the primary 
gamma intensity has a low value and the neutron 
flux at the same time is substantial. Such a situation 
is expected to occur where a large lead shield pre- 
cedes the neutron absorber, and this occurs in this 
illustrative example in nearly all cases. However, 
adequate treatment requires volume integration of 
the shield secondary source radiation contribution 
terms. This is not now done and is considered to be 
a shortcoming upon which little improvement can 
be made by hand computation. 

The crew region is now chosen to have a 15 ft 
diameter. Therefore, 7, = 90 in., r, = 5.25 in., 
h = 14in., a = 360 in. when J; = 2.4 X 104 rem/hr, 
and J; = 1.6 X 10° rem/hr. p; = lead density = 
0.408 Ib/in.*, pp = lithium hydride density = 0.0314 
Ib/in.’, Ay is given in Fig. 28, As is given in Fig. 26, 
Ag is given in Fig. 29, As. is given in Fig. 24, and 
2.(7;) are given in Figs. 30 and 31. Spherical seg- 
ment shield geometry is utilized. For the case when 
there are only two shield variables, 7; and 72, and 
only one constraint equation (Eq. 15), Eq. 7 becomes: 


= (27) 
Substitution of Eqs. 16, 17, 19, and 20 into Eq. 27 
gives 


[ex(C + 7; + T2)?] = 
[21(72)AuAw + X 
+ 71)? + peT2(2C + 27, + T2)] (28) 
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Eq. 28 is solved graphically by choosing 7) and 
seeing what value of 72 makes both sides equal. 
When the left- and right-hand sides of Eq. 28 are 
equal, 7; and 72 give both a minimum weight shield 
and the radiation intensity of Eq. 15. Optimized 
reactor shield weights are presented in Fig. 32 for the 
minimum value of the distance, k, from the reactor 
assembly to the shield. These weights can be seen 
to vary from 4,100 lb for the worst case considered, 
namely J, = 10~° rem/hr at 30 ft, to 500 Ib for the 
best case, where 7, = 10~!rem/hr at 200 ft. Further 
weight minimization of the shielded reactor system 
requires that the weight of the necessary supporting 
structure for various separation distances be added 
to the system weights, and the sum minimized. 
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the Sun on May 12, 1959, Phys. Review Letters, Vol. 3, 
No. 4, p. 183, Aug. 15, 1959. 
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19 Brown, R. R., and D’Arcy, R. G., Observations of Solar 
Flare Radiation at High Latitude during the Period July 10-17, 
1959, Phys. Review Letters, Vol. 3, No. 8, p. 391, Oct. 15, 1959, 

2 Winckler, John R., Non-Relativistic Protons from Solar 
Flares, presented at the Symposium on Upper Atmosphere 
Solar Relations, Chicago, IIl., Dec. 28, 1959. 

21 Ehmert, A., Erbe, H., Pfotzer, G., Anger, C. D., and 
Brown, R. R., Observations of Solar Flare Radiation and Modu- 
lation Effects at Balloon Altitudes, July 1959, Jour. of Geo- 
physical Research, Vol. 65, No. 9, p. 2685, Sept. 1960. 

22 Anderson, K. A., and Enemark, D. C., Observations of Solar 
Cosmic Rays near the North Magnetic Pole, Jour. of Geophysical 
Research, Vol. 65, No. 9, p. 2657, Sept. 1960. 

*3 Anderson, K. A., Arnoldy, R., Hoffman, R., Peterson, L., 
and Winckler, J. R., Observations of Low-Energy Solar Cosmic 
Rays from the Flare of 22 August 1958, Jour. of Geophysical 
Research, Vol. 64, No. 9, p. 1133, Sept. 1959. 

24 Winckler, J. R., and Bhavsar, P. D., Low-Energy Solar 
Cosmic Rays and the Geomagnetic Storm of May 12, 1959, Jour. 
of Geophysical Research, Vol. 65, No. 9, p. 2637, Sept. 1960. 

*% Robey, D. H., Radiation Shield Requirements for Two 
Large Solar Flares, Astronautica Acta, Vol. VI, No. 4, p. 206, 
1960. 

26 National Bureau of Standards, Unusual Solar-Terrestrial 
Event, 1GY Bulletin, No. 40, p. 6, Oct. 1960. 

*7 Carmichael, H., Steljes, J. F., Rose, D. C., and Wilson, 
B. G., Cosmic-Ray Neutron Increase from a Flare on the Far 
Side of the Sun, Phys. Review Letters, Vol. 6, No. 2, p. 49, 
Jan. 15, 1961. 

*8 Covington, Arthur E., and Harvey, Gladys A., 10.7-cm 
Solar Noise Burst of November 20, 1960, Phys. Review Letters, 
Vol. 6, No. 2, p. 51, Jan. 15, 1961. 

29 NBS, Central Radio Propagation Laboratory, Compila- 
tions of Solar-Geophysical Data, Boulder, Colo., 1960. 

30 Kiepenheuer, K. O., Solar Activity, The Sun, I, The Solar 
System, p. 332, ed. by G. P. Kuiper, Univ. of Chicago Press, 
1953. 

31 Pettit, Edison, The Sun and Stellar Radiation, Astro- 
physics, p. 284, ed. by J. A. Hynek, McGraw-Hill Book Co., 
New York, N.Y., 1951. 

82 Ql’, A. I., The Possible Periodicity of Solar Activity, Astro- 
nomicheskii Zhurnal, Vol. 37, No. 2, p. 222, March 1960. 

33Van Allen, James A., The Nature and Intensity of the 
Cosmic Radiation, Physics and Medicine of the Upper At- 
mosphere (book), pp. 254-55, Univ. of N. Mex. Press, 1952. 

* Haber, Heinz, Physical Factors of the Space Environment, 
Space Technology (book), p. 27-16; John Wiley & Sons, 
Inc., New York, N.Y., 1959. 

8 Singer, S. F., The Primary Cosmic Radiation and Its Time 
Variations, Progress in Elementary Particle and Cosmic Ray 
Physics, Vol. IV, p. 288; North Holland Publishing Co., 
Amsterdam, 1958. 

86 Minzner, R., Champion, K., and Pond, H., The AFCRC 
Model Atmosphere, AFCRC-TR-59-267, pp. 77-84, Air Force 
Cambridge Research Center, Aug. 1959. 

37 Parker, E. N., Cosmic Ray Modulation by Solar Wind, 
Phys. Review, Vol. 110, No. 6, p. 1445, June 1958. 

88 Elliott, H., Cosmic Ray Intensity in Interplanetary Space, 
Nature, Vol. 186, No. 4721, p. 299, April 23, 1960. 

39 Fan, C. Y., Meyer, Peter, and Simpson, J. A., Expert- 
ments on the 11- Year Changes of Cosmic-Ray Intensity using a 
Space Probe, Phys. Review Letters, Vol. 5, No. 6, p. 272, Sept. 
15, 1960. 

40 Langham, Wright H., Implication of Space Radiations in 
Manned Space Flight, Aerospace Medicine, p. 410, June 1959. 

‘1 Schaefer, Hermann J., Theory of Pretection of Man in the 
Region of the Primary Cosmic Radiation, Jour. of Aviation 
Medicine, p. 338, Aug. 1954. 

42 Schaefer, H. J., Protection of Humans from Heavy Nuclei 
of Cosmic Radiation in Regions Outside the Atmosphere, Astro- 
nautica Acta, Vol. I, p. 100, 1955. 

(Continued on page 96) 
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Letters 


Hodograph Techniques 


g We are delighted to see that the use 
of the hodograph for trajectory analysis 
is finally being practiced, as evidenced 
by the recent article ‘“The Hodograph 
and Ballistic Missile Trajectory Prob- 
Jems” by Erwin W. Paul, appearing in 
the February 1961 issue of AEROSPACE 
ENGINEERING. 

Your readers may be interested in our 
developments of the hodograph tech- 
nique over the past few years, as fol- 
low: 


(1) Altman, S. P., Orbital Transfer for 
Satellites, presented in Aug. 1959 at 
Fourth Ballistic Missile and Space Tech- 
nology Symposium; Pergamon Press, 
1961. 

(2) Pistiner, J. S., Some Character- 
istics of the Planar Satellite Orbit, ARS 
Journal, Vol. 30, No. 3, March 1960. 

(3) Altman, S. P., and Pistiner, J. S., 
Hodograph Analysis of the Orbital Trans- 
fer Problem for Coplanar, Non-Aligned 
Elliptical Orbits, ARS Preprint No. 
1447-60 presented in Dec. 1960 at 15th 
Annual Meeting of ARS. 


Additional developments of the hodo- 
graph for analysis of orbital transfer and 
rendezvous will shortly be presented at 
forthcoming technical meetings. 


Finally, please note that Dr. Paul’s 
parameters p, d are most simply ob- 
tained from the total energy equation for 
the missile. When the total velocity V 
is expressed in radial and normal ve- 
locity components and the radial dis- 
tance r is expressed as a function of the 
normal velocity component by virtue of 
the constant angular momentum H, 
the total energy equation becomes the 
equation of a circle defined by the hodo- 
graph parameters p, d. 


Samuel P. Altman, MIAS 
Operations Research Department 


Josef S. Pistiner 
Research Department 
The Martin Co. 
Denver, Colo. 


The Editors welcome letters 
from readers, although none can 
be acknowledged. All must be 
signed, but identities will be 
withheld on request. 


Who is this man? 


First, you should know a few things about him: 
He’s responsible, as a man who leads others 
through new frontiers must be; he’s a specialist 
. .. but a specialist with time for creative reverie; 
he welcomes new challenges and grows in learn- 
ing and stature with whatever he faces; he’s 
mature, dedicated, and inquisitive—traits of a 
true man of science. Who is he? He’s the indis- 
pensable human element in the operations of 
one of the Navy's laboratories in California. 
Could he be you? 


U. S. NAVAL ORDNANCE TEST STATION at China 
Lake and Pasadena: Research, development, test- 
ing, and evaluation of missiles, advanced propulsion 
systems, and torpedoes and other undersea 
weapons. 


U. S. NAVAL ORDNANCE LABORATORY at Corona: 
Development of guidance and telemetry systems and 
missile components. Research in IR spectroscopy, 
magnetism and semiconductors, etc. 


U. S. NAVAL RADIOLOGICAL DEFENSE LABORA- 
TORY at San Francisco: One of the nation's major 
research centers on nuclear effects and counter- 
measures. 


U. S. NAVY ELECTRONICS LABORATORY at San 
Diego: One of the Navy's largest organizations en- 
gaged in the research and development of radar, 
sonar, radio, and acoustics. 


PACIFIC MISSILE RANGE and U. S. NAVAL MISSILE 
CENTER at Point Mugu: National launching and in- 
strumentation complex, guided missile test and 
evaluation; astronautics; satellite and space vehicle 
research and development. 


U. S. NAVAL CIVIL ENGINEERING LABORATORY 
at Port Hueneme: Research, development, and 
evaluation of processes, materials, equipment, and 
structures necessary to the design, construction, and 
maintenance of the Navy's shore bases. 


Openings for Aeronautical Engineers, Chemists, 
Civil Engineers, Electronic Engineers, Electronic En- 
gineers (Digital Circuitry & Electro-Acoustic), Mathe- 
maticians (Test Data Processing & Analysis), 
Mechanical Engineers, Operations Research Ana- 
lysts, Physicists. 


The man we want must have an advanced de- 
gree, or a Bachelor’s degree with at least three 
years’ solid experience. He should contact... 
Personnel Coordinator, Dept.E 
U. S. Naval Laboratories in California 
1030 East Green Street 
Pasadena, California 


U. S. NAVAL LABORATORIES 


IN CALIFORNIA 
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DYNAMICS ENGINEER 


WE ARE SEEKING FOR OUR STAFF... 
An experienced senior research engineer to 
study dynamics problems in aeronautics, 
astronautics, and high-performance Naval 
vessels. 


A UNIQUE OPPORTUNITY... 
Diversified and stimulating research rang- 
ing from hydrofoil boats to spacecraft. 
Excellent opportunity for advancement. 
Starting salary open. 

IF YOU HAVE... 
An advanced degree plus experience in air- 
plane dynamics and vibrations, aeroelas- 
ticity, or stability and control, please send 
resume and salary requirements to: 


R. C. Mays, Director of Personnel 


Southwest Research Institute 


8500 Culebra Road 
San Antonio 6, Texas 


48 Young, O. B., and Chen, H. Y., Heavy Primary Cosmic 
Rays at Geomagnetic Latitute of 41° N, Phys. Review, Vol. 115, 
No. 6, p. 1719, Sept. 15, 1959. 


44 Yoshida, S., Ludwig, G. H., and Van Allen, J. A., Dis- 
tribution of Trapped Radiation in the Geomagnetic Field, Jour. 
of Geophysical Research, Vol. 65, No. 3, p. 807, March 1960. 

45 Vernov, S. N., Address to the General Meeting of the 
Academy of Sciences of the USSR, 1959. 

46 Schaefer, Hermann J., Radiation Danger in Space, Astro- 
nautics, p. 36, July 1960. 

47 Schaefer, Hermann J., Further Evaluation of Tissue-Depth 
Doses in Proton Radiation Fields in Space, U.S. Naval School 
of Aviation Medicine Rpt. 17, MR 005.13-1002, May 24, 
1960. 

48 Freden, Stanley C., and White, R. Stephen, Particle Fluxes 
in the Inner Radiation Belt, Jour. of Geophysical Research, 
Vol. 65, No. 5, p. 1877, May 1960. 

49 Fan, C. Y., Meyer, Peter, and Simpson, J. A., Trapped 
and Cosmic Radiation Measurements from Explorer VI, 
EFINS-59-73, Enrico Fermi Institute of Nuclear Studies, 
Chicago, IIl., Dec. 1959. 

50 Krassovsky, V. I., Results of Scientific Investigations Made 
by Soviet Sputniks and Cosmic Rockets, Astronautica Acta, Vol. 
VI, No. 1, p. 37, 1960. 

51 Callery Chemical Co., Lithium Borohydride, Tech. Bulletin 
C-130, Pittsburgh, Pa., Nov. 1, 1958. 

52Hughes, D. J., and Schwartz, R. B., Neutron Cross Sec- 
tions, BNL-325, 2nd Ed., Brookhaven National Laboratory, 
July 1, 1958. 

53 Mittelman, P. S., and Liedtke, R. A., Gamma Rays from 
Thermal Neutron Capture, Nucleonics, May 1955. 

54 Goldstein, H., Fundamental Aspects of Reactor Shielding, 
Addison-Wesley, 1959. 

55 McCool, W. J., and Otis, D. R., Calibration of ORNL-Lid 
Tank Fission Plates, Nucleonics, p. 98, April 1960. 

56 Rockwell, T., I1I (Ed.) Reactor Shielding Design Manual, 
p. 447, D. Van Nostrand Co., Inc., Princeton, N.J., 1956. 

57 Weinberg, A. M., and Wigner, E. P., The Physical Theory 
of Neutron Chain Reactors, Univ. of Chicago Press, 1958. 

88 The Martin Company-Nuclear Div., Proposal for 300-Kwe 
Nuclear Power System, MND-2285, p. III-1, Feb. 1960. 

59 Strzelecki, R., Personal Communication, The Martin 
Company-Nuclear Div., April 1960. 


SPECIAL IAS PUBLICATIONS 


Member Nonmember 

Aeronautical Engineering Index - 1956 and 1957 (each)... 10.00 15.00* 
Aeronautical Engineering Index - 1947 through 1955 (each)... 3.00 5 .00* 
National Midwestern Meeting on Air Logistics Proceedings - 1960 (116 pp.)............... 5.00 10.00* 
Weapons Systems Go - National Midwestern Meeting 

National Symposium on St Techniques Proceedings - 1960 (176 pp.)............ 5.00 10.00* 
Recovery of Space Vehicles Symposium Proceedings - 1960 (104 pp.)..................-. 4.00 8.00* 
Manned Space Stations Symposium Proceedings - 1960 (320 pp.)...................000055 5.00 10.00* 
National Specialists Meeting on Guidance of Aerospace Vehicles 

National Specialists Meeting on Dynamics & Acroelasticity Proceedings - 1958 (136 pp.)..... 3.50 6.00* 
1960 National Telemetering Conference Report....... 4.50 6.50* 
Frontiers of Science & Engineering Symposium - 1959 39 5.00 10.00* 
Fourth Turbine-Powered Air Transportation Meeting Proceedings - 1959.................... 2.00 4.00* 
Second Turbine-Powered Air Transportation Meeting Proceedings - 1955................... 3.50 6.00* 
Tth Anglo-American Aeronautical Conference Proceedings - 1959. 20.00 25 
5th International Aeronautical Conference Proceedings - 1955... 15.00 20.00* 
2nd International Aeronautical Conference Proceedings - 5.00 5.00* 
National Naval Aviation Meeting Proceedings - 1957 (124 pp.)..... ee 3.00 5 .00* 
Index to Books on Selected Technical Subjects in the [AS Library 


* Add $1.00 for orders outside the U.S.A. 


These may be obtained by writing to: Special Publications Dept., IAS, 2 E. 64th St., N.Y. 21, N.Y. 
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Personnel 


This section is for the use of individual members of the Institute 
seeking new connections and eligible organizations offering em- 
ployment to specialists in the aerospace industry. Any member 
or eligible organization may have requirements listed without 
charge by writing to the Secretary of the Institute. 


Wanted 


Graduate Assistants—Staff Members— 
University of Kansas, Dept. of Aeronau- 
tical Engineering, Lawrence, Kan., has 
openings for Graduate Assistants and two 
full-time staff members between now and 
September 1961, to teach courses in aero- 
dynamics (basic, theoretical, compressible 
flow, and hypersonics); mechanics of flight 
(including stability and control, aircraft 
design, and aircraft dynamics); aircraft 
structures; wind tunnel and other aero- 
nautical laboratories. Salary and academic 
rank will depend upon education and experi- 
ence. Write Ammon S. Andes, Chairman, 
giving training, teaching interest, starting 
rank and salary desired. Teaching load will 
be reduced when staff engages in research. 
Plans are complete for new building, money 
is partly appropriated. 


Structures or Vibration Engineer—For 
position in research program involving 
analysis of turbojet and turboprop engine 
vibration and engine-induced vibration in 
airframes. Ph.D. and knowledge of com- 
puter programing highly desirable. Oppor- 
tunities for teaching in area of advanced 
structures. Salary $14,000 to $17,000 de- 
pending upon education and experience. 
Send résumé to Chairman, School of Aero- 
nautical & Space Engineering, University of 
Oklahoma, Norman, Okla. 


Technical Personnel—The following posi- 
tions are now open for the TAPCO Group, 
Thompson Ramo Wooldridge, Inc., 23555 
Euclid Ave., Cleveland 17, Ohio: Space 
Power Development Engineer, Turbo- 
machine Preliminary Designer, Advanced 
Designer, Heat-Transfer Specialist, Ad- 
vanced Propulsion Engineer, Missile Propul- 
sion Specialist, Structural Supervisor, Fluid 
Flow Specialist, Thermodynamicist—Heat 
Exchangers, Senior Aerothermodynamicist, 
Auxiliary Power Specialist, Propulsion Sys- 
tems Preliminary Analyst, and Weapons 
Systems Analyst. Write fully to R. J. 
Theibert, Employment Manager. 


Engineers & Scientists—Arde Associates 
has moved into its new facilities at 100 
Century Rd., Paramus, N.J. Engineers and 
scientists of superior abilities are invited to 
investigate these permanent, rewarding 
Positions. Research Scientists (Ph.D. 
level)—Physicists or physical chemists with 
backgrounds in plasmas and _ solid-state 
Physics. Project Engineers—For direction 
of complete projects from supervision of de- 
sign to customer liaison in aeronautical com- 
ponent field. Product Design Engineers— 
Heavily experienced in environmental control 
system design, hydraulic systems and com- 
Ponents, aircraft control components, heat 
exchangers or valves. Mechanical Design 
Engineers—Minimum 10 years’ experience. 
Should be fully capable of carrying project 
from concept to successful completion. 
Stress Analysts—Minimum 7 years’ aero- 


The number preceding the notice 
represents the Box Number of the 
Institute of the Aerospace Sciences to 


which inquiries should be addressed. 


nautical experience in plates, shells, and 
thermal stresses. To arrange a confidential 
interview, write to E. A. Greenlee, 75 Austin 
St., Newark, N.J. 


Available 


111. Aeronautical Engineer—Age 21. Ex- 
perienced in statistical analysis, in-process 
data systems, and reliability of electronic 
components. Position desired as design 
engineer, reliability engineer, manufacturing 
engineer, liaison engineer, etc. Résumé 
sent upon request. 


110. Aerothermo Consulting Engineer— 
West Coast Area. Available as a private 
consultant and/or engineering services 
specialist for problems on accessory turbo- 
machinery .elements and systems. Expe- 
rienced in design point and off-design per- 
formance prediction analysis, fluid dynamic 
design, test programing and test data 
analysis of turbines, compressors, fans and 
pumps. Analysis and design, including 
thermochemical, heat-transfer and cycle 
analysis, of power machinery units and 
systems for space, airborne, and stationary 
applications. Related experience in dynamic 
analysis of mechanical systems, analog 
solutions of problems, mechanical design of 
rotating elements and systems and analysis, 
and design of electrical and electronic systems 
and controls. Age 37; M.S. and BS., 
Licensed Professional Engineer. 


109. Chairman, Aeronautical (or Me- 
chanical) Engineering Department—Dr. de- 
gree, M.E. (Aero Option). Diversified back- 
ground in aeronautical and mechanical 
engineering education combined with indus- 
trial experience in charge of responsible 
research and development activities. Has 
planned and directed sponsored R&D effort 
of considerable scope with outstanding 
results. Demonstrated ability to plan and 
direct research and stimulate desire for 
research by others. P.E.; Member IAS, 
ASME, and other societies. Has authored 
many reports, papers, patents, etc. 


108. Vice-President or Manager, Product 
Development (or Engineering)—Age 48. 
Advanced M.E. degree (Aero option). 
Extensive experience as Director of engineer- 
ing and development in mechanical and 
structural field (airframe industry, trans- 
portation industry and bearing industry). 
Familiar with markets, manufacturing, man- 
power and money problems, besides engineer- 


ing. Has had to plan and carry out pro- 
grams and activities which have yielded 
tangible results in higher sales and profits. 
Would consider assignment with progressive 
organization in aircraft, transportation equip- 
ment, etc. 


107. Manager, Graduate Engineer, & 
Pilot—with over 10 years’ successful experi- 
ence in weapon and automation systems 
R&D desires managerial position in small, 
technically oriented company. Available in 
June. Résumé upon request. 


106. Senior Test Engineer—B.S. in 
A.E.; age 31. Nine years’ experience in 
aerospace fields including design, flight test, 
and assignment as airborne launch engineer. 
Five years as flight-test supervisor and aero- 
mechanical supervisor participating in more 
than one hundred missile launchings of 
various types. Broad experience covering 
rockets, turbojets, ram-jets, aerodynamics, 
structures, electronics, and parachute sys- 
tems. Desires supervisory or staff position 
with a new and growing test organization 
located in the western states away from big 
cities. Résumé on request. 


105. Aeronautical Engineer—B.A.E., M.- 
S.E.; age 37. Seventeen years’ experience 
in missile, airplane, helicopter, and VTOL 
aerodynamics. Eight years in development, 
four as project aerodynamicist; and nine in 
applied research and teaching in charge of 
prime sponsored programs. Desires re- 
sponsible technical direction position in 
industrial or university research. Complete 
résumé on request. 


104. Management Consultant—Special- 
ist in systems management, analysis, and 
organization. Extensive background in cor- 
poration long-range planning and diversifi- 
cation. Academic, industrial, and military 
background encompasses M.S.A.E., M.B.A. 
candidate, USAF weapon system project 
officer, and 10 years’ experience in aerospace 
industry principally in the management of 
missile and space systems. Associate Fellow 
IAS, Senior Member IRE. Desires part- 
time consulting activities. Fee schedule 
and résumé upon request. 


103. Senior Mechanical Engineer—B.S. 
and M.S. in M.E., 17 years’ experience in 
design, specifications, with strong emphasis 
on heat transfer, thermal expansion for 
reliability. Heat exchanger calculation; fore- 
casting schedules, administering funds, etc., 
on a project basis. Registered Professional 
Engineer in Massachusetts. Résumé on 
request. 


102. Aeronautical Engineer—B.S._ in 
M.E., M.S. in Aeronautics; age 30. Three 
years’ experience in research, system analysis 
for space application. Familiar with foreign 
markets and languages. Desires challenging 
position in international department, over- 
seas assignment preferred. Résumé upon 
request. 
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A word 
about 
our 
Advertisers 


The many companies whose 
advertising appears in this issue 
are interested in you and your 
work. They are anxious to let 
you know what they are doing 


to help you solve research and 


design problems in missile, 
rocket, air, and_ spacecraft 
projects. 

AEROSPACE ENGINEERING’S 


advertising pages keep you 
posted on new and improved 
materials, components, services, 
and supplies useful to your 


professional work. 


To request more information 
on any product or service ad- 
vertised, may we suggest you 
write to the advertiser directly, 
at no obligation to you. It 
would be greatly appreciated if 
you would mention that you 


saw the ad in AEROSPACE 


ENGINEERING. 


ENGINEERING 
Established 1934 
2 E. 64th St., New York 21,N. Y. 
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the 1960 AERO/SPACE ENGINEERING 
CATALOG 


The only publication devoted exclusively to the aerospace 
industry, this CATALOG serves as a valuable buyers’ and 
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and missile parts and materials. It is distributed annually 
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components; Government Agencies; Research Organi- 
zations; etc. 
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From General Motors...'Temperatures Made to Order 


HARRISON GIVES 


ACOOL 


Where heat is a critical problem, you'll find 
Harrison on the job! Harrison heat exchangers 
dependably and efficiently control vital operating 
temperatures on the most advanced military j 
and civilian aircraft— piston, turbo-prop and 
turbo-jet powered. That’s why Kaman BA 
specified Harrison engine and transmission oil 
coolers for its H-43B Huskie turbine helicopter. 7 
Powered by a Lycoming T-53 engine, the | ii 
Huskie passed the most rigid reliability tests |e 
with flying colors. You'll find this same built-in foo 
reliability in 15 basic types of Harrison heat ? “A 
transfer construction—a complete line of designs 
which permits the selectivity that assures the 
right heat exchanger for every application. 

To save time and money on your temperature 
control problems, be sure to call in a Harrison 
Sales Engineer at the design stage. 


Harrison heat exchangers—quality pro 
ucts of General Motors Research 


A 1) 


RES 70 ORDER 


AIRCRAFT, AUTOMOTIVE, MARINE AND INDUSTRIAL HEAT EXCHANGERS 


For an informative 48-page 


brochure on the complete Harrison 
line . . . write to Department 902 


HARRISON RADIATOR DIVISION, GENERAL MOTORS CORPORATION, LOCKPORT, NEW YORK 


P 
| 


total competence in computation and data processing —the breadth, the brains and the backgroU 


What it means: total competence in computation as it applies to the contract team program and airborne, space, surla 
and underwater systems. What it stands for: total capability in all areas from basic research through field service, bat 
by 75 years devoted to computation and data processing. What it's been doing: directing the ALRI team and develo) 
its miniaturized airborne data processors; developing and producing Atlas guidance computers and high-speed compl 
for Polaris; developing high-speed computers for Mauler, U. S. Army’s newest automatic firing air defense syst 
Where it’s going: to the undefined, the unexplored, the unknown... serving the common good and the commong 


Burroughs Corporation i 
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